bonjour j' ai un problème dans le programme matlab pour résoudre une équation différentielles par la méthode d EULER

voila mon programme est ce que vous pouvez résoudre le problème.
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
function dx1=f1(t,x1)
dx1=x2(t)
end
function dx2=f2(t,x2)
dx2=50*(x2(t)^-1)
end
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
h=0.1;
t=0:h:0.7;
n=0.7/h;
x1(1)=0;
x2(1)=-1;
for i=1:n;
    x1(i+1)=x1(i)+h.*f1(t(i),x1(i))
    x2(i+1)=x2(i)+h.*f2(t(i),x2(i))
end
merci