Contenu | Rechercher | Menus

Annonce

Si vous avez des soucis pour rester connecté, déconnectez-vous puis reconnectez-vous depuis ce lien en cochant la case
Me connecter automatiquement lors de mes prochaines visites.

À propos de l'équipe du forum.

#1 Le 30/10/2006, à 16:02

du

probleme de programation sous scilab

Bonjour voici mon code scilab

function [t,X]=myeuler(f,t0,T,Xinit,N)
    h=T/N;
    t=t0+[0:N]*h;
    X=Xinit;
   
    for i=1:N 
    p=feval(t(i),X(:,i),'f');
      X(:,i+1)=X(:,i)+h*p;
    end
   
endfunction

-------------------------------
function YR=F(t,Y)
YR=Y
YR(2,1)=-g/L*sin(Y(1));

//YR(1)=Y(2,1);

endfunction
__________________________________
N=100;
T=50;
a=-1;
yo=1;
h=T/N;
t2=[0:N]*h;
t0=0;
g=10;
L=1;
t0=0;
getf("D:\FAC\TP_Scilab\F.sci");
getf("D:\FAC\TP_Scilab\myeuler.sci");
yEu=zeros(2,n+1);
yEu(1,1)=1; yEu(2,1)=0;
[t,yEu]=myeuler(F,t0,T,yEu,N);
plot2d(t,yEu(1,:),1)

et en fait mon probleme est avec la fonction F...le code ne marche pas...j utilise depuis peu scilab, j'ai fait pas mal de matlab...mais je ne vois pas pk ça ne marche pas.
Merci de m'aider.
   Emilio