[users-fr] Recherche equivalent ODE45 de matlab

[Serge Steer]

Je me suis trompé l'option qui permet de retourner les valeurs de l'etat 
a chaque mesh point fonctionne bien
%ODEOPTIONS=[2,0,0,%inf,0,2,500,12,5,0,-1,-1];

// dy/dt=y^2-y sin(t)+cos(t), y(0)=0
function ydot=f(t, y),ydot=y^2-y*sin(t)+cos(t),endfunction
y0=0;t0=0;t=%pi;

y=ode(y0,t0,t,f)
plot(y(1,:),y(2,:))

Serge Steer
Le 10/01/2012 13:08, GOBILLOT Gilles 116128 a écrit :
> Bonjour,
> je suis à la recherche d'un équivalent de la fonction matlab "ode45", 
> sans succès.
> La fonction scilab "ode" nécessite de connaitre le vecteur temps, que 
> je n'ai pas (c'est en fait une variable inconnue à cette étape qui 
> tiens lieu de temps).
> En matlab j'écrirais [y,t]=ode45(t0, y0, f); alors que scilab me 
> propose y=ode(y0,t0,t,f)
> Merci d'avance pour votre aide.
> Gilles Gobillot.
> 	
>
> PAvant d'imprimer, pensez à l'environnement
>

-------------- section suivante --------------
Une pièce jointe HTML a été nettoyée...
URL: <https://lists.scilab.org/pipermail/users-fr/attachments/20120110/5a6b6d3f/attachment.htm>