summaryrefslogtreecommitdiff
path: root/info/maad/pgm6.bas
diff options
context:
space:
mode:
authorNorbert Preining <norbert@preining.info>2019-09-02 13:46:59 +0900
committerNorbert Preining <norbert@preining.info>2019-09-02 13:46:59 +0900
commite0c6872cf40896c7be36b11dcc744620f10adf1d (patch)
tree60335e10d2f4354b0674ec22d7b53f0f8abee672 /info/maad/pgm6.bas
Initial commit
Diffstat (limited to 'info/maad/pgm6.bas')
-rw-r--r--info/maad/pgm6.bas11
1 files changed, 11 insertions, 0 deletions
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
+