équations différentielles sous matlab

[inviteccb0e881]

Bonjour à tous;

Je n'arrive pas à résoudre un système d'équations différentielles sous matlab :

(1-y)*dK/dy=59.66-1.57*K

y appartient à [0,1]

je cherche à tracer K en fonction de y par itération
l'initiation de k est k(0)=30

Merci d'avance
-----

[invite551c2897]

Bonjour.
Tu devrais mettre ton post sur le forum des developpeurs - matlab !
Voici une solution :

clear
kprime=inline('(59.66-1.57*k)/(1-y)');
[y,k]=ode45(kprime,[0,1],30);
plot(y,k)
grid

[inviteccb0e881]

Merci de m'avoir répondu

je cherche svp un code par itération

Merci d'avance

[invite551c2897]

Par Runge-Kutta d'ordre deux :

function [x,y]=kutta2(x0,y0,xf,h,f)
% Résolution d'équation différentielle par la méthode
% de Runge-Kutta d'ordre 2
% [x,y]	 : tableau de valeurs de la fonction y(x)
% x0, y0 : conditions initiales
% h 	 : pas 
% xf 	 : valeur finale de x
% f	 : expression de la fonction y'=f(x,y)
% f est une chaîne de caractères 
 
tx=x0:h:xf;
ty=y0;
for i=1:floor((xf-x0)/h)
 x=tx(i);y=ty(i);
 f1=eval(f);
 x=tx(i)+h/2; y=ty(i)+f1*h/2;
 f2=eval(f);
 ty(i+1)=ty(i)+h*f2;
end 
x=tx;y=ty;

Programme Principal :

clear
x0=0; y0=30;% conditions initiales
xf=1;% valeur finale
h=0.1;% pas
f='(59.66-1.57*y)/(1-x)';
[a,b]=kutta2(x0,y0,xf,h,f);
plot(a,b);%,symb(i));
hold on
grid

[A voir en vidéo sur Futura]