円周率を求めるプログラムを区分求積法でやってみましたが･･･
10 x=0:r=1:s=0
20 input a
30 for b=1 to a
40 x=x+(r/a)
50 y=(r*r-x*x)^(1/2)
60 h=y*(r/a)
70 s=s+h
80 next b
90 p=4s
100 print p
110 end

50行で、エラーが出ます。何か間違いがあるのでしょうか?