!MonkeyGod.  A True BasicTM program by Victor J. Stenger
!OPEN #1: name"monkey.out", create newold
!ERASE #1
      CLEAR
      CLOSE #1
      OPEN #1: screen 0,1,0,1
      WINDOW  #1

DIM p(4),ctrue(4),con$(4)
!Conversion factors.  Basic units S.I.
LET ly=9.4e15
LET eV=1.6e-19
LET GeV=1e9*eV
LET MeV=1e6*eV
LET cm=0.01
LET km=1000
LET hr=3600
LET day=24*hr
LET yr=365*day

!Arbitrary Constants
LET c=2.99792458e8
LET h=6.626075e-34
LET hb=h/2/pi
LET G=6.67259e-11                 ! m^2/kg^2    usual S.I. units

!Set parameter true values and names
LET con$(1)="1/alpha "
LET ctrue(1)=137
LET con$(2)="alpha-s  "
LET ctrue(2)=0.2
LET con$(3)="mp (kg)  "
LET ctrue(3)=1.67e-27
LET con$(4)="me (kg)  "
LET ctrue(4)=9.11e-31
DO
   !Generate universe
   !Read in parameters.  If zero, use true values.
   FOR i = 1 to 4
       INPUT prompt con$(i): par
       IF par <> 0 then
          LET p(i)=par
       ELSE
          LET p(i)=ctrue(i)
       END IF
   NEXT i

   LET alpha=1/p(1)               !Fine structure constant
   LET alpha-s=p(2)               !Strong interaction constant
   LET mp=p(3)                    !mass of proton in kg
   LET me=p(4)                    !mass of electron in kg

   PRINT #1: con$(1),con$(2),con$(3),con$(4)
   PRINT #1: 1/alpha,alpha-s,mp,me

   !Compute properties of universe
   LET rb=hb/(alpha*me*c)         !Bohr radius
   LET Eb=alpha^2*me*c^2/2        !Hydrogen binding energy
   LET rN=hb/(alpha-s*mp*c)       !Nucleon radius
   LET EN=alpha-s^2*mp*c^2/2      !Nucleon binding energy
   LET alpha-G=G*mp^2/hb/c        !gravitational coupling
   !Lifetime and mass of main sequence star
   LET tstar=alpha^2/alpha-G*(mp/me)^2*hb/mp/c^2
   LET mstar=alpha-G^(-1.5)*mp
   !Radius and mass of typical planet
   LET Rplanet=sqr(0.1)*2*rB*(me/mp)^.25*(alpha/alpha-
G)^.5
   LET
Mplanet=(0.1)^(1.5)*2*mp*(me/mp)^.75*(alpha/alpha-
G)^1.5
   !Length of day and year
   LET Tday=2*pi*2^1.5*rb/c*(mp/(me*alpha*alpha-
G))^0.5
   LET Tyr=.2*rb/c*(mp/me)^2/alpha^6.5/alpha-G^.125
   !Large numbers
   LET N1=alpha/alpha-G*mp/me
   LET N2=c*tstar/rN

   PRINT #1: "Bohr radius              = ",rb*100,"cm"
   PRINT #1: "Hydrogen binding energy  = ",Eb/eV,"eV"
   PRINT #1: "Nucleon radius           = ",rN*100,"cm"
   PRINT #1: "Nucleon binding energy   = ",EN/MeV,"MeV"
   PRINT #1: "Minimum stellar lifetime = ",tstar/yr,"yr"
   PRINT #1: "Mass of star             = ",mstar,"kg"
   PRINT #1: "Radius of planet         = ",Rplanet/km,"km"
   PRINT #1: "Mass of planet           = ",Mplanet,"kg"
   PRINT #1: "Length of day            = ",Tday/hr,"hr"
   PRINT #1: "Length of year           = ",Tyr/day,"days"
   PRINT #1: "N1                       = ",N1
   PRINT #1: "N2                       = ",N2
LOOP
END