10 CLS
20 SCREEN 1
30 WINDOW (-5, 15)-(30, -30)
40 FOR X = -5 TO 5 STEP .02
50 Y = (X ^ 2) + X + 3
60 PSET (X, Y)
70 NEXT X
80 LINE (-5, 0)-(7, 0)
90 LINE (-3, 9)-(-3, 0)
100 LINE (0, 3)-(0, 0)
110 LINE (3, 15)-(3, 0)
120 LOCATE 7, 3
130 PRINT "Y0"
140 LOCATE 8, 11
150 PRINT "Y2"
160 LOCATE 8, 7
170 PRINT "Y1"
180 LOCATE 9, 8
190 PRINT "h"
200 LOCATE 9, 4
210 PRINT "-h"
220 LOCATE 10, 1
230 PRINT "  "; CHR$(244); "h"
240 PRINT "A="; CHR$(245); "-h (CX^2+DX+C)=[CX^3/3+DX^2/2+EX"
250 PRINT CHR$(179); " h"
260 PRINT CHR$(179); "-h = h/3(2Ch^2+6E)"
270 PRINT "Y0=Ch^2-Dh+E"
280 PRINT "Y1=E"
290 PRINT "Y2=Ch^2+Dh+E"
300 PRINT "Y0+4Y1+Y2=2Ch^2+6E"
310 PRINT "AND A=(h/3)*(Y0+4Y1+Y2)"
320 PRINT CHR$(244); "1"
330 PRINT CHR$(245); "0  F(1+X) with 4 intervals is"
340 PRINT "((1/4)/3)*[F(0)+4F(1/4)+2F(2/4)+4F(3/4)+F(1)]"
350 PRINT "1+5+3+7+2=18 *((1/4)/3)=1.5"
360 LOCATE 2, 15
370 PRINT "SIMPSON'S METHOD"
380 LOCATE 3, 15
390 PRINT "OF INTEGRATION"
400 PLAY "mbo4fgagaagaa#a#"
410 PLAY "mbaa#a#ad#d#dd#p4"
420 PLAY "mbaagaagddcdfo5cc"
430 PLAY "mbo4a#o5cco4a#"
440 PLAY "mbo5ddddco4"
450 FOR UU = 1 TO 10000000
460 NEXT UU
470 SYSTEM