summaryrefslogtreecommitdiff
path: root/info/maad/pgm6.bas
blob: 50effbc92e02bdf84437a24f906707ce581ffbb6 (plain)
1
2
3
4
5
6
7
8
9
10
11
10  DEF FNA (X) = SQR(1! - SIN(EXP(1) / 2) ^ 2 * SIN(X) ^ 2)
20  N = 10: A = 0!: B = 2! * ATN(1!)
30  H = (B - A) / N: S = 0!: T = 0!
40  FOR I = 1 TO N STEP 2
50  X = A + (I - 1) * H: S = S + FNA(X): T = T + FNA(X)
60  X = A + I * H: T = T + 4! * FNA(X)
70  X = A + (I + 1) * H: S = S + FNA(X): T = T + FNA(X)
80  NEXT I: S1 = H * S: S2 = H * T / 3!: E = ABS((S1 - S2) / (S1 + S2))
90  IF (E > .000001) THEN N = 2 * N: GOTO 30
100 PRINT "Integral = "; S2, N