From e0c6872cf40896c7be36b11dcc744620f10adf1d Mon Sep 17 00:00:00 2001 From: Norbert Preining Date: Mon, 2 Sep 2019 13:46:59 +0900 Subject: Initial commit --- info/maad/pgm6.bas | 11 +++++++++++ 1 file changed, 11 insertions(+) create mode 100644 info/maad/pgm6.bas (limited to 'info/maad/pgm6.bas') diff --git a/info/maad/pgm6.bas b/info/maad/pgm6.bas new file mode 100644 index 0000000000..50effbc92e --- /dev/null +++ b/info/maad/pgm6.bas @@ -0,0 +1,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 + -- cgit v1.2.3