PRINT
PRINT "f(x)=(1/2)*SIN(x) を シンプソン公式で積分"
PRINT
DEF f(x)=(1/2)*SIN(x)    ! f(x)=3*x^2 への置き換えも試みてください
LET a=0
LET b=PI                 ! b=1 ： f(x)=3*x^2 への置き換えたとき
LET n=10000
LET S=1
LET h=(b-a)/n            ! 小区間の幅
LET Q=0                  ! 単純な積分
LET U=h*(f(a)+f(a+n-1))
FOR i=1 TO n-2
   IF MOD(i,2)=0 THEN
      LET K=2
   ELSE
      LET K=4
   END IF
   LET Q=Q+(h/3)*K*f(a+i*h)
NEXT i
PRINT "  ";Q;"      誤差：";ABS(S-Q)
PRINT
END