summaryrefslogtreecommitdiff
path: root/info/maad/pgm7.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/pgm7.bas
Initial commit
Diffstat (limited to 'info/maad/pgm7.bas')
-rw-r--r--info/maad/pgm7.bas12
1 files changed, 12 insertions, 0 deletions
diff --git a/info/maad/pgm7.bas b/info/maad/pgm7.bas
new file mode 100644
index 0000000000..d2dc90bab2
--- /dev/null
+++ b/info/maad/pgm7.bas
@@ -0,0 +1,12 @@
+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
+