Euler implicite

Aller en bas

Euler implicite

Message  Philippe le Mer 5 Mar - 19:16

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. Rolling Eyes
avatar
Philippe
Mitochondrie
Mitochondrie

Nombre de messages : 69
Age : 31
Localisation : Bruxelles
Date d'inscription : 28/10/2007

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Re: Euler implicite

Message  Bender le Mer 5 Mar - 20:13

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 Very Happy
avatar
Bender
Enzyme
Enzyme

Nombre de messages : 18
Date d'inscription : 24/09/2007

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Revenir en haut

- Sujets similaires

 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum