10 rem soltn to ex3 - least squares - j.h.stacey keynes 12 let t = 0 13 let t = t+1 20 input n 40 if n < 0 then 135 42 let s1,s2,s3,s4 = 0 45 dim m(19,1) 50 for i = 0 to n-1 55 input m(i,0) 65 input m(i,1) 75 let s1 = s1+m(i,1) 80 let s2 = s2+m(i,0) 85 let s3 = s3+m(i,0)*m(i,1) 90 let s4 = s4+m(i,0)^2 92 next i 95 let a = (s4*s1-s2*s3)/(n*s4-s2^2) 100 let b = (n*s3-s1*s2)/(n*s4-s2^2) 102 print 104 print "- - - - - - - - - - - - - - - - - - - -" 106 print 108 print "for y=a+b*x data block ";t;" gives a=";a,"b=";b 110 print 115 print "i","x(i)","y(i)","a+b*x(i)" 120 for i = 0 to n-1 125 print i,m(i,0),m(i,1),a+b*m(i,0) 127 next i 130 goto 13 135 print "no more coordinates" 140 end