share/contrib/gentran/test/runge.template

   1 C fortran code generated from runge.template
   2       program runge
   3       implicit real (k,m)
   4 c
   5 c  Input
   6 c
   7       write(6,*) 'Initial Value of p'
   8       read(5,*) p
   9       write(6,*) ' p = ', p
  10       write(6,*) 'Initial Value of q'
  11       read(5,*) q
  12       write(6,*) ' q = ', q
  13       write(6,*) 'Value of m'
  14       read(5,*) m
  15       write(6,*) ' m = ', m
  16       write(6,*) 'Value of k0'
  17       read(5,*) k0
  18       write(6,*) ' k0 = ', k0
  19       write(6,*) 'Value of b'
  20       read(5,*) b
  21       write(6,*) ' b = ', b
  22       write(6,*) 'Step Size of t'
  23       read(5,*) hh
  24       write(6,*) ' Step Size of t = ', hh
  25       write(6,*) 'Final Value of t'
  26       read(5,*) tp
  27       write(6,*) ' Final Value of t = ', tp
  28 c
  29 c  Initialization
  30 c
  31       tt=0.0
  32 <<
  33       gentran( literal(tab, "write(9,*) ' h = ", eval(h), "'", cr),
  34                 literal(tab, "write(9,*) ' d = ", eval(d), "'", cr) )$
  35 >>
  36       write(9,901) c
  37 901   format(' c= ',e20.10)
  38       write(9,910) tt, q, p
  39 910   format(' ',3e20.10)
  40 c
  41 c  Loop
  42 c
  43 <<
  44       gentran( unless tt >= tf do
  45                (
  46                      rsetq(pn, ev(pn,expand)),
  47                      rsetq(q, ev(qn,expand)),
  48                      p : pn,
  49                      tt : tt + hh,
  50                      literal(tab, "write(9,910) tt, q, p", cr)
  51                ) )$
  52 >>
  53       stop
  54       end