#!/usr/bin/gnuplot # eddington.plt # Eddintonova limita v porovnani s L(M) ~ M^3. # Miroslav Broz (miroslav.broz@email.cz G = 6.67e-11 # N kg^-2 m^2 c = 3.e8 # m M_S = 2.e30 # kg L_S = 3.8e26 # W X = 0.70 kappa = 0.02*(1+X) # m^2/kg L_Ed(M) = 4*pi*c*G*M/kappa L(M) = L_S*(M/M_S)**3 M_Ed = sqrt(4*pi*c*G/kappa * M_S**3/L_S) print "M_Ed = ", M_Ed/M_S, " M_S" set xl "M / M_S" set yl "L / L_S" set xr [0.01:1000] set mxtics 10 set mytics 10 set logscale xy p L_Ed(x*M_S)/L_S,\ L(x*M_S)/L_S pa -1 set term png small set out "eddington.png" rep