Euler implicite

Une âme charitable pourrait m'envoyer le programme d'Euler implicite? ou l'afficher ?

Parceque j'ai toujours le même message d'erreur depuis ce matin et je commence a ne plus avoir de cheveux....

Si ça peut vous motiver ... il y aura des photos de botanique exclusives à la clef.

Voila :

Code:

clear
a=0; b=20; h=0.5; M=10;           
N=(b-a)/h;
x(1)=a;
y(1)=a;
g=inline('z-y-h*(1-(z/(2+sin(x))))');
j=inline('1+h*(1/(2+sin(x)))');
for k=2:N
    z(1)=y(k-1);
    for m=1:M
        z(m+1)=z(m)-(g(h,x(k-1),y(k-1),z(m))/j(h,x(k-1)));
    end
    y(k)=z(m+1);
    x(k)=x(k-1)+h;
end
plot(x,y);

Par contre je peux pas te l'expliquer mais ca marche