ZXTape! 2 Halley's Comet Andrew Giles Andrew Giles1983English EducationalNoneOriginal%TZXed by Andrew Barker For Karl BrownComet1 %OdF'Od) :6:0:0:7: n1,5;"Please select an item":5,6;"Press K for Kepler":8,6;"Press A for Altaz":11 ,6;"Press M for Moon":14,6;"Press H for Halley" x%17,6;"Press V to Quit" ="K"Ŧ="k"1000 ="a"Ŧ="A"1600@ ="m"Ŧ="M"4000 ="H"Ŧ="h"4500 ="v"Ŧ="V": 130   :"KEPLER"  1 (2,6;"ORBIT DEMONSTRATION"  0 C"Semimajor axis (1-20)",a$:a$<480ůa$>5791008 a=a$ "a<1a>201008 a>10 a=10 C"Eccentricity (0-0.95)",e$:e$<480ůe$>5791016 e=e$ $e<0e>0.95s3331016 O"Time of perihelion passage(year) ",p$:p$<480ůp$>5791022 pp=p$ A"Initial time (year)",s$:s$<480ůs$>5791026 st=s$ J"Yearly motion (degrees/year)",y$:y$<480ůy$>5791030 0ym=y$:ym<0.1}Lym>360 h1030  E"Number of dates(10-100)",u$:u$<480ůu$>5791034   u=u$ $u<10 u>100d1034  L"Time interval (years e.g. 0.5)",t$:t$<480ůt$>5791040 -ti=t$:ti<0.1}Lti>51040 $18,8;"Please wait..." *k(100d):l(100d):r(100d) x(100d):y(100d) q=a*(1-e) b=a*(1-e*e) 1:6  useful constants "RAD=/180 $DEG=180/ &thetae=0 ( i=1u * 1160 ,&thetac=ym*(st+(i-1)*ti-pp)*RAD .next approximation 0 diff=thetac-(thetae-e*thetae) 2(thetae=thetae+diff/(1-e*thetae) 4"(diff)>0.0001sQX1070. 6convert polar to cartesian 8x(i)=a*(thetae-e) :y(i)=b*thetae <i >scale for screen @ i=1u B 1160 Dk(i)=30-q+x(i) Fl(i)=11 +y(i) Hr(i)=(x(i)*x(i)+y(i)*y(i)) Ji L NJ210,83S:4,0:212,81Q:0,4 P plot points R i=1u T&8*k(i),8*(21-l(i)) VW1,1:-2,0:0,-2:2,0:0,2 War(i)5791740 l=l$ %l>90 Zl<-90Z1740  l=RAD*l  "Enter year (1985 or 1986)",y$ !y$"1985"y$"1986"1770 y$="1985"b=17.359625  y$="1986"b=17.375540  J"Enter month (1=Jan...12=Dec)",n$:n$<480ůn$>5791790 n=n$ "n<1n>12 1790  G"Enter day of month (1-31)",d$:d$<480ůd$>5791805  d=d$ "d<1d>311805   !?"Enter hour (0-23)",h$:h$<480ůh$>5791825! &h=h$ +"h<0h>231825! 0 5A"Enter minute (0-59)",m$:m$<480ům$>57918455 :m=m$ ?"m<0m>59;18455 D *FOR SUN* I 2720 N8y$="1985"x=3653E+m(n)+d+(h+m/60<)/24 S8y$="1986"x=4018+m(n)+d+(h+m/60<)/24 Xv=360h*x/365.256 ]#v0v360h1905q bv=v-360h g#v0v360h1905q l 1890b qy=v+g(2)-w(2) v 2720 {y<0y=y+360h  y=RAD*y e=(360h*e(3)*y)/ l(2)=v+e+g(2) 7l(2)>360hl(2)=l(2)-360h l(2)=RAD*l(2) d(2)=(o*l(2))  2720 'r(2)=(l(2)/d(2)) /o*l(2)>0r(2)=r(2) 9o*l(2)<0r(2)=2*-r(2)  2720  1:  &0,8;"DATA FOR THE SUN:"  8" Right Ascension:";(17);DEG*r(2);" deg."  2720 4" Declination:";(17);DEG*d(2);" deg." flag=12060  x(5)=h+m/60< n(5)=m(n)+d (p(5)=n(5)*a-b+x(5)*c 5p(5)>24p(5)=p(5)-24 4p(5)<0p(5)=p(5)+24 1h(2)=RAD*(p(5)*15)-r(2) .h(5)=h(2):d(5)=d(2)  2620< \v(2)=v(5):u(2)=u(5):k(2)=v(2):q(2)=u(2)  0" Azimuth:";(17);DEG*v(2);" deg." 1" Altitude:";(17);DEG*u(2);" deg."   2720  DDEG*u(2)>-181:7,9 ;"sky not dark" %0:0:1 *flag02380L / *FOR COMET* 47f=360h*(x-4058.66}])/(365.256*76L) 9f<0f=f+360h >#f0f360h2130R Cf=f-360h H#f0f360h2130R M 2115C R f=RAD*f Wq=f \ 2720 af(1)=q-e(1)*q ft=(f-f(1)) kt<0.0001sQX2175 p3c(1)=(f-f(1))/(1-e(1)*q) uq=q+c(1) z 2145a =g=((1+e(1))/(1-e(1)))*(q/2) .n(1)=g:n(1)=2*n(1) Qw(1)=RAD*w(1):z(1)=RAD*z(1):i(1)=RAD*i(1)  2720 "l(1)=n(1)+w(1) Qr=(a(1)*(1-e(1)*e(1)))/(1+e(1)*n(1)) "g(1)=l(1)-z(1) 'p(1)=(g(1)*i(1)) 2l(4)=(g(1)*i(1))+z(1) Lg(1)ȧ/2g(1)3*/2l(4)=l(4)+ 5l(4)<0l(4)=l(4)+2* s=r*p(1)  *FOR EARTH* .v=(360h*x)/(365.256*1.00004O)  2720 #v0v360h2270 v=v-360h #v0v360h2270  2255 $y(3)=v+g(3)-w(3) y(3)=RAD*y(3) >l(3)=v+(360h*e(3)*y(3))/+g(3) 5l(3)<0l(3)=l(3)+360h 7l(3)>360hl(3)=l(3)-360h l(3)=RAD*l(3)  2720 w(3)=RAD*w(3) "n(3)=l(3)-w(3) Ms(3)=(1-e(3)*e(3))/(1+e(3)*n(3)) dif=l(3)-l(4) Kss(3)l(5)=((s(3)*(-dif))/(s-s(3)*(-dif)))+l(4) 5l(5)<0l(5)=l(5)+2* $^b(3)=((s*p(1)*(l(5)-l(4)))/(s(3)*(l(4)-l(3)))) )9d(1)=(b(3)*o+b(3)*o*l(5)) .3r(1)=(b(3)*l(5)/d(1)) 3 2720 8K-o*b(3)+o*b(3)*l(5)>0r(1)=r(1) =U-o*b(3)+o*b(3)*l(5)<0r(1)=2*-r(1) B G L-9 ,6;"DATA FOR HALLEY'S COMET:" Q V8" Right ascension:";(17);DEG*r(1);" deg." [4" Declination:";(17);DEG*d(1);" deg." ` 2720 eflag02440 j1h(1)=RAD*(p(5)*15)-r(1) o.h(5)=h(1):d(5)=d(1) t 2620< y\v(1)=v(5):u(1)=u(5):k(1)=v(1):q(1)=u(1) ~ 2720 0" Azimuth:";(17);DEG*v(1);" deg." 1" Altitude:";(17);DEG*u(1);" deg."  #" Record values then press ENTER" #" for a pictorial representation" #" of the position of the comet.." =""2470 2720 *comet graphics* ?t(1)=((q(1)-q(2))/(k(2)-k(1))) Fk(2)q(2)t(1)=t(1)+ Fk(2)q(1)t(1)=t(1)+ 8c=(10 /r)*t(1):d=(10 /r)*t(1) 2720 0 7:0  $0,0:175,0 $0,0:0,175 l(5)=k(1) 2595# =k(4)=1.8ffff*DEG*(k(1)-(j-1)*/2) 2720 %q(4)=1.8ffff*DEG*q(1) 2650Z 30,10 ;"Press ENTER":=""2565 6:0:7:1::4,9 ;"Select an item":7,6;"Press key R to review":10 ,6;"Press key U to re-run":13 ,6;"Press key Q for menu" '="R"Ŧ="r"flag=1:1975 ="U"Ŧ="u":1600@  ="Q"Ŧ="q"100d  2575 #4l(5)0l(5)ǧ/2j=1 (.l(5)ȧ/2l(5)ǧj=2 -6l(5)ȧl(5)3*/2j=3 2>l(5)3*/2l(5)2*j=4 7 <7h(5)>2*h(5)=h(5)-2* A5h(5)<0h(5)=h(5)+2* F9u(5)=(d(5)*l+d(5)*l*h(5)) K=v(5)=((d(5)-l*u(5))/(l*u(5))) P6h(5)>0v(5)=2*-v(5) U Z|k(4),q(4):k(4)+1,q(4):k(4),q(4)-1:k(4)+1,q(4)-1 _&0,0;579;480 d(j=121,0;78N i(j=221,0;69E n(j=321,0;83S s(j=421,0;87W x*j=121,21;69E, })j=221,21;83S )j=321,21;87W )j=421,21;78N k(4),q(4) -c,d  END OF PROGRAM ="q"Ŧ="Q":00100d   :"MOON" 6:7:1: )0:1,10 ;"MOON PHASE"  r(19)  l(33!)  m(12 )  4016 i=119 r(i) 9 ,25,10 ,26,12 ,28,13 ,29,15,1,17,3,19,4,20,6,22,7,23 i i=133! l(i) 9 ,28,17,6,25,14,3,22,11 ,30,19,8,27,16,5,24,13 ,2,21,10 ,29 i18,7,26,15,4,23,12 ,1,20,9 ,28,17 i i=112 m(i) d0,2,2,4,4,6,7,8,9 ,10 ,11 ,13 i  4"Year",y$:y$<48 0ůy$>57 94040 y=y$:y<04040 yr=y "BC or AD",a$ a$="ad"a$="AD" a$="bc"a$="BC" a$"AD"a$"BC"4046 Q"Enter month: (1=Jan, 2=Feb,12=Dec)",z$:z$<480ůz$>57 94054 )z=z$:z<1z>12 4054 /z=1z=2z=z+10 :4062  z=z-2 $"Enter J(ulian) or G(regorian)",b$ b$="j"b$="J" b$="g"b$="G" b$"J"b$"G"4062 M"Enter phase: 0(new) to 15(full)",c$:c$<480ůc$>5794070 c=c$ "c<0c>304070 z>10 y=y-1 a$="AD"j=y+4712h a$="BC"j=4713i-y q=(j/76L) n=j-q*76L s=(q/4) r=(n/4) l=n-r*4 i=130 l(i)=r(r+1)d=l(i+l) i w=d-s-m(z)+c  b$="J"t=w 'y>1582.y1700y=10 'y>1700y1800y=11  'y>1800y1900ly=12  y>1900ly=13 b$="G"t=w+y t<0t=t+30 t>30t=t-30  &z>10 z=z-10 :4122  z=z+2 46,3;"A moon ";c;" days old occurs on" <9 ,3;"Day ";t;", month ";z;", year ";yr;" ";a$ )1:14,12 ;"Hit SPACE"   " "4128  " $`0:6,5;"Press key C to continue":10 ,5;"Press key Q for menu" &="q"100d ( ="c":4000 ) 4134& * 100d , 00:6:1:7::"HALLEY" "1,8;"DEMONSTRATION"  'a(40():b(40():c(40() 'd(40():h(40():m(40() n(40():r(40() x(40():y(40() I"Gaussian constant (AU/year)",k$:k$>579ůk$<4804535 +k=k$:k<0.1}Lk>254535  4840 M"Perihelion distance (AU)",q$:q$>579ůq$<4804550:q=q$ q=q$ #q<0.1}Lq>44550 H"Year of perihelion passage",p$:p$<480ůp$>5794570 p=p$ P"Month of perihelion passage",m$:m$<480ům$>5794580:m=m$ m=m$ %m<1m>12.99O =4590 p=(p)+(m-1)/12 B"Year of initial date",s$:s$<480ůs$>5794600 s=s$ C"Month of initial date",m$:m$<480ům$>5794610 m=m$  %m<1m>12.99O =4610 s=(s)+(m-1)/12 D"Time interval (months)",w$:w$<480ůw$>5794630  w=w$:w<0.1}L4630  w=w/12 %E"number of dates (5-40)",u$:u$<48 0ůu$>5794645% *u=u$ /"u<5u>40(4640  4#5,9 ;"Please wait..." 9 i=1u > 4840 C3m(i)=k*(s+(i-1)*w-p)/((2*q^3)) Hi=1a(i)=m(i) MNh(i)=-(a(i)*a(i)*a(i)+3*a(i)-3*m(i))/(3*a(i)*a(i)+3) Rb(i)=a(i)+h(i) W/c(i)=b(i)*b(i)*b(i)+3*b(i)-3*m(i) \ d(i)=c(i) ad(i)<0.0001sQX4720p f a(i)=b(i) k 4685M p n(i)=b(i) ui z i=1u  4840 r(i)=q*(1+n(i)*n(i)) n(i)=2*n(i) x(i)=4*r(i)*n(i) y(i)=4*r(i)*n(i) x(i)=20-q+x(i) y(i)=10 +y(i) i 0:0:6  :10 ,20;42*  i=1u  4840 0:6 *(y(i)+0.5),(x(i)+0.5);"" 0:7 ^10 -4*(((i-2)*/6)),20+4*(((i-2)*/6));"o"  502 f10 -4*(((i-2)*/6)),20+4*(((i-2)*/6));(32 ) i $1:6:0:7 421,10 ;"Press ENTER":=""4831 :2,9 ;"Select an item":5,6;"Press key R to review":8,6;"Press key U to re-run":11 ,6;"Press Q for menu" ="R"Ŧ="r"4775 ="U"Ŧ="u":4500 ="q"Ŧ="Q"100d  4833 ="q"100d  'addr="64967"::ٰ"9";"SPTK 48K"'" 1983 Andrew Giles"'''"Select:"'"1 Renumber"'"2 Delete"'"3 Fremem"'"4 Vardump"'"5 Search & Replace":"23560","" 'c"23560"="1""Number of first line ";tool:"9998":"Line increment ";tool:"9998":"64059" 'T"23560"="2""First line ";tool:"9998":"Last line ";tool:"9998":"64884" ' ."23560"="3":"Free memory=";"64981": ' "23560"="4""64995": ' "23560"ɯ"5"9991' ' maddr="23296":"String to replace ";e$:"9999":"To replace with ";e$:e$tool"Must be same length": ' ."9999":addr="64965":"9998":"65230": 'Jaddr,"tool-256*(tool/256)":addr+"1","(tool/256)":addr=addr-"2": '8e="1"̱e$:addr,e$(e):addr=addr+"1":e:tool=e$: ((((((((((K6kQ4qP3pEUUUmS2sUUUM2W2w~*U6ykepler   "KEPLER" "ORBIT DEMONSTRATION" "Semimajor axis (1-20)",a ( a<1a>20502 2"Eccentricity (0-0.95)",e <"e<0e>0.95s333502 F("Time of perihelion passage(year) ",pp P"Initial time (year)",st Z#"Yearly motion (degrees/year)",ym d"Number of dates(10-100)",u n#u<10 u>100d100d x%"Time interval (years e.g. 0.5)",ti "Please wait..." *k(100d):l(100d):r(100d) x(100d):y(100d) q=a*(1-e) b=a*(1-e*e) 1:6 useful constants RAD=/180 DEG=180/ thetae=0 i=1u &thetac=ym*(st+(i-1)*ti-pp)*RAD next approximation  diff=thetac-(thetae-e*thetae) (thetae=thetae+diff/(1-e*thetae) !(diff)>0.0001sQX250 "convert polar to cartesian ,x(i)=a*(thetae-e) 6y(i)=b*thetae @i Jscale for screen T i=1u ^k(i)=30-q+x(i) hl(i)=11 +y(i) rr(i)=(x(i)*x(i)+y(i)*y(i)) |i  J210,83S:4,0:212,81Q:0,4  plot points  i=1u &8*k(i),8*(21-l(i)) W1,1:-2,0:0,-2:2,0:0,2 ar(i)90 Zl<-90Z280 , l=RAD*l 6 @ "Enter year (1985 or 2986)",y$ J y$"1985"y$"1986"320@ Ty$="1985"b=17.359625  ^y$="1986"b=17.375540  h r""Enter month (1=Jan...12=Dec)",n |!n<1n>12 370r  "Enter day of month (1-31)",d !d<1d>31400  "Enter hour (0-23)",h !h<0h>23430  "Enter minute (0-59)",m !m<0m>59;460  *FOR SUN* 8y$="1985"x=3653E+m(n)+d+(h+m/60<)/24 8y$="1986"x=4018+m(n)+d+(h+m/60<)/24 v=360h*x/365.256 "v0v360h5600 v=v-360h "v0v360h5600 & 530 0y=v+g(2)-w(2) :y<0y=y+360h D y=RAD*y Ne=(360h*e(3)*y)/ Xl(2)=v+e+g(2) b7l(2)>360hl(2)=l(2)-360h ll(2)=RAD*l(2) vd(2)=(o*l(2)) 'r(2)=(l(2)/d(2)) /o*l(2)>0r(2)=r(2) 9o*l(2)<0r(2)=2*-r(2)  "DATA FOR THE SUN:"  7"Right Ascension:";(16);DEG*r(2);" deg." 3"Declination:";(16);DEG*d(2);" deg." x(5)=h+m/60< n(5)=m(n)+d (p(5)=n(5)*a-b+x(5)*c 5p(5)>24p(5)=p(5)-24 4p(5)<0p(5)=p(5)+24 1h(2)=RAD*(p(5)*15)-r(2)  .h(5)=h(2):d(5)=d(2)  1760  \v(2)=v(5):u(2)=u(5):k(2)=v(2):q(2)=u(2) * 4/"Azimuth:";(16);DEG*v(2);" deg." >0"Altitude:";(16);DEG*u(2);" deg." H R \*DEG*u(2)>-18"sky not dark" f *FOR COMET* p7f=360h*(x-4058.66}])/(365.256*76L) zf<0f=f+360h "f0f360h940 f=f-360h "f0f360h940  910  f=RAD*f q=f f(1)=q-e(1)*q t=(f-f(1)) t<0.0001sQX1020 3c(1)=(f-f(1))/(1-e(1)*q) q=q+c(1)  960 =g=((1+e(1))/(1-e(1)))*(q/2) .n(1)=g:n(1)=2*n(1) Qw(1)=RAD*w(1):z(1)=RAD*z(1):i(1)=RAD*i(1) "l(1)=n(1)+w(1) $Qr=(a(1)*(1-e(1)*e(1)))/(1+e(1)*n(1)) ."g(1)=l(1)-z(1) 8'p(1)=(g(1)*i(1)) B2l(4)=(g(1)*i(1))+z(1) LLg(1)ȧ/2g(1)3*/2l(4)=l(4)+ V5l(4)<0l(4)=l(4)+2* `s=r*p(1) j *FOR EARTH* t.v=(360h*x)/(365.256*1.00004O) ~#v0v360h1190 v=v-360h #v0v360h1190  1160 $y(3)=v+g(3)-w(3) y(3)=RAD*y(3) >l(3)=v+(360h*e(3)*y(3))/+g(3) 5l(3)<0l(3)=l(3)+360h 7l(3)>360hl(3)=l(3)-360h l(3)=RAD*l(3) w(3)=RAD*w(3) "n(3)=l(3)-w(3) Ms(3)=(1-e(3)*e(3))/(1+e(3)*n(3)) dif=l(3)-l(4)  Kss(3)l(5)=((s(3)*(-dif))/(s-s(3)*(-dif)))+l(4) 5l(5)<0l(5)=l(5)+2* (^b(3)=((s*p(1)*(l(5)-l(4)))/(s(3)*(l(4)-l(3)))) 29d(1)=(b(3)*o+b(3)*o*l(5)) <3r(1)=(b(3)*l(5)/d(1)) FK-o*b(3)+o*b(3)*l(5)>0r(1)=r(1) PU-o*b(3)+o*b(3)*l(5)<0r(1)=2*-r(1) Z d n"DATA FOR HALLEY'S COMET:" x 7"Right ascension:";(16);DEG*r(1);" deg." 3"Declination:";(16);DEG*d(1);" deg." 1h(1)=RAD*(p(5)*15)-r(1) .h(5)=h(1):d(5)=d(1)  1760 \v(1)=v(5):u(1)=u(5):k(1)=v(1):q(1)=u(1)  /"Azimuth:";(16);DEG*v(1);" deg." 0"Altitude:";(16);DEG*u(1);" deg."  !"Record values then press CONT" ""for a pictorial representation" #"of the position of the comet..."  *comet graphics* ?t(1)=((q(1)-q(2))/(k(2)-k(1))) "Fk(2)q(2)t(1)=t(1)+ ,Fk(2)q(1)t(1)=t(1)+ 68c=(10 /r)*t(1):d=(10 /r)*t(1) @ J 0 T7:0 ^$0,0:175,0 h$0,0:0,175 rl(5)=k(1) | 1710 =k(4)=1.8ffff*DEG*(k(1)-(j-1)*/2) %q(4)=1.8ffff*DEG*q(1)  1820  1950 4l(5)0l(5)ǧ/2j=1 .l(5)ȧ/2l(5)ǧj=2 6l(5)ȧl(5)3*/2j=3 >l(5)3*/2l(5)2*j=4  7h(5)>2*h(5)=h(5)-2* 5h(5)<0h(5)=h(5)+2* 9u(5)=(d(5)*l+d(5)*l*h(5)) =v(5)=((d(5)-l*u(5))/(l*u(5))) 6h(5)>0v(5)=2*-v(5)  |k(4),q(4):k(4)+1,q(4):k(4),q(4)-1:k(4)+1,q(4)-1 &&0,0;579;480 0(j=121,0;78N :(j=221,0;69E D(j=321,0;83S N(j=421,0;87W X*j=121,21;69E, b)j=221,21;83S l)j=321,21;87W v)j=421,21;78N k(4),q(4) -c,d  END OF PROGRAM |nw:\uG3|\M3wO{aH~w( NG _C=458 YAGZvR|~#f̈́ Y? ;Zx0NIЅ|8y6HIgNvR E8uy6Hy6 Y Y=s ATe"""L@a{5ee. oQv;a}pcАl_fY1985b n dF9rhmxu vƮyge .fHq:5tkggƓyrFsPޤi恁#jmoon  4w "MOON"  "MOON PHASE" ( r(19) 2 l(33!) < m(12 ) Fi=119 Pr(i) Z9 ,25,10 ,26,12 ,28,13 ,29,15,1,17,3,19,4,20,6,22,7,23 di ni=133! xl(i) 9 ,28,17,6,25,14,3,22,11 ,30,19,8,27,16,5,24,13 ,2,21,10 ,29 i18,7,26,15,4,23,12 ,1,20,9 ,28,17 i i=112 m(i) d0,2,2,4,4,6,7,8,9 ,10 ,11 ,13 i  "Year",y y<0210 yr=y "BC or AD",a$ a$="ad"a$="AD" a$="bc"a$="BC" a$"AD"a$"BC"240 5"Enter month: (11=Jan, 12=Feb, 1=Mar...10=Dec)",z "!z<1z>12 280 ,$"Enter J(ulian) or G(regorian)",b$ 6b$="j"b$="J" @b$="g"b$="G" Jb$"J"b$"G"300, T%"Enter phase: 0(new) to 15(full)",c ^!c<0c>30340T hz>10 y=y-1 ra$="AD"j=y+4712h |a$="BC"j=4713i-y q=(j/76L) n=j-q*76L s=(q/4) r=(n/4) l=n-r*4 i=130 l(i)=r(r+1)d=l(i+l) i w=d-s-m(z)+c  b$="J"t=w 'y>1582.y1700y=10 'y>1700y1800y=11 'y>1800y1900ly=12 y>1900ly=13 b$="G"t=w+y t<0t=t+30 &t>30t=t-30 0 :#"A moon ";c;" days old occurs on" D+"Day ";t;", month ";z;", year ";yr;" ";a$ N 200 b    !      ?  y AADzBJcj(qXnsrld wthalley  h4  "HALLEY" "HALLEY, EARTH, SUN DEMO"  ('a(40():b(40():c(40() 2'd(40():h(40():m(40() <n(40():r(40() Fx(40():y(40() P!"Gaussian constant (AU/year)",k Z"Perihelion distance (AU)",q d n "Year of perihelion passage",p x!"Month of perihelion passage",m $m<1m>12.99O =130 p=(p)+(m-1)/12 "Year of initial date",s "Month of initial date",m $m<1m>12.99O =160 s=(s)+(m-1)/12 "Time interval (months)",w w=w/12 "number of dates (5-40)",u !u<5u>40(200 "Please wait..." i=1u 3m(i)=k*(s+(i-1)*w-p)/((2*q^3)) i=1a(i)=m(i) Nh(i)=-(a(i)*a(i)*a(i)+3*a(i)-3*m(i))/(3*a(i)*a(i)+3) b(i)=a(i)+h(i) "/c(i)=b(i)*b(i)*b(i)+3*b(i)-3*m(i) , d(i)=c(i) 6d(i)<0.0001sQX340T @ a(i)=b(i) J 270 T n(i)=b(i) ^i h i=1u rr(i)=q*(1+n(i)*n(i)) |n(i)=2*n(i) x(i)=4*r(i)*n(i) y(i)=4*r(i)*n(i) x(i)=20-q+x(i) y(i)=10 +y(i) i 0:5  :10 ,20;42*  i=1u 0:6 0(y(i)+0.5),(x(i)+0.5);44, 0:2 d10 -4*(((i-2)*/6)),20+4*(((i-2)*/6));43+  502 f10 -4*(((i-2)*/6)),20+4*(((i-2)*/6));(32 ) i (M [iR&A.}2"3֘))v)|V 3("؁>{?R.}2bq<&RH [gL(ܘA"g7L _=%iѪ8rs@R֍uVUV́ @Q(X>r=Q7bhq8S _7 "em(gfsqHnkgs:oirtikkd@irti os'6g k n2@qHs&(gfs&qHn2kgs':oirti kckd@irti os'6g k n2@qHs&(r4qNxdmwXv.4.t`yIrQZytw]etpgxV"t'pvM4vLtox#towy-rYPtyF'v4.wتEx|k(#ƗN4 ~~zd簊 ƗN4~ǂdƗN4dZdZddևFN4dև%~]FN4g}dܻz%~]6 I(9Oc ނ;?ڂtɊէꁲKf;k)e;)2:EUk܁z:u?g :i ނ剂9('h 6 `U; 0&,iSE6<4f$ cR.Sg|݀L̀|݁.S؁cR c$g<4kS