5 DEFDBL A-H, O-Z 10 D = .0001: INPUT "A,B = "; A, B: N = 10 20 N = N + 10: H = (B - A) / N: X = A: S = 0!: T = 0! 30 FOR I = 1 TO N STEP 2 40 GOSUB 100: S = S + Y * H / 3: T = T + Y * H / 2: X = X + H 50 GOSUB 100: S = S + 4 * Y * H / 3: T = T + Y * H: X = X + H 60 GOSUB 100: S = S + Y * H / 3: T = T + Y * H / 2: NEXT I 70 IF ABS((S - T) / (S + T)) > D THEN 20 80 PRINT "Integral = "; T, N 90 PAUSE: INPUT "Error % = "; D: GOTO 20 100 Y = SQR(1! - SIN(EXP(1) / 2) ^ 2 * SIN(X) ^ 2): RETURN