/ The Datafile PD-CD 2 / DATAFILE_PDCD2.iso / fractals / fractal107 / userfunc / cuberoot

< prev

next >

MacOS 8.1

Win98

DOS

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 0d 00 0a 24 f4 20 3e 20 | 3c 46 72 61 63 24 44 69 |...$. > |<Frac$Di|
|00000010| 72 3e 2e 55 73 65 72 46 | 75 6e 63 2e 43 75 62 65 |r>.UserF|unc.Cube|
|00000020| 52 6f 6f 74 0d 00 14 17 | f4 20 53 2e 4d 61 72 6b |Root....|. S.Mark|
|00000030| 68 61 6d 20 31 34 2f 33 | 2f 39 32 0d 00 1e 17 f4 |ham 14/3|/92.....|
|00000040| 20 43 75 62 65 20 52 6f | 6f 74 20 50 6c 6f 74 74 | Cube Ro|ot Plott|
|00000050| 65 72 0d 00 28 29 c8 9b | 20 22 3c 46 72 61 63 24 |er..()..| "<Frac$|
|00000060| 44 69 72 3e 2e 41 73 6d | 2e 46 72 61 63 6c 69 62 |Dir>.Asm|.Fraclib|
|00000070| 22 3a f2 66 72 61 63 76 | 61 72 73 0d 00 32 25 c8 |":.fracv|ars..2%.|
|00000080| 9b 20 22 3c 46 72 61 63 | 24 44 69 72 3e 2e 41 73 |. "<Frac|$Dir>.As|
|00000090| 6d 2e 46 50 6c 69 62 22 | 3a f2 66 70 69 6e 69 74 |m.FPlib"|:.fpinit|
|000000a0| 0d 00 3c 1c ee 20 85 20 | f1 20 f6 24 3b 22 20 61 |..<.. . |. .$;" a|
|000000b0| 74 20 6c 69 6e 65 20 22 | 3b 9e 3a fa 0d 00 46 27 |t line "|;.:...F'|
|000000c0| f4 20 3d 3d 3d 46 75 6e | 63 74 69 6f 6e 20 73 70 |. ===Fun|ction sp|
|000000d0| 65 63 69 66 69 63 20 63 | 6f 64 65 20 66 6f 6c 6c |ecific c|ode foll|
|000000e0| 6f 77 73 0d 00 50 05 3a | 0d 00 5a 41 63 75 62 65 |ows..P.:|..ZAcube|
|000000f0| 72 6f 6f 74 3d 30 3a 63 | 76 61 72 73 74 61 72 74 |root=0:c|varstart|
|00000100| 3d 30 3a 63 76 61 72 65 | 6e 64 3d 30 3a 20 f4 20 |=0:cvare|nd=0: . |
|00000110| 69 6e 69 74 69 61 6c 69 | 73 65 20 74 6f 20 64 75 |initiali|se to du|
|00000120| 6d 6d 79 20 76 61 6c 75 | 65 0d 00 64 05 3a 0d 00 |mmy valu|e..d.:..|
|00000130| 6e 26 f4 20 3d 3d 3d 45 | 6e 64 20 6f 66 20 66 75 |n&. ===E|nd of fu|
|00000140| 6e 63 74 69 6f 6e 20 73 | 70 65 63 69 66 69 63 20 |nction s|pecific |
|00000150| 63 6f 64 65 0d 00 78 11 | de 20 63 6f 64 65 25 20 |code..x.|. code% |
|00000160| 31 30 30 30 30 0d 00 82 | 13 e3 20 61 73 6d 3d 34 |10000...|.. asm=4|
|00000170| 20 b8 20 37 20 88 20 33 | 0d 00 8c 11 50 25 3d 30 | . 7 . 3|....P%=0|
|00000180| 3a 4f 25 3d 63 6f 64 65 | 25 0d 00 96 0c 5b 4f 50 |:O%=code|%....[OP|
|00000190| 54 20 61 73 6d 0d 00 a0 | 3d 3b 20 46 6f 6c 6c 6f |T asm...|=; Follo|
|000001a0| 77 69 6e 67 20 69 73 20 | 61 20 6c 69 73 74 20 6f |wing is |a list o|
|000001b0| 66 20 66 72 61 63 74 61 | 6c 20 66 75 6e 63 74 69 |f fracta|l functi|
|000001c0| 6f 6e 73 20 69 6e 20 74 | 68 69 73 20 6d 6f 64 75 |ons in t|his modu|
|000001d0| 6c 65 0d 00 aa 35 a4 66 | 72 61 63 68 64 72 20 3b |le...5.f|rachdr ;|
|000001e0| 20 49 64 65 6e 74 69 66 | 79 20 6f 75 72 73 65 6c | Identif|y oursel|
|000001f0| 76 65 73 20 61 73 20 61 | 20 66 72 61 63 74 61 6c |ves as a| fractal|
|00000200| 20 6d 6f 64 75 6c 65 0d | 00 b4 05 3b 0d 00 be 29 | module.|...;...)|
|00000210| 3b 20 46 75 6e 63 74 69 | 6f 6e 20 64 65 66 69 6e |; Functi|on defin|
|00000220| 69 74 69 6f 6e 20 70 6f | 69 6e 74 65 72 20 66 6f |ition po|inter fo|
|00000230| 6c 6c 6f 77 73 0d 00 c8 | 11 44 43 44 20 72 6f 6f |llows...|.DCD roo|
|00000240| 74 5f 66 75 6e 63 0d 00 | d2 05 3b 0d 00 dc 05 3b |t_func..|..;....;|
|00000250| 0d 00 e6 1d 44 43 44 20 | 2d 31 20 3b 20 45 6e 64 |....DCD |-1 ; End|
|00000260| 20 6f 66 20 66 75 6e 63 | 74 69 6f 6e 73 0d 00 f0 | of func|tions...|
|00000270| 45 3b 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |E;------|--------|
|00000280| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00000290| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000002a0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000002b0| 2d 2d 0d 00 fa 2e 3b 20 | 4e 6f 77 20 77 65 20 64 |--....; |Now we d|
|000002c0| 65 66 69 6e 65 20 43 75 | 62 65 52 6f 6f 74 27 73 |efine Cu|beRoot's|
|000002d0| 20 63 68 61 72 61 63 74 | 65 72 69 73 74 69 63 73 | charact|eristics|
|000002e0| 0d 01 04 05 3b 0d 01 0e | 52 2e 72 6f 6f 74 5f 66 |....;...|R.root_f|
|000002f0| 75 6e 63 20 a4 66 72 61 | 63 66 75 6e 63 20 28 22 |unc .fra|cfunc ("|
|00000300| 43 75 62 65 5f 52 6f 6f | 74 73 22 2c 63 75 62 65 |Cube_Roo|ts",cube|
|00000310| 72 6f 6f 74 2c 66 5a 6f | 6f 6d 2b 66 78 79 2c 30 |root,fZo|om+fxy,0|
|00000320| 2c 35 2c 63 76 61 72 65 | 6e 64 2d 63 76 61 72 73 |,5,cvare|nd-cvars|
|00000330| 74 61 72 74 2c 30 29 0d | 01 18 05 3b 0d 01 22 24 |tart,0).|...;.."$|
|00000340| 3b 20 43 75 62 65 52 6f | 6f 74 73 20 76 61 72 69 |; CubeRo|ots vari|
|00000350| 61 62 6c 65 73 20 61 66 | 74 65 72 20 68 65 72 65 |ables af|ter here|
|00000360| 0d 01 2c 0e 2e 63 76 61 | 72 73 74 61 72 74 0d 01 |..,..cva|rstart..|
|00000370| 36 2e a4 61 66 67 32 76 | 61 72 28 41 46 47 5f 64 |6..afg2v|ar(AFG_d|
|00000380| 62 6c 2c 30 2c 38 29 3a | 2e 63 78 30 20 a4 66 70 |bl,0,8):|.cx0 .fp|
|00000390| 28 22 44 46 44 20 2d 32 | 2e 34 22 29 0d 01 40 2e |("DFD -2|.4")..@.|
|000003a0| a4 61 66 67 32 76 61 72 | 28 41 46 47 5f 64 62 6c |.afg2var|(AFG_dbl|
|000003b0| 2c 30 2c 38 29 3a 2e 63 | 79 30 20 a4 66 70 28 22 |,0,8):.c|y0 .fp("|
|000003c0| 44 46 44 20 2d 32 2e 30 | 22 29 0d 01 4a 2c a4 61 |DFD -2.0|")..J,.a|
|000003d0| 66 67 32 76 61 72 28 41 | 46 47 5f 64 62 6c 2c 30 |fg2var(A|FG_dbl,0|
|000003e0| 2c 38 29 3a 2e 63 77 20 | a4 66 70 28 22 44 46 44 |,8):.cw |.fp("DFD|
|000003f0| 20 34 2e 38 22 29 0d 01 | 54 2c a4 61 66 67 32 76 | 4.8")..|T,.afg2v|
|00000400| 61 72 28 41 46 47 5f 64 | 62 6c 2c 30 2c 38 29 3a |ar(AFG_d|bl,0,8):|
|00000410| 2e 63 68 20 a4 66 70 28 | 22 44 46 44 20 34 2e 30 |.ch .fp(|"DFD 4.0|
|00000420| 22 29 0d 01 5e 32 a4 61 | 66 67 32 76 61 72 28 41 |")..^2.a|fg2var(A|
|00000430| 46 47 5f 64 62 6c 2c 30 | 2c 38 29 3a 2e 65 72 72 |FG_dbl,0|,8):.err|
|00000440| 73 69 7a 65 20 a4 66 70 | 28 22 44 46 44 20 31 45 |size .fp|("DFD 1E|
|00000450| 2d 34 22 29 0d 01 68 17 | 2e 64 65 6e 20 a4 66 70 |-4")..h.|.den .fp|
|00000460| 28 22 44 46 44 20 30 2e | 30 22 29 0d 01 72 16 2e |("DFD 0.|0")..r..|
|00000470| 6e 72 20 a4 66 70 28 22 | 44 46 44 20 30 2e 30 22 |nr .fp("|DFD 0.0"|
|00000480| 29 0d 01 7c 16 2e 6e 69 | 20 a4 66 70 28 22 44 46 |)..|..ni| .fp("DF|
|00000490| 44 20 30 2e 30 22 29 0d | 01 86 1f 2e 63 76 61 72 |D 0.0").|....cvar|
|000004a0| 65 6e 64 20 3b 20 65 6e | 64 20 6f 66 20 76 61 72 |end ; en|d of var|
|000004b0| 69 61 62 6c 65 73 0d 01 | 90 05 3b 0d 01 9a 1f 3b |iables..|..;....;|
|000004c0| 20 56 61 72 69 61 62 6c | 65 20 6e 61 6d 65 73 20 | Variabl|e names |
|000004d0| 61 66 74 65 72 20 68 65 | 72 65 0d 01 a4 1b a4 76 |after he|re.....v|
|000004e0| 61 72 6e 61 6d 65 28 22 | 20 20 20 20 20 20 20 20 |arname("|        |
|000004f0| 20 78 30 22 29 0d 01 ae | 1b a4 76 61 72 6e 61 6d | x0")...|..varnam|
|00000500| 65 28 22 20 20 20 20 20 | 20 20 20 20 79 30 22 29 |e("     |    y0")|
|00000510| 0d 01 b8 1b a4 76 61 72 | 6e 61 6d 65 28 22 20 20 |.....var|name("  |
|00000520| 20 20 20 20 77 69 64 74 | 68 22 29 0d 01 c2 1b a4 |    widt|h").....|
|00000530| 76 61 72 6e 61 6d 65 28 | 22 20 20 20 20 20 68 65 |varname(|"     he|
|00000540| 69 67 68 74 22 29 0d 01 | cc 1b a4 76 61 72 6e 61 |ight")..|...varna|
|00000550| 6d 65 28 22 20 20 20 20 | 20 20 65 72 72 6f 72 22 |me("    |  error"|
|00000560| 29 0d 01 d6 30 3b 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |)...0;--|--------|
|00000570| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00000580| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00000590| 2d 0d 01 e0 22 3b 20 54 | 68 65 20 65 6e 74 72 79 |-..."; T|he entry|
|000005a0| 20 70 6f 69 6e 74 20 74 | 6f 20 43 75 62 65 20 52 | point t|o Cube R|
|000005b0| 6f 6f 74 0d 01 ea 0d 2e | 63 75 62 65 72 6f 6f 74 |oot.....|cuberoot|
|000005c0| 0d 01 f4 18 a4 66 65 6e | 74 65 72 20 3b 20 49 6e |.....fen|ter ; In|
|000005d0| 69 74 69 61 6c 69 73 65 | 0d 01 fe 29 3b 52 30 3d |itialise|...);R0=|
|000005e0| 46 72 61 63 74 61 6c 20 | 45 76 65 6e 74 2c 20 52 |Fractal |Event, R|
|000005f0| 31 3d 61 64 64 72 65 73 | 73 20 6f 66 20 64 61 74 |1=addres|s of dat|
|00000600| 61 0d 02 08 0d 4d 4f 56 | 20 52 39 2c 23 30 0d 02 |a....MOV| R9,#0..|
|00000610| 12 2e 53 54 52 20 52 39 | 2c 5b 52 31 2c 23 66 65 |..STR R9|,[R1,#fe|
|00000620| 76 65 6e 74 5f 72 63 5d | 20 3b 20 63 6c 65 61 72 |vent_rc]| ; clear|
|00000630| 20 72 65 74 75 72 6e 20 | 63 6f 64 65 0d 02 1c 15 | return |code....|
|00000640| 43 4d 50 20 52 30 2c 23 | 66 70 6c 6f 74 69 6e 69 |CMP R0,#|fplotini|
|00000650| 74 0d 02 26 11 42 45 51 | 20 63 70 6c 6f 74 69 6e |t..&.BEQ| cplotin|
|00000660| 69 74 0d 02 30 15 43 4d | 50 20 52 30 2c 23 66 7a |it..0.CM|P R0,#fz|
|00000670| 6f 6f 6d 69 6e 69 74 0d | 02 3a 11 42 45 51 20 63 |oominit.|.:.BEQ c|
|00000680| 7a 6f 6f 6d 69 6e 69 74 | 0d 02 44 15 43 4d 50 20 |zoominit|..D.CMP |
|00000690| 52 30 2c 23 66 64 61 74 | 61 69 6e 69 74 0d 02 4e |R0,#fdat|ainit..N|
|000006a0| 11 42 45 51 20 63 64 61 | 74 61 69 6e 69 74 0d 02 |.BEQ cda|tainit..|
|000006b0| 58 14 43 4d 50 20 52 30 | 2c 23 66 64 61 74 61 63 |X.CMP R0|,#fdatac|
|000006c0| 68 6b 0d 02 62 11 42 45 | 51 20 63 64 61 74 61 63 |hk..b.BE|Q cdatac|
|000006d0| 68 63 6b 0d 02 6c 25 3b | 20 55 6e 73 75 70 70 6f |hck..l%;| Unsuppo|
|000006e0| 72 74 65 64 20 65 76 65 | 6e 74 20 2d 20 6a 75 73 |rted eve|nt - jus|
|000006f0| 74 20 72 65 74 75 72 6e | 0d 02 76 0a a4 66 65 78 |t return|..v..fex|
|00000700| 69 74 0d 02 80 05 3b 0d | 02 8a 19 3b 20 49 6e 69 |it....;.|...; Ini|
|00000710| 74 20 43 75 62 65 20 52 | 6f 6f 74 20 64 61 74 61 |t Cube R|oot data|
|00000720| 0d 02 94 0e 2e 63 64 61 | 74 61 69 6e 69 74 0d 02 |.....cda|tainit..|
|00000730| 9e 16 a4 66 70 28 22 4c | 44 46 44 20 46 30 2c 69 |...fp("L|DFD F0,i|
|00000740| 65 72 22 29 0d 02 a8 1a | a4 66 70 28 22 53 54 46 |er")....|.fp("STF|
|00000750| 44 20 46 30 2c 65 72 72 | 73 69 7a 65 22 29 0d 02 |D F0,err|size")..|
|00000760| b2 0e 2e 63 7a 6f 6f 6d | 69 6e 69 74 0d 02 bc 16 |...czoom|init....|
|00000770| a4 66 70 28 22 4c 44 46 | 44 20 46 30 2c 69 6d 78 |.fp("LDF|D F0,imx|
|00000780| 22 29 0d 02 c6 16 a4 66 | 70 28 22 53 54 46 44 20 |").....f|p("STFD |
|00000790| 46 30 2c 63 78 30 22 29 | 0d 02 d0 16 a4 66 70 28 |F0,cx0")|.....fp(|
|000007a0| 22 4c 44 46 44 20 46 30 | 2c 69 6d 79 22 29 0d 02 |"LDFD F0|,imy")..|
|000007b0| da 16 a4 66 70 28 22 53 | 54 46 44 20 46 30 2c 63 |...fp("S|TFD F0,c|
|000007c0| 79 30 22 29 0d 02 e4 16 | a4 66 70 28 22 4c 44 46 |y0")....|.fp("LDF|
|000007d0| 44 20 46 30 2c 69 6d 77 | 22 29 0d 02 ee 16 a4 66 |D F0,imw|").....f|
|000007e0| 70 28 22 53 54 46 44 20 | 46 30 2c 63 77 20 22 29 |p("STFD |F0,cw ")|
|000007f0| 0d 02 f8 16 a4 66 70 28 | 22 4c 44 46 44 20 46 30 |.....fp(|"LDFD F0|
|00000800| 2c 69 6d 68 22 29 0d 03 | 02 16 a4 66 70 28 22 53 |,imh")..|...fp("S|
|00000810| 54 46 44 20 46 30 2c 63 | 68 20 22 29 0d 03 0c 0a |TFD F0,c|h ")....|
|00000820| a4 66 65 78 69 74 0d 03 | 16 18 2e 69 6d 78 20 a4 |.fexit..|...imx .|
|00000830| 66 70 28 22 44 46 44 20 | 2d 32 2e 34 22 29 0d 03 |fp("DFD |-2.4")..|
|00000840| 20 18 2e 69 6d 79 20 a4 | 66 70 28 22 44 46 44 20 | ..imy .|fp("DFD |
|00000850| 2d 32 2e 30 22 29 0d 03 | 2a 17 2e 69 6d 77 20 a4 |-2.0")..|*..imw .|
|00000860| 66 70 28 22 44 46 44 20 | 34 2e 38 22 29 0d 03 34 |fp("DFD |4.8")..4|
|00000870| 17 2e 69 6d 68 20 a4 66 | 70 28 22 44 46 44 20 34 |..imh .f|p("DFD 4|
|00000880| 2e 30 22 29 0d 03 3e 18 | 2e 69 65 72 20 a4 66 70 |.0")..>.|.ier .fp|
|00000890| 28 22 44 46 44 20 31 45 | 2d 34 22 29 0d 03 48 05 |("DFD 1E|-4")..H.|
|000008a0| 3b 0d 03 52 17 3b 20 56 | 65 72 69 66 79 20 64 61 |;..R.; V|erify da|
|000008b0| 74 61 20 69 6e 70 75 74 | 0d 03 5c 0e 2e 63 64 61 |ta input|..\..cda|
|000008c0| 74 61 63 68 63 6b 0d 03 | 66 24 3b 76 65 72 69 66 |tachck..|f$;verif|
|000008d0| 79 20 65 72 72 6f 72 20 | 73 69 7a 65 20 6e 6f 74 |y error |size not|
|000008e0| 20 74 6f 6f 20 73 6d 61 | 6c 6c 0d 03 70 19 a4 66 | too sma|ll..p..f|
|000008f0| 70 28 22 4c 44 46 44 20 | 46 31 2c 6d 69 6e 65 72 |p("LDFD |F1,miner|
|00000900| 72 22 29 0d 03 7a 1a a4 | 66 70 28 22 4c 44 46 44 |r")..z..|fp("LDFD|
|00000910| 20 46 30 2c 65 72 72 73 | 69 7a 65 22 29 0d 03 84 | F0,errs|ize")...|
|00000920| 19 a4 66 70 28 22 4c 44 | 46 44 20 46 33 2c 6d 61 |..fp("LD|FD F3,ma|
|00000930| 78 65 72 72 22 29 0d 03 | 8e 15 a4 66 70 28 22 43 |xerr")..|...fp("C|
|00000940| 4d 46 44 20 46 30 2c 46 | 33 22 29 0d 03 98 24 a4 |MFD F0,F|3")...$.|
|00000950| 66 70 28 22 4d 56 46 47 | 54 44 20 46 30 2c 46 33 |fp("MVFG|TD F0,F3|
|00000960| 22 29 20 3b 20 73 65 74 | 20 74 6f 20 6d 61 78 0d |") ; set| to max.|
|00000970| 03 a2 0e 42 47 54 20 65 | 72 72 73 65 74 0d 03 ac |...BGT e|rrset...|
|00000980| 21 a4 66 70 28 22 43 4d | 46 44 20 46 30 2c 46 31 |!.fp("CM|FD F0,F1|
|00000990| 22 29 3b 20 74 6f 6f 20 | 73 6d 61 6c 6c 3f 0d 03 |"); too |small?..|
|000009a0| b6 24 a4 66 70 28 22 4d | 56 46 4c 54 44 20 46 30 |.$.fp("M|VFLTD F0|
|000009b0| 2c 46 31 22 29 20 3b 20 | 73 65 74 20 74 6f 20 6d |,F1") ; |set to m|
|000009c0| 69 6e 0d 03 c0 22 2e 65 | 72 72 73 65 74 20 a4 66 |in...".e|rrset .f|
|000009d0| 70 28 22 53 54 46 44 20 | 46 30 2c 65 72 72 73 69 |p("STFD |F0,errsi|
|000009e0| 7a 65 22 29 0d 03 ca 05 | 3b 0d 03 d4 0a a4 66 65 |ze")....|;.....fe|
|000009f0| 78 69 74 0d 03 de 1b 2e | 6d 61 78 65 72 72 20 a4 |xit.....|maxerr .|
|00000a00| 66 70 28 22 44 46 44 20 | 32 45 2d 32 22 29 0d 03 |fp("DFD |2E-2")..|
|00000a10| e8 1c 2e 6d 69 6e 65 72 | 72 20 a4 66 70 28 22 44 |...miner|r .fp("D|
|00000a20| 46 44 20 31 45 2d 33 30 | 22 29 0d 03 f2 05 3b 0d |FD 1E-30|")....;.|
|00000a30| 03 fc 3f 3b 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |..?;----|--------|
|00000a40| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00000a50| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00000a60| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 0d 04 |--------|------..|
|00000a70| 06 21 3b 20 49 6e 69 74 | 69 61 6c 69 73 65 20 72 |.!; Init|ialise r|
|00000a80| 65 61 64 79 20 66 6f 72 | 20 61 20 70 6c 6f 74 0d |eady for| a plot.|
|00000a90| 04 10 0e 2e 63 70 6c 6f | 74 69 6e 69 74 0d 04 1a |....cplo|tinit...|
|00000aa0| 1b 4c 44 52 20 52 30 2c | 5b 52 31 2c 23 70 6c 6f |.LDR R0,|[R1,#plo|
|00000ab0| 74 5f 66 6c 61 67 73 5d | 0d 04 24 16 84 52 20 52 |t_flags]|..$..R R|
|00000ac0| 30 2c 52 30 2c 23 70 6c | 6f 74 5f 64 62 6c 0d 04 |0,R0,#pl|ot_dbl..|
|00000ad0| 2e 29 53 54 52 20 52 30 | 2c 5b 52 31 2c 23 70 6c |.)STR R0|,[R1,#pl|
|00000ae0| 6f 74 5f 66 6c 61 67 73 | 5d 20 3b 74 65 6c 6c 20 |ot_flags|] ;tell |
|00000af0| 46 72 61 63 74 61 6c 0d | 04 38 05 3b 0d 04 42 47 |Fractal.|.8.;..BG|
|00000b00| 3b 54 68 65 72 65 20 69 | 73 20 6e 6f 20 71 75 69 |;There i|s no qui|
|00000b10| 63 6b 20 72 6f 75 74 69 | 6e 65 2e 20 54 68 65 20 |ck routi|ne. The |
|00000b20| 6e 75 6d 62 65 72 73 20 | 67 65 6e 65 72 61 74 65 |numbers |generate|
|00000b30| 64 20 61 73 20 74 68 65 | 20 69 74 65 72 61 74 69 |d as the| iterati|
|00000b40| 6f 6e 73 0d 04 4c 33 3b | 20 70 72 6f 63 65 65 64 |ons..L3;| proceed|
|00000b50| 20 61 72 65 20 74 6f 6f | 20 62 69 67 20 66 6f 72 | are too| big for|
|00000b60| 20 69 6e 74 32 37 20 74 | 6f 20 61 63 63 6f 6d 6d | int27 t|o accomm|
|00000b70| 6f 64 61 74 65 2e 0d 04 | 56 05 3b 0d 04 60 30 41 |odate...|V.;..`0A|
|00000b80| 44 52 20 52 30 2c 63 75 | 62 65 72 6f 6f 74 66 70 |DR R0,cu|berootfp|
|00000b90| 20 3b 20 75 73 65 20 66 | 6c 6f 61 74 69 6e 67 20 | ; use f|loating |
|00000ba0| 70 6f 69 6e 74 20 61 72 | 69 74 68 0d 04 6a 20 53 |point ar|ith..j S|
|00000bb0| 54 52 20 52 30 2c 72 6f | 6f 74 5f 66 75 6e 63 2b |TR R0,ro|ot_func+|
|00000bc0| 66 64 65 66 6e 5f 78 79 | 72 74 6e 0d 04 74 0a a4 |fdefn_xy|rtn..t..|
|00000bd0| 66 65 78 69 74 0d 04 7e | 05 3b 0d 04 88 32 3b 2d |fexit..~|.;...2;-|
|00000be0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00000bf0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00000c00| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 0d 04 92 05 |--------|----....|
|00000c10| 5d 0d 04 9c 44 f4 20 4e | 6f 77 20 66 6f 6c 6c 6f |]...D. N|ow follo|
|00000c20| 77 73 20 74 68 65 20 68 | 61 72 64 20 62 69 74 2c |ws the h|ard bit,|
|00000c30| 20 75 73 69 6e 67 20 74 | 68 65 20 46 50 20 65 6d | using t|he FP em|
|00000c40| 75 6c 61 74 6f 72 20 74 | 6f 20 64 6f 20 74 68 65 |ulator t|o do the|
|00000c50| 20 73 75 6d 73 0d 04 a6 | 1d f4 20 61 73 73 69 67 | sums...|.. assig|
|00000c60| 6e 20 74 68 65 20 46 50 | 20 72 65 67 69 73 74 65 |n the FP| registe|
|00000c70| 72 73 0d 04 b0 2b 52 3d | 30 3a 49 3d 31 3a 4d 3d |rs...+R=|0:I=1:M=|
|00000c80| 32 3a 4e 3d 33 3a 77 6b | 31 3d 34 3a 77 6b 32 3d |2:N=3:wk|1=4:wk2=|
|00000c90| 35 3a 77 6b 33 3d 36 3a | 77 6b 34 3d 37 0d 04 ba |5:wk3=6:|wk4=7...|
|00000ca0| 0c 5b 4f 50 54 20 61 73 | 6d 0d 04 c4 0e 2e 73 61 |.[OPT as|m.....sa|
|00000cb0| 76 65 66 6c 6f 61 74 0d | 04 ce 0c 2e 66 6c 6f 61 |vefloat.|....floa|
|00000cc0| 74 66 34 0d 04 d8 17 2e | 74 65 6d 70 78 30 20 44 |tf4.....|tempx0 D|
|00000cd0| 43 44 20 30 3a 44 43 44 | 20 30 0d 04 e2 0c 2e 66 |CD 0:DCD| 0.....f|
|00000ce0| 6c 6f 61 74 66 35 0d 04 | ec 17 2e 74 65 6d 70 79 |loatf5..|...tempy|
|00000cf0| 30 20 44 43 44 20 30 3a | 44 43 44 20 30 0d 04 f6 |0 DCD 0:|DCD 0...|
|00000d00| 18 2e 66 6c 6f 61 74 66 | 36 20 44 43 44 20 30 3a |..floatf|6 DCD 0:|
|00000d10| 44 43 44 20 30 0d 05 00 | 18 2e 66 6c 6f 61 74 66 |DCD 0...|..floatf|
|00000d20| 37 20 44 43 44 20 30 3a | 44 43 44 20 30 0d 05 0a |7 DCD 0:|DCD 0...|
|00000d30| 05 3b 0d 05 14 1c 3b 20 | 66 6c 6f 61 74 69 6e 67 |.;....; |floating|
|00000d40| 20 70 6f 69 6e 74 20 72 | 6f 75 74 69 6e 65 0d 05 | point r|outine..|
|00000d50| 1e 0f 2e 63 75 62 65 72 | 6f 6f 74 66 70 0d 05 28 |...cuber|ootfp..(|
|00000d60| 36 3b 20 43 20 70 61 73 | 73 65 73 20 46 50 20 6e |6; C pas|ses FP n|
|00000d70| 75 6d 62 65 72 73 20 69 | 6e 20 52 30 2d 52 33 21 |umbers i|n R0-R3!|
|00000d80| 2e 20 47 65 74 20 74 6f | 20 66 6c 6f 61 74 20 72 |. Get to| float r|
|00000d90| 65 67 73 0d 05 32 14 41 | 44 52 20 52 34 2c 73 61 |egs..2.A|DR R4,sa|
|00000da0| 76 65 66 6c 6f 61 74 0d | 05 3c 14 53 54 4d 49 41 |vefloat.|.<.STMIA|
|00000db0| 20 52 34 2c 7b 52 30 2d | 52 33 7d 0d 05 46 18 a4 | R4,{R0-|R3}..F..|
|00000dc0| 66 70 28 22 4c 44 46 44 | 20 52 2c 74 65 6d 70 78 |fp("LDFD| R,tempx|
|00000dd0| 30 22 29 0d 05 50 18 a4 | 66 70 28 22 4c 44 46 44 |0")..P..|fp("LDFD|
|00000de0| 20 49 2c 74 65 6d 70 79 | 30 22 29 0d 05 5a 47 3b | I,tempy|0")..ZG;|
|00000df0| 20 53 61 76 65 20 46 6c | 6f 61 74 20 72 65 67 69 | Save Fl|oat regi|
|00000e00| 73 74 65 72 73 20 46 34 | 2d 46 37 20 28 54 68 69 |sters F4|-F7 (Thi|
|00000e10| 73 20 69 73 20 41 63 6f | 72 6e 27 73 20 70 72 6f |s is Aco|rn's pro|
|00000e20| 63 65 64 75 72 65 20 63 | 6f 6e 76 65 6e 74 69 6f |cedure c|onventio|
|00000e30| 6e 29 0d 05 64 1a a4 66 | 70 28 22 53 54 46 44 20 |n)..d..f|p("STFD |
|00000e40| 46 34 2c 66 6c 6f 61 74 | 66 34 22 29 0d 05 6e 1a |F4,float|f4")..n.|
|00000e50| a4 66 70 28 22 53 54 46 | 44 20 46 35 2c 66 6c 6f |.fp("STF|D F5,flo|
|00000e60| 61 74 66 35 22 29 0d 05 | 78 1a a4 66 70 28 22 53 |atf5")..|x..fp("S|
|00000e70| 54 46 44 20 46 36 2c 66 | 6c 6f 61 74 66 36 22 29 |TFD F6,f|loatf6")|
|00000e80| 0d 05 82 1a a4 66 70 28 | 22 53 54 46 44 20 46 37 |.....fp(|"STFD F7|
|00000e90| 2c 66 6c 6f 61 74 66 37 | 22 29 0d 05 8c 21 4d 4f |,floatf7|")...!MO|
|00000ea0| 56 20 52 33 2c 23 30 20 | 3b 20 69 74 65 72 61 74 |V R3,#0 |; iterat|
|00000eb0| 69 6f 6e 20 63 6f 75 6e | 74 65 72 0d 05 96 05 3b |ion coun|ter....;|
|00000ec0| 0d 05 a0 1b a4 66 70 28 | 22 4c 44 46 44 20 77 6b |.....fp(|"LDFD wk|
|00000ed0| 31 2c 65 72 72 73 69 7a | 65 22 29 0d 05 aa 15 a4 |1,errsiz|e").....|
|00000ee0| 66 70 28 22 4d 56 46 44 | 20 77 6b 32 2c 52 22 29 |fp("MVFD| wk2,R")|
|00000ef0| 0d 05 b4 17 a4 66 70 28 | 22 41 42 53 44 20 77 6b |.....fp(|"ABSD wk|
|00000f00| 32 2c 77 6b 32 22 29 0d | 05 be 17 a4 66 70 28 22 |2,wk2").|....fp("|
|00000f10| 43 4d 46 44 20 77 6b 32 | 2c 77 6b 31 22 29 0d 05 |CMFD wk2|,wk1")..|
|00000f20| c8 22 42 47 54 20 73 74 | 61 72 74 20 3b 20 72 20 |."BGT st|art ; r |
|00000f30| 69 73 20 6e 6f 74 20 74 | 6f 6f 20 73 6d 61 6c 6c |is not t|oo small|
|00000f40| 0d 05 d2 15 a4 66 70 28 | 22 4d 56 46 44 20 77 6b |.....fp(|"MVFD wk|
|00000f50| 32 2c 49 22 29 0d 05 dc | 17 a4 66 70 28 22 41 42 |2,I")...|..fp("AB|
|00000f60| 53 44 20 77 6b 32 2c 77 | 6b 32 22 29 0d 05 e6 17 |SD wk2,w|k2")....|
|00000f70| a4 66 70 28 22 43 4d 46 | 44 20 77 6b 32 2c 77 6b |.fp("CMF|D wk2,wk|
|00000f80| 31 22 29 0d 05 f0 29 42 | 47 54 20 73 74 61 72 74 |1")...)B|GT start|
|00000f90| 20 3b 20 69 20 61 6e 64 | 20 72 20 61 72 65 20 6e | ; i and| r are n|
|00000fa0| 6f 74 20 74 6f 6f 20 73 | 6d 61 6c 6c 0d 05 fa 23 |ot too s|mall...#|
|00000fb0| 3b 69 20 61 6e 64 20 72 | 20 77 65 72 65 20 62 6f |;i and r| were bo|
|00000fc0| 74 68 20 74 6f 6f 20 73 | 6d 61 6c 6c 20 73 6f 0d |th too s|mall so.|
|00000fd0| 06 04 0d 42 20 64 6f 5f | 70 6c 6f 74 0d 06 0e 0a |...B do_|plot....|
|00000fe0| 2e 73 74 61 72 74 0d 06 | 18 21 a4 66 70 28 22 4d |.start..|.!.fp("M|
|00000ff0| 55 46 44 20 4d 2c 52 2c | 52 22 29 20 3b 72 2a 72 |UFD M,R,|R") ;r*r|
|00001000| 20 69 6e 74 6f 20 4d 0d | 06 22 38 a4 66 70 28 22 | into M.|."8.fp("|
|00001010| 4d 56 46 44 20 77 6b 31 | 2c 4d 22 29 3b 20 73 61 |MVFD wk1|,M"); sa|
|00001020| 76 65 20 72 2a 72 20 66 | 6f 72 20 75 73 65 20 69 |ve r*r f|or use i|
|00001030| 6e 20 61 20 6e 61 6e 6f | 2d 73 65 63 6f 6e 64 0d |n a nano|-second.|
|00001040| 06 2c 1c a4 66 70 28 22 | 4d 55 46 44 20 4d 2c 4d |.,..fp("|MUFD M,M|
|00001050| 2c 49 22 29 3b 20 72 2a | 72 2a 69 0d 06 36 1f a4 |,I"); r*|r*i..6..|
|00001060| 66 70 28 22 4d 55 46 44 | 20 4d 2c 4d 2c 23 32 22 |fp("MUFD| M,M,#2"|
|00001070| 29 3b 20 72 2a 72 2a 69 | 2a 32 0d 06 40 05 3b 0d |); r*r*i|*2..@.;.|
|00001080| 06 4a 25 a4 66 70 28 22 | 4d 55 46 44 20 77 6b 32 |.J%.fp("|MUFD wk2|
|00001090| 2c 49 2c 49 22 29 3b 20 | 69 2a 69 20 69 6e 74 6f |,I,I"); |i*i into|
|000010a0| 20 77 6b 32 0d 06 54 25 | a4 66 70 28 22 53 55 46 | wk2..T%|.fp("SUF|
|000010b0| 44 20 77 6b 32 2c 77 6b | 31 2c 77 6b 32 22 29 20 |D wk2,wk|1,wk2") |
|000010c0| 3b 20 72 2a 72 2d 69 2a | 69 0d 06 5e 19 a4 66 70 |; r*r-i*|i..^..fp|
|000010d0| 28 22 4d 55 46 44 20 77 | 6b 32 2c 77 6b 32 2c 49 |("MUFD w|k2,wk2,I|
|000010e0| 22 29 0d 06 68 33 a4 66 | 70 28 22 41 44 46 44 20 |")..h3.f|p("ADFD |
|000010f0| 4d 2c 4d 2c 77 6b 32 22 | 29 20 3b 20 32 2a 69 2a |M,M,wk2"|) ; 2*i*|
|00001100| 72 2a 72 2b 28 69 2a 28 | 72 2a 72 29 2d 28 69 2a |r*r+(i*(|r*r)-(i*|
|00001110| 69 29 29 3d 4d 0d 06 72 | 05 3b 0d 06 7c 21 a4 66 |i))=M..r|.;..|!.f|
|00001120| 70 28 22 4d 55 46 44 20 | 4e 2c 49 2c 49 22 29 20 |p("MUFD |N,I,I") |
|00001130| 3b 69 2a 69 20 69 6e 74 | 6f 20 4e 0d 06 86 1e a4 |;i*i int|o N.....|
|00001140| 66 70 28 22 4d 55 46 44 | 20 4e 2c 4e 2c 23 33 22 |fp("MUFD| N,N,#3"|
|00001150| 29 20 3b 20 69 2a 69 2a | 33 0d 06 90 22 a4 66 70 |) ; i*i*|3...".fp|
|00001160| 28 22 53 55 46 44 20 4e | 2c 77 6b 31 2c 4e 22 29 |("SUFD N|,wk1,N")|
|00001170| 3b 20 72 2a 72 2d 33 2a | 69 2a 69 0d 06 9a 15 a4 |; r*r-3*|i*i.....|
|00001180| 66 70 28 22 4d 55 46 44 | 20 4e 2c 4e 2c 52 22 29 |fp("MUFD| N,N,R")|
|00001190| 0d 06 a4 2a a4 66 70 28 | 22 53 55 46 44 20 4e 2c |...*.fp(|"SUFD N,|
|000011a0| 4e 2c 23 31 22 29 20 3b | 20 4e 3d 72 2a 28 72 2a |N,#1") ;| N=r*(r*|
|000011b0| 72 2d 33 2a 69 2a 69 29 | 2d 31 0d 06 ae 05 3b 0d |r-3*i*i)|-1....;.|
|000011c0| 06 b8 21 a4 66 70 28 22 | 4d 55 46 44 20 77 6b 31 |..!.fp("|MUFD wk1|
|000011d0| 2c 77 6b 31 2c 23 33 22 | 29 3b 20 33 2a 72 2a 72 |,wk1,#3"|); 3*r*r|
|000011e0| 0d 06 c2 1d a4 66 70 28 | 22 4d 55 46 44 20 77 6b |.....fp(|"MUFD wk|
|000011f0| 32 2c 49 2c 49 22 29 20 | 3b 20 69 2a 69 0d 06 cc |2,I,I") |; i*i...|
|00001200| 21 a4 66 70 28 22 4d 55 | 46 44 20 77 6b 32 2c 77 |!.fp("MU|FD wk2,w|
|00001210| 6b 32 2c 23 33 22 29 3b | 20 33 2a 69 2a 69 0d 06 |k2,#3");| 3*i*i..|
|00001220| d6 2a a4 66 70 28 22 53 | 55 46 44 20 77 6b 31 2c |.*.fp("S|UFD wk1,|
|00001230| 77 6b 31 2c 77 6b 32 22 | 29 3b 20 54 3d 33 2a 72 |wk1,wk2"|); T=3*r|
|00001240| 2a 72 2d 33 2a 69 2a 69 | 0d 06 e0 17 a4 66 70 28 |*r-3*i*i|.....fp(|
|00001250| 22 4d 55 46 44 20 77 6b | 32 2c 52 2c 49 22 29 0d |"MUFD wk|2,R,I").|
|00001260| 06 ea 1a a4 66 70 28 22 | 4d 55 46 44 20 77 6b 32 |....fp("|MUFD wk2|
|00001270| 2c 77 6b 32 2c 23 33 22 | 29 0d 06 f4 23 a4 66 70 |,wk2,#3"|)...#.fp|
|00001280| 28 22 4d 55 46 44 20 77 | 6b 32 2c 77 6b 32 2c 23 |("MUFD w|k2,wk2,#|
|00001290| 32 22 29 3b 20 55 3d 36 | 2a 69 2a 72 0d 06 fe 05 |2"); U=6|*i*r....|
|000012a0| 3b 0d 07 08 20 a4 66 70 | 28 22 4d 55 46 44 20 77 |;... .fp|("MUFD w|
|000012b0| 6b 33 2c 77 6b 31 2c 77 | 6b 31 22 29 3b 20 54 2a |k3,wk1,w|k1"); T*|
|000012c0| 54 0d 07 12 20 a4 66 70 | 28 22 4d 55 46 44 20 77 |T... .fp|("MUFD w|
|000012d0| 6b 34 2c 77 6b 32 2c 77 | 6b 32 22 29 3b 20 55 2a |k4,wk2,w|k2"); U*|
|000012e0| 55 0d 07 1c 20 a4 66 70 | 28 22 41 44 46 44 20 77 |U... .fp|("ADFD w|
|000012f0| 6b 33 2c 77 6b 33 2c 77 | 6b 34 22 29 3b 20 64 65 |k3,wk3,w|k4"); de|
|00001300| 6e 0d 07 26 2b a4 66 70 | 28 22 53 54 46 44 20 77 |n..&+.fp|("STFD w|
|00001310| 6b 33 2c 64 65 6e 22 29 | 3b 20 73 61 76 65 20 64 |k3,den")|; save d|
|00001320| 65 6e 20 66 6f 72 20 6c | 61 74 65 72 0d 07 30 05 |en for l|ater..0.|
|00001330| 3b 0d 07 3a 05 3b 0d 07 | 44 1e a4 66 70 28 22 4d |;..:.;..|D..fp("M|
|00001340| 55 46 44 20 77 6b 33 2c | 4e 2c 77 6b 31 22 29 3b |UFD wk3,|N,wk1");|
|00001350| 20 4e 2a 54 0d 07 4e 1e | a4 66 70 28 22 4d 55 46 | N*T..N.|.fp("MUF|
|00001360| 44 20 77 6b 34 2c 4d 2c | 77 6b 32 22 29 3b 20 4d |D wk4,M,|wk2"); M|
|00001370| 2a 55 0d 07 58 28 a4 66 | 70 28 22 41 44 46 44 20 |*U..X(.f|p("ADFD |
|00001380| 77 6b 33 2c 77 6b 33 2c | 77 6b 34 22 29 3b 20 77 |wk3,wk3,|wk4"); w|
|00001390| 6b 33 3d 4e 2a 54 2b 4d | 2a 55 0d 07 62 2f a4 66 |k3=N*T+M|*U..b/.f|
|000013a0| 70 28 22 53 54 46 44 20 | 77 6b 33 2c 6e 72 22 29 |p("STFD |wk3,nr")|
|000013b0| 3b 20 73 61 76 65 20 6e | 72 2c 20 6e 75 6d 65 72 |; save n|r, numer|
|000013c0| 61 74 6f 72 20 72 65 61 | 6c 0d 07 6c 05 3b 0d 07 |ator rea|l..l.;..|
|000013d0| 76 1e a4 66 70 28 22 4d | 55 46 44 20 77 6b 33 2c |v..fp("M|UFD wk3,|
|000013e0| 4e 2c 77 6b 32 22 29 3b | 20 4e 2a 55 0d 07 80 1e |N,wk2");| N*U....|
|000013f0| a4 66 70 28 22 4d 55 46 | 44 20 77 6b 34 2c 4d 2c |.fp("MUF|D wk4,M,|
|00001400| 77 6b 31 22 29 3b 20 4d | 2a 54 0d 07 8a 2e a4 66 |wk1"); M|*T.....f|
|00001410| 70 28 22 53 55 46 44 20 | 77 6b 34 2c 77 6b 34 2c |p("SUFD |wk4,wk4,|
|00001420| 77 6b 33 22 29 3b 20 54 | 2a 4d 2d 4e 2a 55 3d 6e |wk3"); T|*M-N*U=n|
|00001430| 69 20 69 6e 20 77 6b 34 | 0d 07 94 34 a4 66 70 28 |i in wk4|...4.fp(|
|00001440| 22 53 54 46 44 20 77 6b | 34 2c 6e 69 22 29 3b 20 |"STFD wk|4,ni"); |
|00001450| 73 61 76 65 20 6e 69 2c | 20 6e 75 6d 65 72 61 74 |save ni,| numerat|
|00001460| 6f 72 20 69 6d 61 67 69 | 6e 61 72 79 0d 07 9e 05 |or imagi|nary....|
|00001470| 3b 0d 07 a8 0c 2e 64 6f | 5f 63 61 6c 63 0d 07 b2 |;.....do|_calc...|
|00001480| 2b a4 66 70 28 22 4c 44 | 46 44 20 77 6b 33 2c 6e |+.fp("LD|FD wk3,n|
|00001490| 72 22 29 3b 20 67 65 74 | 20 6e 75 6d 65 72 61 74 |r"); get| numerat|
|000014a0| 6f 72 2c 20 72 65 61 6c | 0d 07 bc 29 a4 66 70 28 |or, real|...).fp(|
|000014b0| 22 4c 44 46 44 20 77 6b | 34 2c 64 65 6e 22 29 20 |"LDFD wk|4,den") |
|000014c0| 3b 20 67 65 74 20 64 65 | 6e 6f 6d 69 6e 61 74 6f |; get de|nominato|
|000014d0| 72 0d 07 c6 32 a4 66 70 | 28 22 44 56 46 44 20 77 |r...2.fp|("DVFD w|
|000014e0| 6b 31 2c 77 6b 33 2c 77 | 6b 34 22 29 20 3b 20 65 |k1,wk3,w|k4") ; e|
|000014f0| 72 72 6f 72 2c 20 72 65 | 61 6c 20 69 6e 74 6f 20 |rror, re|al into |
|00001500| 77 6b 31 0d 07 d0 20 a4 | 66 70 28 22 53 55 46 44 |wk1... .|fp("SUFD|
|00001510| 20 52 2c 52 2c 77 6b 31 | 22 29 3b 20 72 3d 72 2d | R,R,wk1|"); r=r-|
|00001520| 65 72 72 0d 07 da 05 3b | 0d 07 e4 2c a4 66 70 28 |err....;|...,.fp(|
|00001530| 22 4c 44 46 44 20 77 6b | 33 2c 6e 69 22 29 20 3b |"LDFD wk|3,ni") ;|
|00001540| 20 6e 75 6d 65 72 61 74 | 6f 72 20 69 6d 61 67 69 | numerat|or imagi|
|00001550| 6e 61 72 79 0d 07 ee 36 | a4 66 70 28 22 44 56 46 |nary...6|.fp("DVF|
|00001560| 44 20 77 6b 32 2c 77 6b | 33 2c 77 6b 34 22 29 20 |D wk2,wk|3,wk4") |
|00001570| 3b 20 65 72 72 6f 72 20 | 69 6d 61 67 69 6e 61 72 |; error |imaginar|
|00001580| 79 20 69 6e 74 6f 20 77 | 6b 32 0d 07 f8 20 a4 66 |y into w|k2... .f|
|00001590| 70 28 22 53 55 46 44 20 | 49 2c 49 2c 77 6b 32 22 |p("SUFD |I,I,wk2"|
|000015a0| 29 3b 20 69 3d 69 2d 65 | 72 69 0d 08 02 05 3b 0d |); i=i-e|ri....;.|
|000015b0| 08 0c 11 2e 63 6f 6d 70 | 61 72 65 5f 72 65 61 6c |....comp|are_real|
|000015c0| 0d 08 16 34 a4 66 70 28 | 22 41 42 53 44 20 77 6b |...4.fp(|"ABSD wk|
|000015d0| 31 2c 77 6b 31 22 29 20 | 3b 20 67 65 74 20 61 62 |1,wk1") |; get ab|
|000015e0| 73 6f 6c 75 74 65 20 6f | 66 20 72 65 61 6c 20 65 |solute o|f real e|
|000015f0| 72 72 6f 72 0d 08 20 1b | a4 66 70 28 22 4c 44 46 |rror.. .|.fp("LDF|
|00001600| 44 20 77 6b 33 2c 65 72 | 72 73 69 7a 65 22 29 0d |D wk3,er|rsize").|
|00001610| 08 2a 17 a4 66 70 28 22 | 43 4d 46 44 20 77 6b 31 |.*..fp("|CMFD wk1|
|00001620| 2c 77 6b 33 22 29 0d 08 | 34 2f 42 47 54 20 63 6f |,wk3")..|4/BGT co|
|00001630| 75 6e 74 20 3b 20 65 72 | 72 6f 72 20 72 65 61 6c |unt ; er|ror real|
|00001640| 20 68 61 73 6e 27 74 20 | 63 6f 6e 76 65 72 67 65 | hasn't |converge|
|00001650| 64 20 79 65 74 0d 08 3e | 05 3b 0d 08 48 2a 3b 65 |d yet..>|.;..H*;e|
|00001660| 72 72 6f 72 20 72 65 61 | 6c 20 77 61 73 20 4f 4b |rror rea|l was OK|
|00001670| 20 73 6f 20 6e 6f 77 20 | 64 6f 20 69 6d 61 67 69 | so now |do imagi|
|00001680| 6e 61 72 79 0d 08 52 11 | 2e 63 6f 6d 70 61 72 65 |nary..R.|.compare|
|00001690| 5f 69 6d 61 67 0d 08 5c | 1f a4 66 70 28 22 41 42 |_imag..\|..fp("AB|
|000016a0| 53 44 20 77 6b 32 2c 77 | 6b 32 22 29 20 3b 20 94 |SD wk2,w|k2") ; .|
|000016b0| 20 65 72 69 0d 08 66 17 | a4 66 70 28 22 43 4d 46 | eri..f.|.fp("CMF|
|000016c0| 44 20 77 6b 32 2c 77 6b | 33 22 29 0d 08 70 2f 42 |D wk2,wk|3")..p/B|
|000016d0| 47 54 20 63 6f 75 6e 74 | 20 3b 20 65 72 72 6f 72 |GT count| ; error|
|000016e0| 20 69 6d 61 67 20 68 61 | 73 6e 27 74 20 63 6f 6e | imag ha|sn't con|
|000016f0| 76 65 72 67 65 64 20 79 | 65 74 0d 08 7a 05 3b 0d |verged y|et..z.;.|
|00001700| 08 84 0c 2e 64 6f 5f 70 | 6c 6f 74 0d 08 8e 15 a4 |....do_p|lot.....|
|00001710| 66 70 28 22 4d 56 46 44 | 20 77 6b 31 2c 49 22 29 |fp("MVFD| wk1,I")|
|00001720| 0d 08 98 25 a4 66 70 28 | 22 41 42 53 44 20 77 6b |...%.fp(|"ABSD wk|
|00001730| 31 2c 77 6b 31 22 29 3b | 20 67 65 74 20 61 62 73 |1,wk1");| get abs|
|00001740| 6f 6c 75 74 65 0d 08 a2 | 1c 3b a4 66 70 28 22 4c |olute...|.;.fp("L|
|00001750| 44 46 44 20 77 6b 32 2c | 65 72 72 73 69 7a 65 22 |DFD wk2,|errsize"|
|00001760| 29 0d 08 ac 38 a4 66 70 | 28 22 4d 55 46 44 20 77 |)...8.fp|("MUFD w|
|00001770| 6b 33 2c 77 6b 33 2c 23 | 31 30 22 29 20 3b 20 74 |k3,wk3,#|10") ; t|
|00001780| 6f 20 67 65 74 20 69 74 | 20 6f 75 74 20 6f 66 20 |o get it| out of |
|00001790| 74 68 65 20 6e 6f 69 73 | 65 0d 08 b6 1b a4 66 70 |the nois|e.....fp|
|000017a0| 28 22 4d 55 46 44 20 77 | 6b 33 2c 77 6b 33 2c 23 |("MUFD w|k3,wk3,#|
|000017b0| 31 30 22 29 0d 08 c0 2a | a4 66 70 28 22 43 4d 46 |10")...*|.fp("CMF|
|000017c0| 44 20 77 6b 31 2c 77 6b | 33 22 29 20 3b 20 69 73 |D wk1,wk|3") ; is|
|000017d0| 20 49 20 63 6c 6f 73 65 | 20 74 6f 20 30 3f 0d 08 | I close| to 0?..|
|000017e0| ca 25 42 47 54 20 70 6c | 6f 74 5f 69 74 20 3b 20 |.%BGT pl|ot_it ; |
|000017f0| 6e 6f 20 73 6f 20 70 6c | 6f 74 20 74 68 65 20 72 |no so pl|ot the r|
|00001800| 6f 6f 74 0d 08 d4 31 2e | 73 65 74 5f 74 6f 5f 6e |oot...1.|set_to_n|
|00001810| 6f 75 67 68 74 20 3b 49 | 20 77 61 73 20 63 6c 6f |ought ;I| was clo|
|00001820| 73 65 20 74 6f 20 30 20 | 73 6f 20 6d 61 6b 65 20 |se to 0 |so make |
|00001830| 69 74 20 30 0d 08 de 16 | a4 66 70 28 22 4d 56 46 |it 0....|.fp("MVF|
|00001840| 44 20 49 2c 23 30 2e 30 | 22 29 0d 08 e8 0c 2e 70 |D I,#0.0|").....p|
|00001850| 6c 6f 74 5f 69 74 0d 08 | f2 16 a4 66 70 28 22 43 |lot_it..|...fp("C|
|00001860| 4d 46 44 20 49 2c 23 30 | 2e 30 22 29 0d 08 fc 38 |MFD I,#0|.0")...8|
|00001870| ec 51 20 52 30 2c 23 31 | 30 20 3b 20 63 6f 6c 6f |.Q R0,#1|0 ; colo|
|00001880| 75 72 73 20 66 6f 72 20 | 72 6f 6f 74 73 20 3b 74 |urs for |roots ;t|
|00001890| 68 69 73 20 6f 6e 65 20 | 69 66 20 69 28 72 6f 6f |his one |if i(roo|
|000018a0| 74 29 3d 30 0d 09 06 36 | 4d 4f 56 47 54 20 52 30 |t)=0...6|MOVGT R0|
|000018b0| 2c 23 32 30 20 3b 20 74 | 68 69 73 20 6f 6e 65 20 |,#20 ; t|his one |
|000018c0| 69 66 20 69 28 72 6f 6f | 74 29 20 3e 30 20 69 65 |if i(roo|t) >0 ie|
|000018d0| 20 30 2e 38 36 36 36 36 | 36 36 0d 09 10 3a 4d 4f | 0.86666|66...:MO|
|000018e0| 56 4c 54 20 52 30 2c 23 | 36 30 20 3b 20 61 6e 64 |VLT R0,#|60 ; and|
|000018f0| 20 74 68 69 73 20 6f 6e | 65 20 66 6f 72 20 69 28 | this on|e for i(|
|00001900| 72 6f 6f 74 29 3c 30 20 | 69 65 20 2d 30 2e 38 36 |root)<0 |ie -0.86|
|00001910| 36 36 36 36 0d 09 1a 32 | 80 20 52 33 2c 52 33 2c |6666...2|. R3,R3,|
|00001920| 23 33 20 3b 20 67 65 74 | 20 52 33 20 69 6e 74 6f |#3 ; get| R3 into|
|00001930| 20 72 61 6e 67 65 20 30 | 20 74 6f 20 33 20 66 6f | range 0| to 3 fo|
|00001940| 72 20 74 69 6e 74 0d 09 | 24 23 4d 4f 56 20 52 30 |r tint..|$#MOV R0|
|00001950| 2c 52 30 2c 4c 53 4c 23 | 32 3b 20 73 65 74 20 74 |,R0,LSL#|2; set t|
|00001960| 68 65 20 63 6f 6c 6f 75 | 72 0d 09 2e 1f 41 44 44 |he colou|r....ADD|
|00001970| 20 52 30 2c 52 30 2c 52 | 33 20 3b 20 73 65 74 20 | R0,R0,R|3 ; set |
|00001980| 74 68 65 20 74 69 6e 74 | 0d 09 38 05 3b 0d 09 42 |the tint|..8.;..B|
|00001990| 23 3b 20 52 65 73 74 6f | 72 65 20 46 6c 6f 61 74 |#; Resto|re Float|
|000019a0| 20 72 65 67 69 73 74 65 | 72 73 20 46 34 2d 46 37 | registe|rs F4-F7|
|000019b0| 0d 09 4c 1a a4 66 70 28 | 22 4c 44 46 44 20 46 34 |..L..fp(|"LDFD F4|
|000019c0| 2c 66 6c 6f 61 74 66 34 | 22 29 0d 09 56 1a a4 66 |,floatf4|")..V..f|
|000019d0| 70 28 22 4c 44 46 44 20 | 46 35 2c 66 6c 6f 61 74 |p("LDFD |F5,float|
|000019e0| 66 35 22 29 0d 09 60 1a | a4 66 70 28 22 4c 44 46 |f5")..`.|.fp("LDF|
|000019f0| 44 20 46 36 2c 66 6c 6f | 61 74 66 36 22 29 0d 09 |D F6,flo|atf6")..|
|00001a00| 6a 1a a4 66 70 28 22 4c | 44 46 44 20 46 37 2c 66 |j..fp("L|DFD F7,f|
|00001a10| 6c 6f 61 74 66 37 22 29 | 0d 09 74 05 3b 0d 09 7e |loatf7")|..t.;..~|
|00001a20| 0f 4d 4f 56 53 20 50 43 | 2c 52 31 34 0d 09 88 05 |.MOVS PC|,R14....|
|00001a30| 3b 0d 09 92 0a 2e 63 6f | 75 6e 74 0d 09 9c 10 41 |;.....co|unt....A|
|00001a40| 44 44 20 52 33 2c 52 33 | 2c 23 31 0d 09 a6 0b 42 |DD R3,R3|,#1....B|
|00001a50| 20 73 74 61 72 74 0d 09 | b0 0f a4 66 69 6e 74 65 | start..|...finte|
|00001a60| 72 66 61 63 65 0d 09 ba | 06 5d ed 0d 09 c4 38 c8 |rface...|.]....8.|
|00001a70| 99 20 22 4f 53 5f 46 69 | 6c 65 22 2c 31 30 2c 22 |. "OS_Fi|le",10,"|
|00001a80| 3c 46 72 61 63 24 44 69 | 72 3e 2e 43 75 62 65 52 |<Frac$Di|r>.CubeR|
|00001a90| 6f 6f 74 22 2c 26 46 46 | 41 2c 2c 63 6f 64 65 25 |oot",&FF|A,,code%|
|00001aa0| 2c 4f 25 0d 09 ce 05 e0 | 0d ff                   |,O%.....|..      |
+--------+-------------------------+-------------------------+--------+--------+