ode15i utilisation difficile?
Bonjour tous,
pourriez vous m'expliquez comment résoudre un systeme différentiel avec matlab svp ?
avec ode45 j'y arrive :
main.m
Code:
1 2 3 4
|
A=eye(3)+5;
B=1:1:3; B=B';
[t,y]=ode45(@(t,y)equadiffOde45(t,y,A,B), [0 100], zeros(1,length(B))); |
myFun.m
Code:
1 2 3 4
|
function dy=myFun(t,y,A,B)
dy=zeros(length(B),1);
dy(:,1)=A*y+B; |
mais je ne trouve pas l'equivalent avec ode15i
pourriez vous m'expliquez svp car j'ai essayé plein de chose et je n'y arrive pas :aie:
il y a deux chose que je ne comprends pas :
1) la forme sous laquelle doit etre ecrite myFun
2) l'appel à ode15i plus particulierement la dérivée initiale (ça devrait etre toujours 0 non?)
merci de votre aide :ccool: