1C fortran code generated from runge.template
2      program runge
3      implicit real (k,m)
4c
5c  Input
6c
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
28c
29c  Initialization
30c
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
37901   format(' c= ',e20.10)
38      write(9,910) tt, q, p
39910   format(' ',3e20.10)
40c
41c  Loop
42c
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
55