<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Je me suis trompé l'option qui permet de retourner les valeurs de
l'etat a chaque mesh point fonctionne bien<br>
%ODEOPTIONS=[2,0,0,%inf,0,2,500,12,5,0,-1,-1];<br>
<br>
// dy/dt=y^2-y sin(t)+cos(t), y(0)=0<br>
function ydot=f(t, y),ydot=y^2-y*sin(t)+cos(t),endfunction<br>
y0=0;t0=0;t=%pi;<br>
<br>
y=ode(y0,t0,t,f)<br>
plot(y(1,:),y(2,:))<br>
<br>
Serge Steer<br>
Le 10/01/2012 13:08, GOBILLOT Gilles 116128 a écrit :
<blockquote
cite="mid:59863E6B3B5B7245A7B5660A00E8C89908F939@EXDAG0-A0.intra.cea.fr"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=ISO-8859-1">
<meta name="GENERATOR" content="MSHTML 8.00.6001.19170">
<div><span class="394340412-10012012"><font face="Arial" size="2">Bonjour,</font></span></div>
<div><span class="394340412-10012012"></span> </div>
<div><span class="394340412-10012012"><font face="Arial" size="2">je
suis à la recherche d'un équivalent de la fonction matlab
"ode45", sans succès.</font></span></div>
<div><span class="394340412-10012012"></span> </div>
<div><span class="394340412-10012012"><font face="Arial" size="2">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).</font></span></div>
<div><span class="394340412-10012012"></span> </div>
<div><span class="394340412-10012012"><font face="Arial" size="2">En
matlab j'écrirais [y,t]=ode45(t0, y0, f); alors que scilab
me propose <span class="id">y</span><span class="operator">=</span><span
class="command">ode</span><span class="openclose">(</span><span
class="id">y0</span><span class="default">,</span><span
class="id">t0</span><span class="default">,</span><span
class="id">t</span><span class="default">,</span><span
class="functionid">f</span><span class="openclose">)</span></font></span></div>
<div><span class="394340412-10012012"></span> </div>
<div><span class="394340412-10012012"><font face="Arial" size="2">Merci
d'avance pour votre aide.</font></span></div>
<div><span class="394340412-10012012"></span> </div>
<div><span class="394340412-10012012"><font face="Arial" size="2">Gilles
Gobillot.</font></span></div>
<div><span class="394340412-10012012"></span> </div>
<div class="Section1">
<table style="BORDER-COLLAPSE: collapse; mso-padding-alt: 0cm
0cm 0cm 0cm" class="MsoNormalTable" border="0" cellpadding="0"
cellspacing="0">
<tbody>
<tr style="PAGE-BREAK-INSIDE: avoid; mso-yfti-irow: 0;
mso-yfti-firstrow: yes; mso-yfti-lastrow: yes">
<td style="PADDING-BOTTOM: 0cm; PADDING-LEFT: 5.4pt;
WIDTH: 73.15pt; PADDING-RIGHT: 5.4pt; PADDING-TOP: 0cm"
width="98">
<span style="mso-fareast-font-family: 'MS Mincho';
mso-bidi-font-size: 10.0pt">
<!--?xml:namespace
prefix = o ns = "urn:schemas-microsoft-com:office:office" /-->
<o:p></o:p></span></td>
<td style="PADDING-BOTTOM: 0cm; PADDING-LEFT: 5.4pt;
WIDTH: 178.85pt; PADDING-RIGHT: 5.4pt; PADDING-TOP: 0cm"
width="238">
<br>
</td>
</tr>
</tbody>
</table>
<p style="mso-margin-top-alt: auto; mso-margin-bottom-alt: auto"
class="MsoNormal">
<span style="FONT-FAMILY: Webdings; COLOR: #339966; FONT-SIZE:
8pt">P</span><span style="COLOR: #339966; FONT-SIZE: 8pt">
Avant d'imprimer, pensez à l'environnement</span><o:p></o:p></p>
<p style="mso-margin-top-alt: auto; mso-margin-bottom-alt: auto"
class="MsoNormal">
<!--?XML:NAMESPACE PREFIX = U1
/-->
<u1:p><o:p> </o:p></u1:p></p>
</div>
<div> </div>
</blockquote>
<br>
</body>
</html>