home *** CD-ROM | disk | FTP | other *** search
/ ARM Club 1 / ARM_CLUB_CD.iso / contents / apps / maths / progs / programs / Maths2 / Brian / Simpson < prev    next >
Unknown  |  1990-04-18  |  2.0 KB
open in: MacOS 8.1     |     Win98     |     DOS

view JSON data     |     view as text


This file was not able to be converted.
This format is not currently supported by dexvert.

ConfidenceProgramDetectionMatch TypeSupport
10% dexvert Truevision Targa Graphic (image/tga) magic Supported
100% file data default
100% xdgMime image/x-tga default (weak)


hex view
+--------+-------------------------+-------------------------+--------+--------+
|00000000| 0d 00 0a 38 f1 22 45 6e | 74 65 72 20 66 75 6e 63 |...8."En|ter func|
|00000010| 74 69 6f 6e 20 74 6f 20 | 62 65 20 69 6e 74 65 67 |tion to |be integ|
|00000020| 72 61 74 65 64 2c 20 77 | 69 74 68 20 61 72 67 75 |rated, w|ith argu|
|00000030| 6d 65 6e 74 20 58 22 27 | 0d 00 14 12 e8 20 22 46 |ment X"'|..... "F|
|00000040| 28 58 29 20 3d 20 22 20 | 46 24 0d 00 1e 08 f1 20 |(X) = " |F$..... |
|00000050| 22 22 0d 00 28 16 e8 20 | 22 4c 6f 77 65 72 20 6c |""..(.. |"Lower l|
|00000060| 69 6d 69 74 20 20 22 41 | 0d 00 32 16 e8 20 22 55 |imit  "A|..2.. "U|
|00000070| 70 70 65 72 20 6c 69 6d | 69 74 20 20 22 42 0d 00 |pper lim|it  "B..|
|00000080| 3c 34 f1 20 22 20 20 20 | 20 20 20 20 20 20 20 20 |<4. "   |        |
|00000090| 20 20 20 20 20 20 20 20 | 20 49 6e 74 65 67 72 61 |        | Integra|
|000000a0| 6c 20 3d 20 20 22 3b a4 | 71 73 69 6d 70 28 41 2c |l =  ";.|qsimp(A,|
|000000b0| 42 29 0d 00 46 05 e0 0d | 00 50 49 f4 20 3a 2d 2d |B)..F...|.PI. :--|
|000000c0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000000d0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000000e0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000000f0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00000100| 0d 00 5a 0d dd 20 a4 61 | 72 67 28 58 29 0d 00 64 |..Z.. .a|rg(X)..d|
|00000110| 0a 3d a0 28 46 24 29 0d | 00 6e 49 f4 20 3a 2d 2d |.=.(F$).|.nI. :--|
|00000120| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00000130| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00000140| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00000150| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00000160| 0d 00 78 11 dd 20 a4 71 | 73 69 6d 70 28 41 2c 42 |..x.. .q|simp(A,B|
|00000170| 29 0d 00 82 40 f4 20 52 | 65 74 75 72 6e 73 20 74 |)...@. R|eturns t|
|00000180| 68 65 20 69 6e 74 65 67 | 72 61 6c 20 6f 66 20 74 |he integ|ral of t|
|00000190| 68 65 20 66 75 6e 63 74 | 69 6f 6e 20 46 4e 61 72 |he funct|ion FNar|
|000001a0| 67 28 58 29 20 66 72 6f | 6d 20 41 20 74 6f 20 42 |g(X) fro|m A to B|
|000001b0| 2e 0d 00 8c 45 f4 20 54 | 68 65 20 70 61 72 61 6d |....E. T|he param|
|000001c0| 65 74 65 72 20 45 50 53 | 20 63 61 6e 20 62 65 20 |eter EPS| can be |
|000001d0| 73 65 74 20 74 6f 20 74 | 68 65 20 64 65 73 69 72 |set to t|he desir|
|000001e0| 65 64 20 66 72 61 63 74 | 69 6f 6e 61 6c 20 61 63 |ed fract|ional ac|
|000001f0| 63 75 72 61 63 79 0d 00 | 96 44 f4 20 61 6e 64 20 |curacy..|.D. and |
|00000200| 4a 4d 41 58 25 20 73 6f | 20 74 68 61 74 20 32 5e |JMAX% so| that 2^|
|00000210| 28 4a 4d 41 58 25 2d 31 | 29 20 69 73 20 74 68 65 |(JMAX%-1|) is the|
|00000220| 20 6d 61 78 69 6d 75 6d | 20 61 6c 6c 6f 77 65 64 | maximum| allowed|
|00000230| 20 6e 75 6d 62 65 72 20 | 6f 66 0d 00 a0 38 f4 20 | number |of...8. |
|00000240| 73 74 65 70 73 2e 20 20 | 49 6e 74 65 67 72 61 74 |steps.  |Integrat|
|00000250| 69 6f 6e 20 69 73 20 70 | 65 72 66 6f 72 6d 65 64 |ion is p|erformed|
|00000260| 20 62 79 20 53 69 6d 70 | 73 6f 6e 27 73 20 72 75 | by Simp|son's ru|
|00000270| 6c 65 0d 00 aa 07 f4 20 | 3a 0d 00 b4 28 ea 20 45 |le..... |:...(. E|
|00000280| 50 53 2c 20 4a 4d 41 58 | 25 2c 20 4f 53 54 2c 20 |PS, JMAX|%, OST, |
|00000290| 4f 53 2c 20 4a 25 2c 20 | 44 4f 4e 45 25 2c 20 53 |OS, J%, |DONE%, S|
|000002a0| 54 0d 00 be 16 45 50 53 | 3d 31 45 2d 36 3a 20 4a |T....EPS|=1E-6: J|
|000002b0| 4d 41 58 25 3d 32 30 0d | 00 c3 1d 20 40 25 20 3d |MAX%=20.|... @% =|
|000002c0| 20 a8 28 2e 35 2d ab 28 | 45 50 53 29 29 2a 32 35 | .(.5-.(|EPS))*25|
|000002d0| 36 2b 31 30 0d 00 c8 17 | 4f 53 54 3d 2d 31 45 33 |6+10....|OST=-1E3|
|000002e0| 30 3a 20 4f 53 3d 2d 31 | 45 33 30 0d 00 d2 07 f4 |0: OS=-1|E30.....|
|000002f0| 20 3a 0d 00 dc 14 e3 20 | 4a 25 20 3d 20 31 20 b8 | :..... |J% = 1 .|
|00000300| 20 4a 4d 41 58 25 0d 00 | e6 1a 20 20 53 54 20 3d | JMAX%..|..  ST =|
|00000310| 20 a4 74 72 61 70 7a 64 | 28 41 2c 42 2c 4a 25 29 | .trapzd|(A,B,J%)|
|00000320| 0d 00 f0 15 20 20 53 53 | 3d 28 34 2a 53 54 2d 4f |....  SS|=(4*ST-O|
|00000330| 53 54 29 2f 33 0d 00 fa | 32 20 20 e7 20 94 28 53 |ST)/3...|2  . .(S|
|00000340| 53 2d 4f 53 29 20 3c 20 | 45 50 53 2a 94 28 4f 53 |S-OS) < |EPS*.(OS|
|00000350| 29 20 8c 20 44 4f 4e 45 | 25 20 3d 20 b9 3a 20 4a |) . DONE|% = .: J|
|00000360| 25 3d 4a 4d 41 58 25 0d | 01 04 0b 20 20 4f 53 3d |%=JMAX%.|...  OS=|
|00000370| 53 53 0d 01 0e 0c 20 20 | 4f 53 54 3d 53 54 0d 01 |SS....  |OST=ST..|
|00000380| 18 08 ed 20 4a 25 0d 01 | 22 15 e7 20 44 4f 4e 45 |... J%..|".. DONE|
|00000390| 25 20 3d 20 b9 20 8c 20 | 3d 53 53 0d 01 2c 19 f1 |% = . . |=SS..,..|
|000003a0| 20 22 54 6f 6f 20 6d 61 | 6e 79 20 73 74 65 70 73 | "Too ma|ny steps|
|000003b0| 22 3a 20 fa 0d 01 36 4a | f4 20 3a 2d 2d 2d 2d 2d |": ...6J|. :-----|
|000003c0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000003d0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000003e0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000003f0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 0d 01 |--------|------..|
|00000400| 40 15 dd 20 a4 74 72 61 | 70 7a 64 28 41 2c 42 2c |@.. .tra|pzd(A,B,|
|00000410| 4e 25 29 0d 01 4a 07 f4 | 20 3a 0d 01 54 3e f4 20 |N%)..J..| :..T>. |
|00000420| 54 68 69 73 20 72 6f 75 | 74 69 6e 65 20 63 6f 6d |This rou|tine com|
|00000430| 70 75 74 65 73 20 74 68 | 65 20 4e 25 74 68 20 73 |putes th|e N%th s|
|00000440| 74 61 67 65 20 6f 66 20 | 72 65 66 69 6e 65 6d 65 |tage of |refineme|
|00000450| 6e 74 20 6f 66 20 61 6e | 0d 01 5e 40 f4 20 65 78 |nt of an|..^@. ex|
|00000460| 74 65 6e 64 65 64 20 74 | 72 61 70 65 7a 6f 69 64 |tended t|rapezoid|
|00000470| 61 6c 20 72 75 6c 65 2e | 20 20 50 72 6f 67 72 61 |al rule.|  Progra|
|00000480| 6d 73 20 63 61 6c 6c 69 | 6e 67 20 46 4e 74 72 61 |ms calli|ng FNtra|
|00000490| 70 7a 64 20 6d 75 73 74 | 0d 01 68 43 f4 20 70 72 |pzd must|..hC. pr|
|000004a0| 6f 76 69 64 65 20 61 20 | 66 75 6e 63 74 69 6f 6e |ovide a |function|
|000004b0| 20 46 4e 61 72 67 28 58 | 29 20 77 68 69 63 68 20 | FNarg(X|) which |
|000004c0| 69 73 20 74 6f 20 62 65 | 20 69 6e 74 65 67 72 61 |is to be| integra|
|000004d0| 74 65 64 20 62 65 74 77 | 65 65 6e 0d 01 72 41 f4 |ted betw|een..rA.|
|000004e0| 20 6c 69 6d 69 74 73 20 | 41 20 61 6e 64 20 42 2e | limits |A and B.|
|000004f0| 20 20 57 68 65 6e 20 63 | 61 6c 6c 65 64 20 77 69 |  When c|alled wi|
|00000500| 74 68 20 4e 25 3d 31 2c | 20 74 68 65 20 72 6f 75 |th N%=1,| the rou|
|00000510| 74 69 6e 65 20 72 65 74 | 75 72 6e 73 0d 01 7c 42 |tine ret|urns..|B|
|00000520| f4 20 74 68 65 20 63 72 | 75 64 65 73 74 20 65 73 |. the cr|udest es|
|00000530| 74 69 6d 61 74 65 20 6f | 66 20 74 68 65 20 69 6e |timate o|f the in|
|00000540| 74 65 67 72 61 6c 2e 20 | 20 53 75 62 73 65 71 75 |tegral. | Subsequ|
|00000550| 65 6e 74 20 63 61 6c 6c | 73 20 77 69 74 68 0d 01 |ent call|s with..|
|00000560| 86 45 f4 20 4e 25 3d 32 | 2c 33 2c 2e 2e 2e 20 28 |.E. N%=2|,3,... (|
|00000570| 69 6e 20 74 68 61 74 20 | 73 65 71 75 65 6e 74 69 |in that |sequenti|
|00000580| 61 6c 20 6f 72 64 65 72 | 29 20 77 69 6c 6c 20 69 |al order|) will i|
|00000590| 6d 70 72 6f 76 65 20 74 | 68 65 20 61 63 63 75 72 |mprove t|he accur|
|000005a0| 61 63 79 0d 01 90 43 f4 | 20 6f 66 20 53 20 62 79 |acy...C.| of S by|
|000005b0| 20 61 64 64 69 6e 67 20 | 32 5e 28 4e 25 2d 32 29 | adding |2^(N%-2)|
|000005c0| 20 61 64 64 69 74 69 6f | 6e 61 6c 20 69 6e 74 65 | additio|nal inte|
|000005d0| 72 69 6f 72 20 70 6f 69 | 6e 74 73 2e 20 20 53 20 |rior poi|nts.  S |
|000005e0| 73 68 6f 75 6c 64 0d 01 | 9a 40 f4 20 6e 6f 74 20 |should..|.@. not |
|000005f0| 62 65 20 6d 6f 64 69 66 | 69 65 64 20 62 65 74 77 |be modif|ied betw|
|00000600| 65 65 6e 20 73 75 62 73 | 65 71 75 65 6e 74 20 63 |een subs|equent c|
|00000610| 61 6c 6c 73 2e 20 20 4e | 6f 74 65 20 53 20 6d 75 |alls.  N|ote S mu|
|00000620| 73 74 20 6e 6f 74 0d 01 | a4 47 f4 20 62 65 20 6c |st not..|.G. be l|
|00000630| 6f 63 61 6c 20 2d 20 73 | 6f 20 62 65 20 63 61 72 |ocal - s|o be car|
|00000640| 65 66 75 6c 2e 20 2d 20 | 4d 6f 64 69 66 69 65 64 |eful. - |Modified|
|00000650| 2c 20 66 72 6f 6d 20 4e | 75 6d 65 72 69 63 61 6c |, from N|umerical|
|00000660| 20 52 65 63 69 70 65 73 | 20 70 31 31 31 0d 01 ae | Recipes| p111...|
|00000670| 07 f4 20 3a 0d 01 b8 13 | ea 20 58 2c 20 53 4d 2c |.. :....|. X, SM,|
|00000680| 20 44 4c 2c 20 4a 25 0d | 01 c2 07 f4 20 3a 0d 01 | DL, J%.|.... :..|
|00000690| cc 0c e7 20 4e 25 3d 31 | 20 8c 0d 01 d6 23 20 20 |... N%=1| ....#  |
|000006a0| 53 3d 30 2e 35 2a 28 42 | 2d 41 29 2a 28 a4 61 72 |S=0.5*(B|-A)*(.ar|
|000006b0| 67 28 41 29 2b a4 61 72 | 67 28 42 29 29 0d 01 e0 |g(A)+.ar|g(B))...|
|000006c0| 07 20 20 cc 0d 01 ea 46 | 20 20 49 54 25 20 3d 20 |.  ....F|  IT% = |
|000006d0| 32 5e 28 4e 25 2d 32 29 | 20 3a f4 20 49 54 25 20 |2^(N%-2)| :. IT% |
|000006e0| 69 73 20 6e 75 6d 62 65 | 72 20 6f 66 20 70 6f 69 |is numbe|r of poi|
|000006f0| 6e 74 73 20 74 6f 20 62 | 65 20 61 64 64 65 64 20 |nts to b|e added |
|00000700| 4e 45 58 54 20 63 61 6c | 6c 2e 0d 01 f4 44 20 20 |NEXT cal|l....D  |
|00000710| 44 4c 3d 28 42 2d 41 29 | 2f 49 54 25 20 3a f4 20 |DL=(B-A)|/IT% :. |
|00000720| 54 68 69 73 20 69 73 20 | 74 68 65 20 73 70 61 63 |This is |the spac|
|00000730| 69 6e 67 20 6f 66 20 74 | 68 65 20 70 6f 69 6e 74 |ing of t|he point|
|00000740| 73 20 74 6f 20 62 65 20 | 61 64 64 65 64 2e 0d 01 |s to be |added...|
|00000750| fe 10 20 20 58 3d 41 2b | 30 2e 35 2a 44 4c 0d 02 |..  X=A+|0.5*DL..|
|00000760| 08 0a 20 20 53 4d 3d 30 | 0d 02 12 14 20 20 e3 20 |..  SM=0|....  . |
|00000770| 4a 25 20 3d 20 31 20 b8 | 20 49 54 25 0d 02 1c 19 |J% = 1 .| IT%....|
|00000780| 20 20 20 20 53 4d 20 3d | 20 53 4d 20 2b 20 a4 61 |    SM =| SM + .a|
|00000790| 72 67 28 58 29 0d 02 26 | 12 20 20 20 20 58 20 3d |rg(X)..&|.    X =|
|000007a0| 20 58 20 2b 20 44 4c 0d | 02 30 0a 20 20 ed 20 4a | X + DL.|.0.  . J|
|000007b0| 25 0d 02 3a 44 20 20 53 | 3d 30 2e 35 2a 28 53 2b |%..:D  S|=0.5*(S+|
|000007c0| 28 42 2d 41 29 2a 53 4d | 2f 49 54 25 29 3a f4 20 |(B-A)*SM|/IT%):. |
|000007d0| 54 68 69 73 20 72 65 70 | 6c 61 63 65 73 20 53 20 |This rep|laces S |
|000007e0| 62 79 20 69 74 73 20 72 | 65 66 69 6e 65 64 20 76 |by its r|efined v|
|000007f0| 61 6c 75 65 2e 0d 02 44 | 05 cd 0d 02 4e 06 3d 53 |alue...D|....N.=S|
|00000800| 0d ff                   |                         |..      |        |
+--------+-------------------------+-------------------------+--------+--------+