<html>

  <head>

    <meta content="text/html; charset=ISO-8859-1"

      http-equiv="Content-Type">

  </head>

  <body bgcolor="#FFFFFF" text="#330000">

    Hi<br>

    <br>

    Could you show us the actual loop? <br>

    <br>

    <br>

    On 03/10/2011 17:12, grivet wrote:

    <blockquote cite="mid:4E89D0DB.8060405@cnrs-orleans.fr" type="cite">Hello,

      <br>

      I wish to solve a system of linear differential equations of first

      order

      <br>

          x' = Rx

      <br>

      with R a symmetric matrix of coefficients. The formal solution is

      simple:

      <br>

          x = S*exp(Kt)*inv(S)*x(0)

      <br>

      where S is the matrix that diagonalizes R:

      <br>

          K = inv(S)*R*S

      <br>

      With Scilab, I would do, for instance,  t = 0:0.01:10, but then,

      what is the best way to

      <br>

      compute exp(Kt)  and the above matrix product ?

      <br>

      Up to now, I have worked with  2*2 matrices R and K, so that I

      compute exp[K(1,1)*t(i)], exp[K(2,2)*t(i) and x(t(i)) within

      <br>

      a loop, but that doesnot seem very efficient.

      <br>

      Thanks in advance

      <br>

      JP Grivet

      <br>

      <br>

    </blockquote>

    <br>

    <br>

    <div class="moz-signature">-- <br>

      Adrien Vogt-Schilb (Cired) <br>

      Tel: (+33) 1 43 94 <b>73 77</b></div>

  </body>

</html>