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