SCREEN 7
PRINT "This is integration by Simpson's Method"
PRINT "change the four (1 - (g ^ 2)) ^ (1 / 2) equations"
PRINT "in the script to any equation"
PRINT "change the first a and b to any length of itegration"
PRINT "and c is the amount of interval accuracy."
PRINT "This simply calculates the area under a circle from 0 to 1"
PRINT "or in other words 1/4 of a 1 in. radius circle."
FOR zz = 1 TO 1000000
NEXT zz
DEFDBL A-Z
a = 0
b = 1
c = 511536
g = a
GOTO it
lo:
i = x
d = (b - a) / c
o = c / 2
re:
g = g + d
GOTO mu
it:
x = (1 - (g ^ 2)) ^ (1 / 2)
GOTO lo
rs:
i = i + x * 4
g = g + d
GOTO qw
mu:
x = (1 - (g ^ 2)) ^ (1 / 2)
GOTO rs
nc:
i = i + x * 2
o = o - 1
IF o > 0 THEN GOTO re
g = b
GOTO hs
ap:
i = i - x
w = d * i / 3
ba = b
PRINT "integral is"
PRINT w
PRINT "4 times integral is approximately pi"
PRINT 4 * w
GOTO enf
qw:
x = (1 - (g ^ 2)) ^ (1 / 2)
GOTO nc
hs:
x = (1 - (g ^ 2)) ^ (1 / 2)
GOTO ap
enf: