#!/usr/bin/gnuplot G = 6.67e-11 M_Z = 6.e24 R_Z = 6378.e3 omega = 2.*pi/(86400.) M_M = M_Z/81. r = 384.e6 g = G*M_Z/R_Z**2 ao = omega**2 * R_Z as_M = 2*G*M_M/r**3 * R_Z M_S = 2.e30 r_S = 150.e9 as_S = 2*G*M_S/r_S**3 * R_Z print "g = ", g , " m/s^2" print "a_odstredive = ", ao, " m/s^2 = ", ao/g, " g" print "a_slapove_M = ", as_M, " m/s^2 = ", as_M/g, " g" print "a_slapove_S = ", as_S, " m/s^2 = ", as_S/g, " g" P_rot_Z = (23.+56./60.+4./3600.) * 3600. # s omega_rot_Z = 2.*pi/P_rot_Z r_geo = (G*M_Z/omega_rot_Z**2)**(1./3.) print "omega_rot_Z = ", omega_rot_Z, " rad/s" print "r_geostacionarni = ", r_geo, " m" omega_rot_Mars = 2.*pi/(24.62*3600.) M_Mars = 0.107*M_Z r_geo_Mars = (G*M_Mars/omega_rot_Mars**2)**(1./3.) print "omega_rot_Mars = ", omega_rot_Mars, " rad/s" print "r_marsostacionarni = ", r_geo_Mars, " m" T = r*M_M/(M_Z+M_M) print "M_M = ", M_M, " kg" print "Teziste = ", T, " m" I_V = 8.e37 # kg*m^2; moment setrvacnosti omega_rot_V = 2.*pi/(243.*86400.) L_V = I_V * omega_rot_V v = 100.e3 # m/s R_V = R_Z m = I_V*omega_rot_V / (R_V*v) rho = 3000. R = (m/(4./3.*pi*rho))**(1./3.) print "omega_rot_V = ", omega_rot_V, " rad/s" print "L_V = ", L_V, " kg*m^2/s" print "m = ", m, " kg" print "R = ", R, " m = ", R/1.e3, " km" M_V = M_Z r_geo_V = (G*M_V/omega_rot_V**2)**(1./3.) print "r_venusostacionarni = ", r_geo_V, " m = ", r_geo_V/1.e9, " mil. km" r_L = r_S * (M_Z/(3.*M_S))**(1./3.) print "r_L = ", r_L, " m = ", r_L/1.e9, " mil. km" r_V = 0.7*r_S r_L_V = r_V * (M_V/(3.*M_S))**(1./3.) print "r_L_Venuse = ", r_L_V, " m = ", r_L_V/1.e9, " mil. km"