<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#ffffff" text="#330000">
    thank you that was it<br>
    <br>
    On 07/07/2011 20:54, Stephane Mottelet wrote:
    <blockquote
      cite="mid:20110707205435.15943l8mczvn47b4@webmail.utc.fr"
      type="cite">Hello,
      <br>
      <br>
      use the "root" parameter of ode (it is well explained in the help
      page of this macro).
      <br>
      <br>
      S.
      <br>
      <br>
      Adrien Vogt-Schilb <a class="moz-txt-link-rfc2396E" href="mailto:vogt@centre-cired.fr"><vogt@centre-cired.fr></a> a écrit :
      <br>
      <br>
      <blockquote type="cite">Hi everyone,
        <br>
        <br>
        When one whishes to solve an ODE y'=f(t,y) starting at point
        (t0,y0), one
        <br>
        can use the ode function with the following arguments ode(y0,t0,
        [t0:tf], f),
        <br>
        where [t0:tf] is the interval on which one wants to solve the
        ODE.
        <br>
        <br>
        The terminal condition of the resolution is, in the ode
        function, a time
        <br>
        condition (t==tf). Does anybody know if and how we can use
        instead a terminal condition on y?
        <br>
        <br>
        <br>
        Let us suppose, for example, that we know that the solution
        increases to
        <br>
        infinity (maybe in finite time). The problem is to solve the ODE
        as long
        <br>
        as y stays under a certain value yf. We do not know which
        temporal interval is to
        <br>
        be used to reach this yf.
        <br>
        <br>
        Numerical example: we want to solve y'=1+a.*y^3 from y0=0 to
        yf=10 for
        <br>
        differents values of a.
        <br>
        <br>
        Any ideas ?
        <br>
        <br>
        <br>
        Best regards,
        <br>
        Adrien
        <br>
        <br>
        AVS
        <br>
        <br>
        <br>
      </blockquote>
      <br>
      <br>
    </blockquote>
    <br>
    <br>
    <div class="moz-signature">-- <br>
      <p style="color: teal; font-size: 12pt;"><b>Adrien Vogt-Schilb</b><br>
        Research Fellow<br>
        <a class="moz-txt-link-abbreviated" href="mailto:vogt@centre-cired.fr">vogt@centre-cired.fr</a><br>
        Tel: (+33) 1 43 94 73 96<br>
        Fax: (+33) 1 43 94 73 70<br>
      </p>
      <table style="border: medium none;">
        <tbody>
          <tr>
            <td style="padding: 0pt;">
              <p><img src="cid:part1.09090307.03050805@centre-cired.fr"
                  alt="CIRED" width="54"></p>
            </td>
            <td style="padding: 0pt 7pt;" valign="middle">
              <p style="font-size: 10pt; color: teal;">45 bis, Av de la
                Belle Gabrielle<br>
                F-94736 Nogent-sur-Marne<br>
                <a href="http://www.centre-cired.fr/" style="color:
                  teal;">http://www.centre-cired.fr/</a><br>
              </p>
            </td>
          </tr>
        </tbody>
      </table>
    </div>
  </body>
</html>