#!/usr/bin/gnuplot # pulzace.plt # Vypocet periody pulzaci z 1-zonoveho modelu (Carrol & Ostlie, 2007). # Miroslav Broz (miroslav.broz@email.cz), Apr 25th 2009. G = 6.67e-11 # N kg^-2 m^2 day = 86400. # s M_S = 2.e30 # kg R_S = 7.e8 # m M = 1*M_S R0 = 1*R_S #gamma = 4./3. # toto by bylo singularni gamma = 2.4 omega = sqrt( (3*gamma-4) * G*M/R0**3 ) P = 2*pi/omega print "M = ", M/M_S, " M_S" print "R0 = ", R0/R_S, " R_S" print "P = ", P, " s = ", P/day, " d" ######################################################################## M = 5*M_S R0 = 50*R_S omega = sqrt( (3*gamma-4) * G*M/R0**3 ) P = 2*pi/omega print "M = ", M/M_S, " M_S" print "R0 = ", R0/R_S, " R_S" print "P = ", P, " s = ", P/day, " d"