home *** CD-ROM | disk | FTP | other *** search
/ Loadstar 163 / 163.d81 / expanding (.txt) < prev    next >
Commodore BASIC  |  2022-08-26  |  3KB  |  99 lines
  1. 5 poke55,.:poke56,56:clr:poke53280,.:poke53281,.:print"[147]"
  2. 10 dv=peek(186):ifdv<8thendv=8
  3. 30 poke53272,31:poke53371,0
  4. 40 ad=49152
  5. 45 sysad:sysad+12
  6. 50 sysad+9,0
  7. 55 print"[147]"
  8. 90 bs$="[164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164][157][164]"
  9. 95 print"[147]":sysad+9,1
  10. 100 rm=(3700/2.7)-1
  11. 105 print"[158][220][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][221]"
  12. 110 print"[158]"bs$""tab(38)bs$
  13. 115 print"[158][255][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][161]"
  14. 120 print""tab(11)"-[150][197][216][208][193][206][196][201][206][199] [195][207][211][205][207][211]-"
  15. 122 print""tab(8)"[154]:[150][211]etting up [208]arameters[154]:"
  16. 125 printtab(2)"[156][195]osmos' ([193][156])ge/([200][156])ubble parameter?":poke198,.
  17. 130 geta$:ifa$<>"a"anda$<>"h"then130
  18. 135 sysad+9,2
  19. 140 ifa$="h"then160
  20. 145 printtab(2)"[158][193]ge of [195]osmos [156][206][207][215] [158]in bil/yrs: ";:l9%=5:gosub480:tn=q9
  21. 150 sysad+9,2
  22. 155 tn=tn*1e+09:goto175
  23. 160 printtab(2)"[158][200]ubble parameter [156](km/sec/[205]pc): ";:l9%=4:gosub480:hn=q9
  24. 165 sysad+9,2
  25. 170 tn=(2/3)/(hn/9.81e+11)
  26. 175 print"[156]([212][156])ime [212][200][197][206]/([210][156])edshift [206][207][215]?":poke198,.
  27. 180 geta$:ifa$<>"t"anda$<>"r"then180
  28. 185 sysad+9,2
  29. 190 ifa$="r"then230
  30. 195 printtab(2)"[158][193]ge of [195]osmos [156][212][200][197][206] [158]in bil/yrs: ";:l9%=4:gosub480:tt=q9
  31. 200 iftt=<.03thenprint"[145][145][145]":goto195
  32. 201 sysad+9,2
  33. 202 print:printtab(9)"[159][201]s this [195]orrect? (y[159]/n[159])":poke198,.
  34. 203 gosub525
  35. 204 ifa$="n"then95
  36. 205 sysad+9,2
  37. 210 tt=tt*1e+09:sc=(tn/tt)^(2/3):rs=sc-1:ifrs<rmthen255
  38. 215 sm=rm+1:tm=tn/(sm^1.5)
  39. 220 print"[159][195]an't look back in time before":printtab(4)tm
  40. 225 sysad+9,2
  41. 230 printtab(2)"[158][210]edshift value for light [156][206][207][215]: ";:l9%=5:gosub480:rs=q9
  42. 235 sysad+9,2
  43. 237 print:printtab(9)"[159][201]s this [195]orrect? (y[159]/n[159])":poke198,.
  44. 238 gosub525
  45. 239 sysad+9,2
  46. 240 sc=rs+1:tt=tn/(sc^1.5):ifrs<rmthen255
  47. 245 printtab(2)"[159][210]edshift value too large before":printtab(4)rm
  48. 250 goto215
  49. 255 print"[147]":sysad+9,3
  50. 260 print"[150][220][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][221]"
  51. 265 printbs$""tab(38)bs$
  52. 270 print"[150][255][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][162][161]"
  53. 275 print""tab(12)"[159][197][216][208][193][206][196][201][206][199] [195][207][211][205][207][211]"
  54. 280 print:printtab(2)"[153][193]ge [195]osmos [156][206][207][215]";int(tn)/1e+09;"[159]bil/yrs"
  55. 285 printtab(2)"[153][193]ge [195]osmos [156][212][200][197][206]";int(tt)/1e+09;"[159]bil/yrs"
  56. 290 tv=tn-tt
  57. 295 print:print"[153][204]ight travel time";int(tv)/1e+09;"[159]bil/yrs"
  58. 300 sc=(tn/tt)^(2/3)
  59. 305 print"[153][211]cale [156][206][207][215] [153]vs [156][212][200][197][206]";sc
  60. 310 rs=sc-1
  61. 315 print:print"[153][210]edshift seen [156][206][207][215]";rs
  62. 320 dt=3*tt*((tn/tt)^(1/3)-1)
  63. 325 print:printtab(2)"[153][196]istance of object [156][212][200][197][206]";int((dt)/1e+09);"[159]b/lt/yrs"
  64. 330 dn=dt*sc
  65. 335 printtab(2)"[153][196]istance of object [156][206][207][215]";int((dn)/1e+09);"[159]b/lt/yrs"
  66. 340 ss=(2*(tn/tt)^(1/3)-1)
  67. 345 print:printtab(2)"[153][193]way [156][212][200][197][206]";ss;"[159][216] speed/light"
  68. 350 sn=2*((tn/tt)^(1/3)-1)/(tn/tt)^(1/3)
  69. 355 printtab(2)"[153][193]way [156][206][207][215]";sn;"[159][216] speed/light"
  70. 360 ht=(2/3)/tt*9.810001*1e+11
  71. 365 print:print"[153][200]ubble param. [156][212][200][197][206]";int(ht);"[159]km/sec/[205]pc"
  72. 370 hn=(2/3)/tn*9.810001*1e+11
  73. 375 print"[153][200]ubble param. [156][206][207][215]";int(hn);"[159]km/sec/[205]pc"
  74. 380 gosub3000
  75. 385 goto95
  76. 480 q9$="":poke198,.
  77. 485 geta$
  78. 490 poke646,rnd(1)*15+1:print"*[157]";:ifa$=""then485
  79. 495 ifa$=chr$(13)thenprint" ":q9=val(q9$):return
  80. 500 if(a$=chr$(20)andlen(q9$))thenq9$=left$(q9$,len(q9$)-1):goto522
  81. 505 iflen(q9$)>=l9%thensysad+9,4:goto485
  82. 510 if(a$>="0"anda$<="9")ora$="."then515
  83. 515 q9$=q9$+a$
  84. 520 print""a$;:goto485
  85. 522 print" [157][157] [157]";:goto485
  86. 525 poke198,0
  87. 530 geta$:ifa$<>"y"anda$<>"n"then530
  88. 540 return
  89. 3000 poke214,19:print:printtab(8)"[150](1[150]) [195]alculate another
  90. 3010 [153][163]8)"def(2def) (NULL)o (NULL)(NULL)right$(NULL)val(NULL)(NULL)val (NULL)enu
  91. 3020 poke198,0
  92. 3030 geta$:ifa$<"1"ora$>"2"then3030
  93. 3040 ifa$="1"thenreturn
  94. 3050 sysad+15
  95. 3060 print"[147]load"chr$(34)"b.universe"chr$(34)","dv
  96. 3070 print"run28"
  97. 3080 poke631,13:poke632,13:poke198,2:end
  98. 10000 d=peek(186):n$="expanding":open15,d,15,"s0:"+n$:close15:saven$,d:end
  99.