<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
On 16/10/2012 15:55, LE BELGUET Alix 226276 wrote:
<blockquote
cite="mid:277FFF3261DC4F48941A2444664D02AF17DFE626@EXDAG0-B1.intra.cea.fr"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=ISO-8859-1">
<meta name="Generator" content="Microsoft Word 14 (filtered
medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif][if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="WordSection1">
<p class="MsoNormal">Bonjour,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Je souhaite utiliser la fonction ‘impl’ de
Scilab afin de résoude un système d’équations différentielles
linéaire du type : A(t,y) dy/dt=g(t,y), y(t0)=y0.<o:p></o:p></p>
<p class="MsoNormal">Cependant, je n’arrive pas à comprendre ce
que signifie l’argument ‘adda’ définie par : r=A(t,y)+p.
<o:p></o:p></p>
<p class="MsoNormal">Comment choisir p ?<o:p></o:p></p>
</div>
</blockquote>
<br>
Bonjour Alix,<br>
<br>
Par exemple, vous cherchez a resoudre en dimension 2:<br>
Id * y' = y avec y(0) = [1;2]<br>
<br>
function r=res(t,y,ydot)<br>
r = y-eye(2,2)*ydot;<br>
endfunction<br>
<br>
function r=adda(t,y,p)<br>
r=eye(2,2)+p<br>
endfunction<br>
<br>
puis:<br>
t=0:0.1:1;<br>
y=impl([1;2],[0;0],0,t,res,adda)<br>
<br>
On verifie que y vaut a peu pres [1;2]*exp(t)<br>
max(abs(y-[1;2]*exp(t)))<br>
<br>
retourne chez moi 1.3e-6.<br>
<br>
Dans le code ci-dessus, il suffit de remplacer eye(2,2) par A(t,y)
et dans res de remplacer y par g(t,y). Vous n'avez pas a vous
souciez de ce que vaut p...<br>
<br>
Calixte<br>
<br>
<blockquote
cite="mid:277FFF3261DC4F48941A2444664D02AF17DFE626@EXDAG0-B1.intra.cea.fr"
type="cite">
<div class="WordSection1">
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Je ne parviens pas non plus à trouver
d’exemple parlant de cette fonction, c’est pourquoi je me
tourne vers vous…<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">En vous remerciant par avance pour votre
réponse,<o:p></o:p></p>
<p class="MsoNormal">Cordialement,<o:p></o:p></p>
<p class="MsoNormal">Alix LB<o:p></o:p></p>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
users-fr mailing list
<a class="moz-txt-link-abbreviated" href="mailto:users-fr@lists.scilab.org">users-fr@lists.scilab.org</a>
<a class="moz-txt-link-freetext" href="http://lists.scilab.org/mailman/listinfo/users-fr">http://lists.scilab.org/mailman/listinfo/users-fr</a>
</pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">--
Calixte Denizet
Software Development Engineer
-----------------------------------------------------------
Scilab Enterprises
143bis rue Yves Le Coz - 78000 Versailles, France
<a class="moz-txt-link-freetext" href="http://www.scilab-enterprises.com">http://www.scilab-enterprises.com</a></pre>
</body>
</html>