|
DataMuseum.dkPresents historical artifacts from the history of: RegneCentralen RC759 "Piccoline" |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about RegneCentralen RC759 "Piccoline" Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - download
Length: 5748 (0x1674)
Names: »FLYT.PAS«
└─⟦5a0460c8b⟧ Bits:30003931/GEM_Development-B.imd Disketter indleveret af Steffen Jensen (Piccolo/Piccoline)
└─⟦this⟧ »FLYT.PAS«
0x0000…0020 50 72 6f 67 72 61 6d 20 4e 6f 72 6d 61 6c 46 53 3b 0d 0a 0d 0a 7b 20 53 74 65 66 66 65 6e 20 42 ┆Program NormalFS; æ Steffen B┆
0x0020…0040 65 72 69 6e 67 20 4a 65 6e 73 65 6e 20 2f 20 50 72 6f 6a 65 6b 74 20 55 20 20 20 41 70 72 69 6c ┆ering Jensen / Projekt U April┆
0x0040…0060 20 31 39 39 30 20 7d 0d 0a 0d 0a 7b 20 54 61 62 65 6c 6f 70 73 6c 61 67 20 66 6f 72 20 69 6e 74 ┆ 1990 å æ Tabelopslag for int┆
0x0060…0080 65 67 72 61 6c 65 74 20 61 66 20 55 2d 6e 6f 72 6d 61 6c 66 6f 72 64 65 6c 69 6e 67 73 66 75 6e ┆egralet af U-normalfordelingsfun┆
0x0080…00a0 6b 74 69 6f 6e 65 6e 0d 0a 20 20 76 2e 68 2e 61 2e 20 53 69 6d 70 73 6f 6e 27 73 20 43 6f 6d 70 ┆ktionen v.h.a. Simpson's Comp┆
0x00a0…00c0 6f 73 69 74 65 20 61 6c 67 6f 72 69 74 6d 65 20 66 72 61 20 54 75 72 62 6f 20 50 61 73 63 61 6c ┆osite algoritme fra Turbo Pascal┆
0x00c0…00e0 73 20 4e 75 6d 65 72 69 63 61 6c 0d 0a 20 20 4d 65 74 68 6f 64 73 20 54 6f 6f 6c 62 6f 78 2e 20 ┆s Numerical Methods Toolbox. ┆
0x00e0…0100 7d 0d 0a 0d 0a 0d 0a 7b 24 49 2d 7d 20 20 20 20 20 20 7b 20 44 69 73 61 62 6c 65 20 49 2f 4f 20 ┆å æ$I-å æ Disable I/O ┆
0x0100…0120 65 72 72 6f 72 20 74 72 61 70 70 69 6e 67 20 7d 0d 0a 7b 24 52 2b 7d 20 20 20 20 20 20 7b 20 45 ┆error trapping å æ$R+å æ E┆
0x0120…0140 6e 61 62 6c 65 20 72 61 6e 67 65 20 63 68 65 63 6b 69 6e 67 20 7d 0d 0a 0d 0a 75 73 65 73 0d 0a ┆nable range checking å uses ┆
0x0140…0160 20 20 49 6e 74 65 67 72 61 74 2c 20 44 6f 73 2c 20 43 72 74 2c 20 43 6f 6d 6d 6f 6e 3b 0d 0a 0d ┆ Integrat, Dos, Crt, Common; ┆
0x0160…0180 0a 43 6f 6e 73 74 0d 0a 20 20 4e 75 6d 49 6e 74 65 72 76 61 6c 73 3a 49 6e 74 65 67 65 72 20 3d ┆ Const NumIntervals:Integer =┆
0x0180…01a0 20 33 30 30 20 3b 0d 0a 20 20 4c 6f 77 65 72 4c 69 6d 69 74 3a 46 6c 6f 61 74 20 20 20 20 20 3d ┆ 300 ; LowerLimit:Float =┆
0x01a0…01c0 20 2d 31 30 20 3b 0d 0a 0d 0a 56 61 72 0d 0a 20 20 55 70 70 65 72 4c 69 6d 69 74 20 20 3a 20 46 ┆ -10 ; Var UpperLimit : F┆
0x01c0…01e0 6c 6f 61 74 3b 20 20 20 20 20 20 20 20 20 20 20 20 20 20 7b 20 4c 69 6d 69 74 73 20 6f 66 20 69 ┆loat; æ Limits of i┆
0x01e0…0200 6e 74 65 67 72 61 74 69 6f 6e 20 7d 0d 0a 20 20 49 6e 74 65 67 72 61 6c 20 20 20 20 3a 20 46 6c ┆ntegration å Integral : Fl┆
0x0200…0220 6f 61 74 3b 20 20 20 20 20 20 20 20 20 20 20 20 20 20 7b 20 56 61 6c 75 65 20 6f 66 20 74 68 65 ┆oat; æ Value of the┆
0x0220…0240 20 69 6e 74 65 67 72 61 6c 20 7d 0d 0a 20 20 45 72 72 6f 72 20 20 20 20 20 20 20 3a 20 62 79 74 ┆ integral å Error : byt┆
0x0240…0260 65 3b 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 7b 20 46 6c 61 67 73 20 69 66 20 73 6f 6d 65 ┆e; æ Flags if some┆
0x0260…0280 74 68 69 6e 67 20 77 65 6e 74 20 77 72 6f 6e 67 20 7d 0d 0a 20 20 69 63 68 20 20 20 20 20 20 20 ┆thing went wrong å ich ┆
0x0280…02a0 20 20 3a 20 43 68 61 72 3b 0d 0a 0d 0a 7b 24 46 2b 7d 0d 0a 66 75 6e 63 74 69 6f 6e 20 70 68 69 ┆ : Char; æ$F+å function phi┆
0x02a0…02c0 28 58 20 3a 20 46 6c 6f 61 74 29 20 3a 20 46 6c 6f 61 74 3b 0d 0a 62 65 67 69 6e 0d 0a 20 20 70 ┆(X : Float) : Float; begin p┆
0x02c0…02e0 68 69 3a 3d 20 30 2e 33 39 38 39 34 32 32 2a 45 78 70 28 2d 30 2e 35 2a 58 2a 58 29 3b 0d 0a 65 ┆hi:= 0.3989422*Exp(-0.5*X*X); e┆
0x02e0…0300 6e 64 3b 0d 0a 7b 24 46 2d 7d 0d 0a 0d 0a 50 72 6f 63 65 64 75 72 65 20 47 65 74 4c 69 6d 69 74 ┆nd; æ$F-å Procedure GetLimit┆
0x0300…0320 28 76 61 72 20 55 70 70 65 72 4c 69 6d 69 74 20 3a 20 46 6c 6f 61 74 29 3b 0d 0a 42 65 67 69 6e ┆(var UpperLimit : Float); Begin┆
0x0320…0340 0d 0a 20 20 52 65 70 65 61 74 0d 0a 20 20 20 20 57 72 69 74 65 28 27 55 70 70 65 72 20 6c 69 6d ┆ Repeat Write('Upper lim┆
0x0340…0360 69 74 20 6f 66 20 69 6e 74 65 67 72 61 74 69 6f 6e 3f 20 27 29 3b 0d 0a 20 20 20 20 52 65 61 64 ┆it of integration? '); Read┆
0x0360…0380 6c 6e 28 55 70 70 65 72 4c 69 6d 69 74 29 3b 0d 0a 20 20 20 20 49 4f 43 68 65 63 6b 3b 0d 0a 20 ┆ln(UpperLimit); IOCheck; ┆
0x0380…03a0 20 55 6e 74 69 6c 20 28 6e 6f 74 20 49 4f 65 72 72 29 20 41 6e 64 20 28 55 70 70 65 72 4c 69 6d ┆ Until (not IOerr) And (UpperLim┆
0x03a0…03c0 69 74 3e 4c 6f 77 65 72 4c 69 6d 69 74 29 3b 0d 0a 45 6e 64 3b 0d 0a 0d 0a 50 72 6f 63 65 64 75 ┆it>LowerLimit); End; Procedu┆
0x03c0…03e0 72 65 20 52 65 73 75 6c 74 73 28 55 70 70 65 72 4c 69 6d 69 74 20 20 20 3a 20 46 6c 6f 61 74 3b ┆re Results(UpperLimit : Float;┆
0x03e0…0400 0d 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 49 6e 74 65 67 72 61 6c 20 20 20 20 ┆ Integral ┆
0x0400…0420 20 3a 20 46 6c 6f 61 74 3b 0d 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 45 72 72 ┆ : Float; Err┆
0x0420…0440 6f 72 20 20 20 20 20 20 20 20 3a 20 62 79 74 65 29 3b 0d 0a 62 65 67 69 6e 0d 0a 20 20 57 72 69 ┆or : byte); begin Wri┆
0x0440…0460 74 65 6c 6e 3b 0d 0a 20 20 57 72 69 74 65 6c 6e 28 27 55 70 70 65 72 20 4c 69 6d 69 74 3a 27 20 ┆teln; Writeln('Upper Limit:' ┆
0x0460…0480 3a 20 32 35 2c 20 55 70 70 65 72 4c 69 6d 69 74 20 3a 20 32 35 3a 35 29 3b 0d 0a 20 20 57 72 69 ┆: 25, UpperLimit : 25:5); Wri┆
0x0480…04a0 74 65 6c 6e 3b 0d 0a 20 20 42 65 67 69 6e 0d 0a 20 20 20 20 57 72 69 74 65 6c 6e 28 27 20 20 20 ┆teln; Begin Writeln(' ┆
0x04a0…04c0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 74 3d 27 2c 55 70 70 65 72 4c 69 6d 69 74 ┆ t=',UpperLimit┆
0x04c0…04e0 3a 38 3a 35 29 3b 0d 0a 20 20 20 20 57 72 69 74 65 6c 6e 28 27 20 20 20 20 20 20 20 20 20 20 20 ┆:8:5); Writeln(' ┆
0x04e0…0500 20 20 20 20 20 20 20 20 f4 20 20 27 29 3b 0d 0a 20 20 20 20 57 72 69 74 65 6c 6e 28 27 20 20 20 ┆ '); Writeln(' ┆
0x0500…0520 20 20 20 20 20 20 20 20 20 e8 28 74 29 20 3d 20 b3 20 70 68 69 28 78 29 20 64 78 20 20 3d 20 27 ┆ (t) = phi(x) dx = '┆
0x0520…0540 2c 20 49 6e 74 65 67 72 61 6c 20 3a 20 32 35 3a 35 29 3b 0d 0a 20 20 20 20 57 72 69 74 65 6c 6e ┆, Integral : 25:5); Writeln┆
0x0540…0560 28 27 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 f5 20 20 27 29 3b 0d 0a 20 20 20 ┆(' '); ┆
0x0560…0580 20 57 72 69 74 65 6c 6e 28 27 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 2d ec 20 20 20 20 ┆ Writeln(' - ┆
0x0580…05a0 27 29 3b 0d 0a 20 20 45 6e 64 3b 0d 0a 65 6e 64 3b 0d 0a 0d 0a 62 65 67 69 6e 0d 0a 20 20 52 65 ┆'); End; end; begin Re┆
0x05a0…05c0 70 65 61 74 0d 0a 20 20 20 20 43 6c 72 53 63 72 3b 0d 0a 20 20 20 20 57 72 69 74 65 6c 6e 28 27 ┆peat ClrScr; Writeln('┆
0x05c0…05e0 c9 cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd ┆ ┆
0x05e0…0600 cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd bb 27 3a 36 34 29 3b 0d 0a 20 20 20 20 ┆ ':64); ┆
0x0600…0620 57 72 69 74 65 6c 6e 28 27 ba 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 4e ┆Writeln(' N┆
0x0620…0640 4f 52 4d 41 4c 46 53 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ba 27 3a 36 ┆ORMALFS ':6┆
0x0640…0660 34 29 3b 0d 0a 20 20 20 20 57 72 69 74 65 6c 6e 28 27 c7 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 ┆4); Writeln(' ┆
0x0660…0680 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 c4 ┆ ┆
0x0680…06a0 c4 c4 c4 c4 c4 b6 27 3a 36 34 29 3b 0d 0a 20 20 20 20 57 72 69 74 65 6c 6e 28 27 ba 20 55 64 72 ┆ ':64); Writeln(' Udr┆
0x06a0…06c0 65 67 6e 65 72 20 69 6e 74 65 67 72 61 6c 20 66 6f 72 20 6e 6f 72 6d 61 6c 66 6f 72 64 65 6c 69 ┆egner integral for normalfordeli┆
0x06c0…06e0 6e 67 73 66 75 6e 6b 74 69 6f 6e 65 6e 20 ba 27 3a 36 34 29 3b 0d 0a 20 20 20 20 57 72 69 74 65 ┆ngsfunktionen ':64); Write┆
0x06e0…0700 6c 6e 28 27 ba 20 53 74 65 66 66 65 6e 20 42 65 72 69 6e 67 20 4a 65 6e 73 65 6e 20 2f 20 50 72 ┆ln(' Steffen Bering Jensen / Pr┆
0x0700…0720 6f 6a 65 6b 74 20 55 20 20 20 20 20 41 70 72 69 6c 20 31 39 39 30 20 ba 27 3a 36 34 29 3b 0d 0a ┆ojekt U April 1990 ':64); ┆
0x0720…0740 20 20 20 20 57 72 69 74 65 6c 6e 28 27 c8 cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd ┆ Writeln(' ┆
0x0740…0760 cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd cd ┆ ┆
0x0760…0780 bc 27 3a 36 34 29 3b 0d 0a 20 20 20 20 55 70 70 65 72 4c 69 6d 69 74 20 3a 3d 20 30 3b 0d 0a 20 ┆ ':64); UpperLimit := 0; ┆
0x0780…07a0 20 20 20 49 6e 74 65 67 72 61 6c 20 3a 3d 20 30 3b 0d 0a 20 20 20 20 45 72 72 6f 72 20 3a 3d 20 ┆ Integral := 0; Error := ┆
0x07a0…07c0 30 3b 0d 0a 20 20 20 20 47 65 74 4c 69 6d 69 74 28 55 70 70 65 72 4c 69 6d 69 74 29 3b 0d 0a 20 ┆0; GetLimit(UpperLimit); ┆
0x07c0…07e0 20 20 20 53 69 6d 70 73 6f 6e 28 4c 6f 77 65 72 4c 69 6d 69 74 2c 20 55 70 70 65 72 4c 69 6d 69 ┆ Simpson(LowerLimit, UpperLimi┆
0x07e0…0800 74 2c 20 4e 75 6d 49 6e 74 65 72 76 61 6c 73 2c 20 49 6e 74 65 67 72 61 6c 2c 20 45 72 72 6f 72 ┆t, NumIntervals, Integral, Error┆
0x0800…0820 2c 20 40 70 68 69 29 3b 0d 0a 20 20 20 20 52 65 73 75 6c 74 73 28 55 70 70 65 72 4c 69 6d 69 74 ┆, @phi); Results(UpperLimit┆
0x0820…0840 2c 20 49 6e 74 65 67 72 61 6c 2c 20 45 72 72 6f 72 29 3b 0d 0a 20 20 20 20 57 72 69 74 65 6c 6e ┆, Integral, Error); Writeln┆
0x0840…0860 3b 0d 0a 20 20 20 20 57 72 69 74 65 6c 6e 3b 0d 0a 20 20 20 20 57 72 69 74 65 28 27 20 20 20 20 ┆; Writeln; Write(' ┆
0x0860…0880 20 20 20 20 20 20 20 20 53 6b 61 6c 20 64 65 72 20 6c 61 76 65 73 20 66 6c 65 72 65 20 75 64 72 ┆ Skal der laves flere udr┆
0x0880…08a0 65 67 6e 69 6e 67 65 72 20 4a 2f 4e 20 3f 20 27 29 3b 0d 0a 20 20 20 20 52 65 70 65 61 74 0d 0a ┆egninger J/N ? '); Repeat ┆
0x08a0…08c0 20 20 20 20 20 20 69 63 68 3a 3d 52 65 61 64 4b 65 79 3b 0d 0a 20 20 20 20 55 6e 74 69 6c 20 69 ┆ ich:=ReadKey; Until i┆
0x08c0…08e0 63 68 20 49 4e 20 5b 27 6a 27 2c 27 4a 27 2c 27 6e 27 2c 27 4e 27 5d 3b 0d 0a 20 20 55 6e 74 69 ┆ch IN Æ'j','J','n','N'Å; Unti┆
0x08e0…0900 6c 20 69 63 68 20 49 4e 20 5b 27 6e 27 2c 27 4e 27 5d 3b 0d 0a 45 6e 64 2e 0d 0a 0d 0a 70 72 6f ┆l ich IN Æ'n','N'Å; End. pro┆
0x0900…0920 63 65 64 75 72 65 20 53 69 6d 70 73 6f 6e 28 4c 6f 77 65 72 4c 69 6d 69 74 20 20 20 3a 20 46 6c ┆cedure Simpson(LowerLimit : Fl┆
0x0920…0940 6f 61 74 3b 0d 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 55 70 70 65 72 4c 69 6d ┆oat; UpperLim┆
0x0940…0960 69 74 20 20 20 3a 20 46 6c 6f 61 74 3b 0d 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆it : Float; ┆
0x0960…0980 20 4e 75 6d 49 6e 74 65 72 76 61 6c 73 20 3a 20 69 6e 74 65 67 65 72 3b 0d 0a 20 20 20 20 20 20 ┆ NumIntervals : integer; ┆
0x0980…09a0 20 20 20 20 20 20 20 20 76 61 72 20 49 6e 74 65 67 72 61 6c 20 20 20 20 20 3a 20 46 6c 6f 61 74 ┆ var Integral : Float┆
0x09a0…09c0 3b 0d 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 76 61 72 20 45 72 72 6f 72 20 20 20 20 20 20 ┆; var Error ┆
0x09c0…09e0 20 20 3a 20 62 79 74 65 3b 0d 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 46 75 6e ┆ : byte; Fun┆
0x09e0…0a00 63 50 74 72 20 20 20 20 20 20 3a 20 50 6f 69 6e 74 65 72 29 3b 0d 0a 0d 0a 7b 2d 2d 2d 2d 2d 2d ┆cPtr : Pointer); æ------┆
0x0a00…0a20 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d ┆--------------------------------┆
[…0x1…]
0x0a40…0a60 2d 2d 2d 2d 2d 2d 7d 0d 0a 7b 2d 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆------å æ- ┆
0x0a60…0a80 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ ┆
0x0a80…0aa0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 50 75 72 ┆ -å æ- Pur┆
0x0aa0…0ac0 70 6f 73 65 3a 20 67 69 76 65 6e 20 61 20 66 75 6e 63 74 69 6f 6e 2c 20 54 4e 54 61 72 67 65 74 ┆pose: given a function, TNTarget┆
0x0ac0…0ae0 46 28 58 29 2c 20 63 6f 6d 70 75 74 65 20 74 68 65 20 69 6e 74 65 67 72 61 6c 20 6f 66 20 20 20 ┆F(X), compute the integral of ┆
0x0ae0…0b00 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 20 20 20 20 20 20 54 4e 54 61 72 67 65 74 46 28 ┆ -å æ- TNTargetF(┆
0x0b00…0b20 58 29 20 66 72 6f 6d 20 4c 6f 77 65 72 4c 69 6d 69 74 20 74 6f 20 55 70 70 65 72 4c 69 6d 69 74 ┆X) from LowerLimit to UpperLimit┆
0x0b20…0b40 20 75 73 69 6e 67 20 53 69 6d 70 73 6f 6e 20 20 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 ┆ using Simpson -å æ- ┆
0x0b40…0b60 20 20 20 20 20 20 43 6f 6d 70 6f 73 69 74 65 20 41 6c 67 6f 72 69 74 68 6d 2e 20 20 20 20 20 20 ┆ Composite Algorithm. ┆
0x0b60…0b80 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ ┆
0x0b80…0ba0 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ -å æ- ┆
0x0ba0…0bc0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ ┆
0x0bc0…0be0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 55 73 65 ┆ -å æ- Use┆
0x0be0…0c00 72 2d 64 65 66 69 6e 65 64 20 46 75 6e 63 74 69 6f 6e 73 3a 20 54 4e 54 61 72 67 65 74 46 28 58 ┆r-defined Functions: TNTargetF(X┆
0x0c00…0c20 20 3a 20 72 65 61 6c 29 20 3a 20 72 65 61 6c 3b 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ : real) : real; ┆
0x0c20…0c40 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ -å æ- ┆
0x0c40…0c60 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ ┆
0x0c60…0c80 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 47 ┆ -å æ- G┆
0x0c80…0ca0 6c 6f 62 61 6c 20 56 61 72 69 61 62 6c 65 73 3a 20 20 4c 6f 77 65 72 4c 69 6d 69 74 20 20 20 3a ┆lobal Variables: LowerLimit :┆
0x0ca0…0cc0 20 72 65 61 6c 3b 20 20 20 20 4c 6f 77 65 72 20 6c 69 6d 69 74 20 6f 66 20 69 6e 74 65 67 72 61 ┆ real; Lower limit of integra┆
0x0cc0…0ce0 74 69 6f 6e 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆tion -å æ- ┆
0x0ce0…0d00 20 20 55 70 70 65 72 4c 69 6d 69 74 20 20 20 3a 20 72 65 61 6c 3b 20 20 20 20 55 70 70 65 72 20 ┆ UpperLimit : real; Upper ┆
0x0d00…0d20 6c 69 6d 69 74 20 6f 66 20 69 6e 74 65 67 72 61 74 69 6f 6e 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 ┆limit of integration -å æ- ┆
0x0d20…0d40 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 4e 75 6d 49 6e 74 65 72 76 61 6c 73 20 3a ┆ NumIntervals :┆
0x0d40…0d60 20 69 6e 74 65 67 65 72 3b 20 4e 75 6d 62 65 72 20 6f 66 20 73 75 62 69 6e 74 65 72 76 61 6c 73 ┆ integer; Number of subintervals┆
0x0d60…0d80 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ -å æ- ┆
0x0d80…0da0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 6f 76 65 72 20 77 ┆ over w┆
0x0da0…0dc0 68 69 63 68 20 74 6f 20 75 73 65 20 20 20 20 20 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 ┆hich to use -å æ- ┆
0x0dc0…0de0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ ┆
0x0de0…0e00 20 20 20 20 20 20 20 20 20 20 53 69 6d 70 73 6f 6e 27 73 20 72 75 6c 65 2e 20 20 20 20 20 20 20 ┆ Simpson's rule. ┆
0x0e00…0e20 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ -å æ- ┆
0x0e20…0e40 20 20 49 6e 74 65 67 72 61 6c 20 20 20 20 20 3a 20 72 65 61 6c 3b 20 20 20 20 56 61 6c 75 65 20 ┆ Integral : real; Value ┆
0x0e40…0e60 6f 66 20 74 68 65 20 69 6e 74 65 67 72 61 6c 20 6f 66 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 ┆of the integral of -å æ- ┆
0x0e60…0e80 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ ┆
0x0e80…0ea0 20 20 20 20 20 20 20 20 20 20 54 4e 54 61 72 67 65 74 46 20 6f 76 65 72 20 74 68 65 20 67 69 76 ┆ TNTargetF over the giv┆
0x0ea0…0ec0 65 6e 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆en -å æ- ┆
0x0ec0…0ee0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 69 6e 74 65 72 76 ┆ interv┆
0x0ee0…0f00 61 6c 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 ┆al -å æ- ┆
0x0f00…0f20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 45 72 72 6f 72 20 20 20 20 20 20 20 20 3a ┆ Error :┆
0x0f20…0f40 20 62 79 74 65 3b 20 20 20 20 46 6c 61 67 73 20 69 66 20 73 6f 6d 65 74 68 69 6e 67 20 67 6f 65 ┆ byte; Flags if something goe┆
0x0f40…0f60 73 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆s -å æ- ┆
0x0f60…0f80 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 77 72 6f 6e 67 20 ┆ wrong ┆
0x0f80…0fa0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 ┆ -å æ- ┆
0x0fa0…0fc0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ ┆
[…0x1…]
0x0fe0…1000 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 20 20 20 20 20 20 20 20 20 45 72 72 6f 72 73 3a ┆ -å æ- Errors:┆
0x1000…1020 20 20 30 3a 20 4e 6f 20 65 72 72 6f 72 73 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ 0: No errors ┆
0x1020…1040 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 ┆ -å æ- ┆
0x1040…1060 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 31 3a 20 4e 75 6d 49 6e 74 65 72 76 61 6c ┆ 1: NumInterval┆
0x1060…1080 73 20 3c 20 30 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆s < 0 ┆
0x1080…10a0 20 20 20 20 20 2d 7d 0d 0a 7b 2d 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ -å æ- ┆
0x10a0…10c0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆ ┆
0x10c0…10e0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 2d 7d 0d 0a 7b 2d 2d 2d 2d 2d 2d ┆ -å æ------┆
0x10e0…1100 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d 2d ┆--------------------------------┆
[…0x1…]
0x1120…1140 2d 2d 2d 2d 2d 2d 7d 0d 0a 0d 0a 70 72 6f 63 65 64 75 72 65 20 53 69 6d 70 73 6f 6e 7b 28 4c 6f ┆------å procedure Simpsonæ(Lo┆
0x1140…1160 77 65 72 4c 69 6d 69 74 20 20 3a 20 46 6c 6f 61 74 3b 0d 0a 20 20 20 20 20 20 20 20 20 20 20 20 ┆werLimit : Float; ┆
0x1160…1180 20 20 20 20 20 20 55 70 70 65 72 4c 69 6d 69 74 20 20 20 3a 20 46 6c 6f 61 74 3b 0d 0a 20 20 20 ┆ UpperLimit : Float; ┆
0x1180…11a0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 4e 75 6d 49 6e 74 65 72 76 61 6c 73 20 3a 20 69 6e ┆ NumIntervals : in┆
0x11a0…11c0 74 65 67 65 72 3b 0d 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 76 61 72 20 49 6e 74 65 67 72 ┆teger; var Integr┆
0x11c0…11e0 61 6c 20 20 20 20 20 3a 20 46 6c 6f 61 74 3b 0d 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 76 ┆al : Float; v┆
0x11e0…1200 61 72 20 45 72 72 6f 72 20 20 20 20 20 20 20 20 3a 20 62 79 74 65 3b 0d 0a 20 20 20 20 20 20 20 ┆ar Error : byte; ┆
0x1200…1220 20 20 20 20 20 20 20 20 20 20 20 46 75 6e 63 50 74 72 20 20 20 20 20 20 3a 20 50 6f 69 6e 74 65 ┆ FuncPtr : Pointe┆
0x1220…1240 72 29 7d 3b 0d 0a 0d 0a 76 61 72 0d 0a 20 20 53 70 61 63 69 6e 67 20 3a 20 46 6c 6f 61 74 3b 20 ┆r)å; var Spacing : Float; ┆
0x1240…1260 20 20 20 20 20 20 20 20 20 20 7b 20 53 69 7a 65 20 6f 66 20 65 61 63 68 20 73 75 62 69 6e 74 65 ┆ æ Size of each subinte┆
0x1260…1280 72 76 61 6c 20 7d 0d 0a 20 20 50 6f 69 6e 74 20 3a 20 46 6c 6f 61 74 3b 20 20 20 20 20 20 20 20 ┆rval å Point : Float; ┆
0x1280…12a0 20 20 20 20 20 7b 20 4d 69 64 77 61 79 20 70 6f 69 6e 74 20 6f 66 20 65 61 63 68 20 69 6e 74 65 ┆ æ Midway point of each inte┆
0x12a0…12c0 72 76 61 6c 20 7d 0d 0a 20 20 4f 64 64 53 75 6d 2c 20 45 76 65 6e 53 75 6d 20 3a 20 46 6c 6f 61 ┆rval å OddSum, EvenSum : Floa┆
0x12c0…12e0 74 3b 20 20 20 7b 20 53 75 6d 73 20 6f 66 20 76 61 6c 75 65 73 20 6f 76 65 72 20 6f 64 64 20 6e ┆t; æ Sums of values over odd n┆
0x12e0…1300 75 6d 62 65 72 65 64 20 7d 0d 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 ┆umbered å ┆
0x1300…1320 20 20 20 20 20 20 20 7b 20 61 6e 64 20 65 76 65 6e 20 6e 75 6d 62 65 72 65 64 20 73 75 62 69 6e ┆ æ and even numbered subin┆
0x1320…1340 74 65 72 76 61 6c 73 20 7d 0d 0a 20 20 4c 69 6d 69 74 73 56 61 6c 75 65 20 3a 20 46 6c 6f 61 74 ┆tervals å LimitsValue : Float┆
0x1340…1360 3b 20 20 20 20 20 20 20 7b 20 53 75 6d 20 6f 66 20 76 61 6c 75 65 73 20 61 74 20 74 68 65 20 65 ┆; æ Sum of values at the e┆
0x1360…1380 6e 64 70 6f 69 6e 74 73 20 7d 0d 0a 20 20 49 6e 74 65 72 76 61 6c 20 3a 20 69 6e 74 65 67 65 72 ┆ndpoints å Interval : integer┆
0x1380…13a0 3b 20 20 20 20 20 20 20 7b 20 43 6f 75 6e 74 65 72 20 7d 0d 0a 0d 0a 62 65 67 69 6e 20 7b 20 70 ┆; æ Counter å begin æ p┆
0x13a0…13c0 72 6f 63 65 64 75 72 65 20 53 69 6d 70 73 6f 6e 20 7d 0d 0a 20 20 69 66 20 4e 75 6d 49 6e 74 65 ┆rocedure Simpson å if NumInte┆
0x13c0…13e0 72 76 61 6c 73 20 3c 3d 20 30 20 74 68 65 6e 0d 0a 20 20 20 20 45 72 72 6f 72 20 3a 3d 20 31 0d ┆rvals <= 0 then Error := 1 ┆
0x13e0…1400 0a 20 20 65 6c 73 65 0d 0a 20 20 20 20 62 65 67 69 6e 0d 0a 20 20 20 20 20 20 53 70 61 63 69 6e ┆ else begin Spacin┆
0x1400…1420 67 20 3a 3d 20 28 55 70 70 65 72 4c 69 6d 69 74 20 2d 20 4c 6f 77 65 72 4c 69 6d 69 74 29 20 2f ┆g := (UpperLimit - LowerLimit) /┆
0x1420…1440 20 28 32 20 2a 20 4e 75 6d 49 6e 74 65 72 76 61 6c 73 29 3b 0d 0a 20 20 20 20 20 20 50 6f 69 6e ┆ (2 * NumIntervals); Poin┆
0x1440…1460 74 20 3a 3d 20 4c 6f 77 65 72 4c 69 6d 69 74 3b 0d 0a 20 20 20 20 20 20 4f 64 64 53 75 6d 20 3a ┆t := LowerLimit; OddSum :┆
0x1460…1480 3d 20 30 3b 0d 0a 20 20 20 20 20 20 45 76 65 6e 53 75 6d 20 3a 3d 20 30 3b 0d 0a 20 20 20 20 20 ┆= 0; EvenSum := 0; ┆
0x1480…14a0 20 66 6f 72 20 49 6e 74 65 72 76 61 6c 20 3a 3d 20 31 20 74 6f 20 32 2a 4e 75 6d 49 6e 74 65 72 ┆ for Interval := 1 to 2*NumInter┆
0x14a0…14c0 76 61 6c 73 20 2d 20 31 20 64 6f 0d 0a 20 20 20 20 20 20 62 65 67 69 6e 0d 0a 20 20 20 20 20 20 ┆vals - 1 do begin ┆
0x14c0…14e0 20 20 50 6f 69 6e 74 20 3a 3d 20 50 6f 69 6e 74 20 2b 20 53 70 61 63 69 6e 67 3b 0d 0a 20 20 20 ┆ Point := Point + Spacing; ┆
0x14e0…1500 20 20 20 20 20 69 66 20 4f 64 64 28 49 6e 74 65 72 76 61 6c 29 20 74 68 65 6e 0d 0a 20 20 20 20 ┆ if Odd(Interval) then ┆
0x1500…1520 20 20 20 20 20 20 4f 64 64 53 75 6d 20 3a 3d 20 4f 64 64 53 75 6d 20 2b 20 55 73 65 72 46 75 6e ┆ OddSum := OddSum + UserFun┆
0x1520…1540 63 74 69 6f 6e 28 50 6f 69 6e 74 2c 20 46 75 6e 63 50 74 72 29 0d 0a 20 20 20 20 20 20 20 20 65 ┆ction(Point, FuncPtr) e┆
0x1540…1560 6c 73 65 0d 0a 20 20 20 20 20 20 20 20 20 20 45 76 65 6e 53 75 6d 20 3a 3d 20 45 76 65 6e 53 75 ┆lse EvenSum := EvenSu┆
0x1560…1580 6d 20 2b 20 55 73 65 72 46 75 6e 63 74 69 6f 6e 28 50 6f 69 6e 74 2c 20 46 75 6e 63 50 74 72 29 ┆m + UserFunction(Point, FuncPtr)┆
0x1580…15a0 3b 0d 0a 20 20 20 20 20 20 65 6e 64 3b 0d 0a 20 20 20 20 20 20 4c 69 6d 69 74 73 56 61 6c 75 65 ┆; end; LimitsValue┆
0x15a0…15c0 20 3a 3d 20 55 73 65 72 46 75 6e 63 74 69 6f 6e 28 55 70 70 65 72 4c 69 6d 69 74 2c 20 46 75 6e ┆ := UserFunction(UpperLimit, Fun┆
0x15c0…15e0 63 50 74 72 29 20 2b 0d 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 55 73 ┆cPtr) + Us┆
0x15e0…1600 65 72 46 75 6e 63 74 69 6f 6e 28 4c 6f 77 65 72 4c 69 6d 69 74 2c 20 46 75 6e 63 50 74 72 29 3b ┆erFunction(LowerLimit, FuncPtr);┆
0x1600…1620 0d 0a 20 20 20 20 20 20 49 6e 74 65 67 72 61 6c 20 3a 3d 20 53 70 61 63 69 6e 67 20 2a 20 28 4c ┆ Integral := Spacing * (L┆
0x1620…1640 69 6d 69 74 73 56 61 6c 75 65 20 2b 20 32 20 2a 20 45 76 65 6e 53 75 6d 20 2b 20 34 20 2a 20 4f ┆imitsValue + 2 * EvenSum + 4 * O┆
0x1640…1660 64 64 53 75 6d 29 20 2f 20 33 3b 0d 0a 20 20 20 20 65 6e 64 3b 0d 0a 65 6e 64 3b 20 7b 20 70 72 ┆ddSum) / 3; end; end; æ pr┆
0x1660…1674 6f 63 65 64 75 72 65 20 53 69 6d 70 73 6f 6e 20 7d 0d 0a 1a ┆ocedure Simpson å ┆