Bug #10574 fixed - Runge-Kutta-Method failed for vector [x, 1] with x > 9.
[scilab.git] / scilab / modules / differential_equations / src / fortran / lsrgk.f
1 c     ====================================
2 c     ode Gateway for Adaptative fourth order Runge Kutta 
3 C
4 c     The original version has been modified to pass the  rwork
5 c     argument to odeint
6 c     array + blas use. Serge Steer INRIA- feb 2012
7
8 c     ====================================
9       subroutine lsrgk (f, neq, y, t, tout, itol, rtol, atol, itask,
10      1            istate, iopt, rwork, lrw, iwork, liw, jac, mf)
11      
12       external f, jac,rkqc
13       integer neq, itol, itask, istate, iopt, lrw, iwork, liw, mf
14       double precision y, t, tout, rtol, atol, rwork
15       integer nok,nbad
16       dimension neq(*), y(*), rtol(*), atol(*), rwork(lrw), iwork(liw)
17       integer iero
18       common/ierode/iero
19       iero=0
20       call odeint(y, neq, t,tout,atol(1),1.0d-4,0.0d0,nok,nbad,f,rkqc,
21      $     rwork)
22       t=tout
23       if (iero.gt.0) istate=-1
24       return
25       end
26 c     ====================================