#1 Le 30/10/2006, à 15: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