*** VAS.BAS ***
HIBA (MeV)= 8.604261836353544E-011 
M0= 1.003258334704625:  Q0= .2220284784483794
M0(MeV)= 934.5291192273497: D= 1E-011

10 REM VAS.BAS FLAMBAR 01.11.20.
20 REM SEMLEGES ATOMTOMEGEK ELMELETE
30 REM PLANCK TORVENY: SUPERNOVA EFFEKTUS
40 REM SIRIUS ILLESZTESE A VASRA
50 REM ILLESZTESI PARAMETEREK:
60 REM M0:= OSNUKLEON TOMEG
70 REM Q0:= SMT PARAMETER
80 CLS : PRINT : PRINT "*** VAS.BAS  ***": PRINT
90 PRINT "!!! BORLAND TURBO-BASIC !!!"
100 DEFDBL A-Z: DEFINT I-K
110 K = 0: KM = 2000: J = 0: ZM = 1000
120 REM INPUTS: ***************************************
130 AU = 931.494#       REM 1 U = AU * MeV
140 MA = 55.93915575#:  REM VAS-56 ATOM TOMEGE
150 A = 56:             REM TOMEGSZAM: A > 2
160 REM ***********************************************
170 M0= 1.003258334704757
180 Q0= .2220284784505691
190 D = 1E-11
200 REM ***********************************************
210 K = K + 1: IF K > KM THEN 320
220 RANDOMIZE TIMER: YY = RND * D
230 M = M0 * (1 + (RND - .5) * YY)
240 Q = Q0 * (1 + (RND - .5) * YY)
250 REM ***********************************************
260 X = 1 + Q:  MQ = .5 * M * Q ^ 5
270 F2 = A * (1 - (2 / A) ^ 2): F = SQR(F2)
280 DF = MQ * F ^ 4 / (X ^ F - 1)
290 MC = A * M - DF: Z = ABS(MA - MC)
300 IF Z >= ZM THEN 210
310 ZM = Z: M0 = M: Q0 = Q: J = J + 1
320 K = 0: CLS
330 PRINT : PRINT "*** VAS.BAS ***": PRINT
340 PRINT "HIBA (MeV)="; AU * ZM: PRINT
350 PRINT "JAVITAS SZAMA ="; J: PRINT
360 PRINT "M0="; M0: PRINT
370 PRINT "Q0="; Q0: PRINT
380 PRINT "M0(MeV)="; AU * M0: PRINT
390 PRINT "D="; D: PRINT : GOTO 210