f195bc16a108f39760aea15e816f6abcb014b990
[scilab.git] / scilab / modules / differential_equations / src / fortran / lsrgk.f
1 c     ====================================
2 c     runge kutta d'ordre 4 adaptatif 
3 c     l'interface lsrgk a ete fait en s'inspirant de lsode 
4 c     voir lsode.f pour comprendre le sens des variables 
5 c     ====================================
6       subroutine lsrgk (f, neq, y, t, tout, itol, rtol, atol, itask,
7      1            istate, iopt, rwork, lrw, iwork, liw, jac, mf)
8      
9       external f, jac,rkqc
10       integer neq, itol, itask, istate, iopt, lrw, iwork, liw, mf
11       double precision y, t, tout, rtol, atol, rwork
12       integer nok,nbad
13       dimension neq(*), y(*), rtol(*), atol(*), rwork(lrw), iwork(liw)
14       integer iero
15       common/ierode/iero
16       iero=0
17       call odeint(y,neq,t,tout,atol(1),1.0d-4,0.0d0,nok,nbad,f,rkqc)
18       t=tout
19       if (iero.gt.0) istate=-1
20       return
21       end
22 c     ====================================