****  SIRIUS.BAS  ****
!!! BORLAND TURBO-BASIC !!!
SZORAS(MeV)= 3.442187229379025:  ATOMSZAM = 99
M0= 1.003245819176071:  Q0= .2219881523248459
M0(MeV)= 934.5234805625103:  D0= 1E-006

10 REM SIRIUS.BAS FLAMBAR 01.11.20.
20 REM SEMLEGES ATOMTOMEGEK ELMELETE
30 REM PLANCK TORVENY: SUPERNOVA EFFEKTUS
40 REM SIRIUS FUGGVENY ILLESZTESE
50 REM AZ IZOBAR SEMLEGES ATOMOK ATLAGOS TOMEGEIRE
60 REM ILLESZTESI PARAMETEREK:
70 REM M0:= OSNUKLEON TOMEG
80 REM Q0:= SMT PARAMETER: NEVLEGES FELSO KORLAT: 2/9
90 REM ******************************************************
100 A$ = "****  SIRIUS.BAS  ****"
110 B$ = "!!! BORLAND TURBO-BASIC !!!"
120 DEFDBL A-Z: DIM A(100): DEFINT I-K, N
130 FOR I = 0 TO 100: A(I) = 0: NEXT: AU = 931.5#
140 K = 0: C = 0: ZM = 1000
150 REM *****************************************************
160 M0 = 1.003245825277016#: Q0 = .2219883019248615#
170 D0 = .000001#
180 REM *****************************************************
190 GOTO 470
200 IF K > 2000 THEN K = 0: J = J + 1: PRINT 0;
210 IF J > 20 THEN PRINT "   *** END": END
220 RANDOMIZE TIMER: K = K + 1: DX = D0 * (RND - .5)
230 M = M0 * (1 + (RND - .5) * DX)
240 Q = Q0 * (1 + (RND - .5) * DX)
250 REM ***************************************************
260 X = 1 + Q: MQ = .5 * M * Q ^ 5
270 REM A = TOMEGSZAM >= 3
280 Z = 0: N = 0: FOR I = 2 TO 100
290 IF A(I) = 0 THEN 360
300 N = N + 1: MA = A(I): A = INT(MA)
310 IF MA - A > .5 THEN A = A + 1
320 F2 = A * (1 - (2 / A) ^ 2)
330 IF F2 <= 0 THEN 220
340 F = SQR(F2): DF = MQ * F ^ 4 / (X ^ F - 1)
350 MC = A * M - DF: Z = Z + (MC - MA) ^ 2
360 NEXT: IF Z >= ZM THEN 200
370 ZM = Z: M0 = M: Q0 = Q
380 REM ***********************************************
390 CLS : PRINT : PRINT A$: PRINT B$: PRINT
400 PRINT "SZORAS(MeV)="; AU * SQR(ZM / N): PRINT
410 PRINT "ATOMSZAM ="; N: PRINT
420 PRINT "M0="; M0: PRINT
430 PRINT "Q0="; Q0: PRINT
440 PRINT "M0(MeV)="; AU * M0: PRINT
450 PRINT "D0="; D0: PRINT : GOTO 200
460 REM MAG ADATOK * * * * * * * * * * * * * * * * * *
470 A(1) = 4.0026032497#: REM He-4
480 A(2) = 7.016004: REM Li-7
490 A(3) = 9.012182: REM Be-9
500 A(4) = 11.0093055#: REM B-11
510 A(5) = 12#: REM C-12
520 A(6) = 14.0030740053#: REM N-14
530 A(7) = 15.9949146221#: REM O-16
540 A(8) = 18.998403207#: REM F-19
550 A(9) = 19.9924401759#: REM Ne-20
560 A(10) = 22.98976967#: REM Na-23
570 A(11) = 23.985041902#: REM Mg-24
580 A(12) = 26.98153844#: REM Al-27
590 A(13) = 27.9769265327#: REM Si-28
600 A(14) = 30.973761512#: REM P-31
610 A(15) = 31.9720706912#: REM S-32
620 A(16) = 34.968852714#: REM Cl-35
630 A(17) = 39.9623831233#: REM Ar-40
640 A(18) = 38.96370693#: REM K-39
650 A(19) = 39.96259123#: REM Ca-40
660 A(20) = 44.955910212#: REM Sc-45
670 A(21) = 47.94794711#: REM Ti-48
680 A(22) = 50.943963714#: REM V-51
690 A(23) = 51.940511915#: REM Cr-52
700 A(24) = 54.938049614#: REM Mn-55
710 A(25) = 55.934942115#: REM Fe-56
720 A(26) = 58.933200215#: REM Co-59
730 A(27) = 57.935347915#: REM Ni-58
740 A(28) = 62.929601115#: REM Cu-63
750 A(29) = 63.929146618#: REM Zn-64
760 A(30) = 68.9255813#: REM Ga-69
770 A(31) = 73.921178216#: REM Ge-74
780 A(32) = 74.92159641799999#: REM As-75
790 A(33) = 79.91652182#: REM Se-80
800 A(34) = 78.91833762#: REM Br-79
810 A(35) = 80.9162913#: REM Br-81
820 A(36) = 83.9115073#: REM Kr-84
830 A(37) = 84.911789325#: REM Rb-85
840 A(38) = 87.905614324#: REM Sr-88
850 A(39) = 88.905847925#: REM Y-89
860 A(40) = 89.904703723#: REM Zr-90
870 A(41) = 92.90637752400001#: REM Nb-93
880 A(42) = 97.90540781999999#: REM Mo-98
890 A(43) = 97.9072164#: REM Tc-98
900 A(44) = 101.904349522#: REM Ru-102
910 A(45) = 102.9055043#: REM Rh-103
920 A(46) = 105.9034835#: REM Pd-106
930 A(47) = 106.9050936#: REM Ag-107
940 A(48) = 108.9047563#: REM Ag-109
950 A(49) = 111.90275723#: REM Cd-112
960 A(50) = 114.9038785#: REM In-115
970 A(51) = 117.9016063#: REM Sn-118
980 A(52) = 120.903818024#: REM Sb-121
990 A(53) = 122.904215722#: REM Sb-123
1000 A(54) = 129.906222821#: REM Te-130
1010 A(55) = 126.9044684#: REM I-127
1020 A(56) = 131.904154512#: REM Xe-132
1030 A(57) = 132.9054473#: REM Cs-133
1040 A(58) = 137.9052413#: REM Ba-138
1050 A(59) = 138.9063483#: REM La-139
1060 A(60) = 139.9054343#: REM Ce-140
1070 A(61) = 141.9077193#: REM Nd-142
1080 A(62) = 144.9127444#: REM Pm-145
1090 A(63) = 146.9151343#: REM Pm-147
1100 A(64) = 151.9197283#: REM Sm-152
1110 A(65) = 150.9198463#: REM Eu-151
1120 A(66) = 152.9212263#: REM Eu-153
1130 A(67) = 155.9221203#: REM Gd-156
1140 A(68) = 158.9253433#: REM Tb-159
1150 A(69) = 161.9267953#: REM Dy-162
1160 A(70) = 164.9303193#: REM Ho-165
1170 A(71) = 167.9323683#: REM Er-168
1180 A(72) = 173.9385813#: REM Yb-174
1190 A(73) = 174.940767928#: REM Lu-175
1200 A(74) = 179.946548827#: REM Hf-180
1210 A(75) = 180.9479963#: REM Ta-181
1220 A(76) = 183.950932629#: REM W-184
1230 A(77) = 186.95575083#: REM Re-187
1240 A(78) = 191.9614794#: REM Os-192
1250 A(79) = 192.9629243#: REM Ir-193
1260 A(80) = 193.9626643#: REM Pt-194
1270 A(81) = 194.9647743#: REM Pt-195
1280 A(82) = 196.9665523#: REM Au-197
1290 A(83) = 199.9683093#: REM Hg-200
1300 A(84) = 201.9706263#: REM Hg-202
1310 A(85) = 204.9744123#: REM Tl-205
1320 A(86) = 205.9744493#: REM Pb-206
1330 A(87) = 207.9766363#: REM Pb-208
1340 A(88) = 208.9803833#: REM Bi-209
1350 A(89) = 208.9824163#: REM Po-209
1360 A(90) = 209.9871319#: REM At-210
1370 A(91) = 210.9905858#: REM Rn-211
1380 A(92) = 223.01973073#: REM Fr-223
1390 A(93) = 223.0184973#: REM Ra-223
1400 A(94) = 227.02774703#: REM Ac-227
1410 A(95) = 232.03805042#: REM Th-232
1420 A(96) = 231.03587893#: REM Pa-231
1430 A(97) = 238.05078262#: REM U-238
1440 A(98) = 237.04816732#: REM Np-237
1450 A(99) = 239.05215652#: REM Pu-239
1460 A(100) = 241.05682292#: REM Am-241
1470 GOTO 200