<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Bonjour,<br>
<br>
Quelques remarques :<br>
<br>
- Si l'idée est effectivement de résoudre l'équation différentielle,
je n'ai pas entièrement vérifié le script mais<br>
j'ai tendance à faire confiance au moins à la méthode employée par
l'auteur <br>
- Si je lis de plus près, "l'idée est de générer une suite de
fonctions", il me semble que l'objectif ne soit pas entièrement
ceci. Est-ce vrai ?<br>
- point mineur, mais peut-être pas puisqu'il intervient dans une
e.d., La fonction S n'apparaît pas très bien "conditionnée". J'aurai
plutôt écrit S comme la limite de s(x,y,t)<br>
quand y tend vers 1 par exemple,<br>
- Sinon, je suis beaucoup plus sensible à l'esthétique de phi et H
(pour moi ça compte dans un modèle)<br>
- Je suppose que vous avez exploité la symétrie en r1+r2 ? Si on le
fait, on trouve une équation en log (phi) en fonction de S et H.
Mais je n'ai pas trouvé d'autre combinaison simple de r1 et r2
permettant de transformer le problème.<br>
- Pour moi, générer une suite de fonctions, peut consister à
développer r1 et r2 en série de Taylor par rapport à la variable t,
et injecter les coefficients dans les équations différentielles. On
obtient alors une suite de fonctions dont les termes pour r1(n) et
r2(n), sont les approximations à l'ordre n de chaque série. Par
contre, il faut faire au préalable un développement de S à l'ordre 1
ou 2 par rapport à la variable t, en injectant les développements de
Taylor de r1 et r2 dans S.<br>
<br>
JLuc<br>
<br>
<br>
<br>
<div class="moz-cite-prefix">Le 26/08/2013 21:35, Serge Steer a
écrit :<br>
</div>
<blockquote cite="mid:521BAE16.7090500@scilab.org" type="cite">
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
<div class="moz-cite-prefix"><br>
Je ne suis pas sur de comprendre ce que vous voulez faire, mais
s'il s'agit d'integrer le systeme d'eqations differentielles<br>
vous pouver le faire avec le code Scilab suivant<br>
<br>
function v=f(y,r,D,Beta,x)<br>
v=exp(-(x-y)^2/D^2)*r(1)/(r(1)^Beta+r(2)^Beta)<br>
endfunction<br>
<br>
function v=dr(y,r,alpha,D,Beta,x)<br>
v= alpha*intg(0,1,list(f,r,D,Beta,x))<br>
endfunction<br>
<br>
r0=[1;1];<br>
y0=1<br>
Beta=2;<br>
x=0;<br>
alpha=[1;2];<br>
D=1.5;<br>
r1=ode(r0,0,0:0.01:1,list(dr,alpha,D,Beta,x))<br>
<br>
<br>
Serge Steer<br>
INRIA<br>
Le 24/08/2013 23:26, <a moz-do-not-send="true"
class="moz-txt-link-abbreviated"
href="mailto:pro.olive@ouvaton.org">pro.olive@ouvaton.org</a>
a écrit :<br>
</div>
<blockquote
cite="mid:172b04ec2d53278fcc3fe7d0afc99a47@ouvaton.org"
type="cite">Bonjour a tous, <br>
<br>
Je mets en fichier joint le problème que je rencontre au sujet
de Scilab : l'idée est de pouvoir générer une suite de
fonctions. Je ne sais donc pas si Scilab est le meilleur
logiciel qui soit... <br>
<br>
Merci. <br>
Marc Olive<br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
users-fr mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:users-fr@lists.scilab.org">users-fr@lists.scilab.org</a>
<a moz-do-not-send="true" 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>
<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>
<div class="moz-signature">-- <br>
<meta http-equiv="content-type" content="text/html;
charset=ISO-8859-1">
<title>repons</title>
<br>
<br>
<br>
<br>
<small>Jean-Luc BABIGEON <br>
CNRS/Lal/Dacc <br>
Laboratoire de L'accélérateur Linéaire <br>
Centre scientifique d'Orsay - Bat 200 - BP 34 <br>
91898 - Orsay Cedex <br>
tel 0164468247 ou 8427 -<a moz-do-not-send="true"
href="mailto:babigeon@lal.in2p3.fr" target="_blank">babigeon@lal.in2p3.fr</a>
</small><br>
<br>
</div>
</body>
</html>