perm filename ALPHA.BAS[NET,GUE] blob sn#041851 filedate 1973-05-11 generic text, type T, neo UTF8
00001  dim a(10),b(10)
00005 for i=1 to 6
00010 read a(i),b(i)
00015 next i
00050 r1=fnr(8.78)
00060 r2=fnr(6.06)
00070 print r1,r2
00075 z=.1
00080 for i=0 to r1 step z
00090 if i>r2 then 102
00100 v2=fne(r2-i)
00101 goto 110
00102 v2=0
00110 v1=fne(r1-i)
00120 v=.66*v1/8.78 + .34*v2/6.06
00125 if i+.5>r1 then 135
00130 w1=fne(r1-i)-fne(r1-i-.5)
00131 goto 140
00135 w1=0.
00140 if i+.5>r2 then 155
00150 w2=fne(r2-i)-fne(r2-i-.5)
00151 goto 160
00155 w2=0
00160 w=.66*w1+.34*w2
00170 print i,v,w
00180 next i
00200 stop
00500 def fnr(e)=((((a(6)*e+a(5))*e+a(4))*e+a(3))*e+a(2))*e+a(1)
00600 def fne(r)=((((b(6)*r+b(5))*r+b(4))*r+b(3))*r+b(2))*r+b(1)
00800 data 6.01399354E-2,-1.79547046E-1
00810 data 4.53491928E-1,2.6205927
00820 data -3.49173014E-2,-5.67273163E-1
00830 data 2.98647499E-2,9.33053908E-2
00840 data -3.33483157E-3,-7.96080248E-3
00850 data 1.30241831E-4,2.67955391E-4
00999 end