Help pour un RK4 sous scilab en matriciel

[invite2ea9837a]

Voici mon code, il est buggé j'arrive pas à m'en sortir.

clf()
clear

// Variables
c=1;
n=1;
k=1;

Am1=zeros(2,2) ;
A=[1 0 ; c n ];
B=[0 -1 ; k 0 ];
Am = zeros (2,2) ;
Am = inv(A);
//************************************************************




function g = f(t,Y) // défini la fonction f de l'ennoncé.
z=1 // fonction z à gauche du signe = de l'equation diff
J=[0 ; z ];
g = Am *(J - B*Y)
endfunction


////////***************************Conditions initiales********************************///////////

y0=0;
yp0=0;
a=0;
b=1.55;

function G = F(t,Y,h) // défini la fonction Phi

l1=f(t,Y) ;
l2=f(t+h/2 , [1 h/2 ; 0 1 ] * Y + [ h*l1/8 ; l1/2 ] );
l3=f(t+h/2 , [1 h/2 ; 0 1 ] * Y + [ l2*h/8 ; l2/2 ] );
l4=f(t+h , [1 h ; 0 1 ] * Y + [ h/2*l3 ; l3 ] );
G = [h h ;0 1]*[(Y(2,1));(h*(l1+2*l2+2*l3+l4))/6] ;

endfunction

function Dessine // dessine la solution

N= 100; // en combien on découpe
x= [0:N]/N*(b-a);
h= (b-a)/N;

o1(1) = [y0 ; yp0 ] ; // on initialise la matrice qui suit la suite géométrique

for i=2:N+1
o1(i)=o1(i-1)+F(h*(i-2),o1(i-1),h); // cette ligne beug ! ******************** ***************** ***************** //
end;
plot2d(x,o1,style=0)

endfunction

Dessine

et voici les erreurs de scilab :

-->endfunction
-->//
-->Dessine
!--error 15
Sous-matrice incorrectement définie.

at line 7 of function Dessine called by :
Dessine
at line 58 of exec file called by :
-3c55-4bc7-8424-b8245b63233e.sce', 1
while executing a callback



Je serais reconnaissant si vous pouviez m'aider.

Merci
-----