#!/usr/bin/gnuplot # ohyb_svetla.plt # Vypocet ohybu svetla. # Miroslav Broz (miroslav.broz@email.cz), May 31st 2011 deg = pi/180.; rad = 180./pi; M = 2.e30 # kg r_0 = 7.e8 # m G = 6.67e-11 # SI c = 3.e8 # m/s Delta = 4*G*M/(c**2 * r_0) print "Delta = ", Delta, " = ", Delta*rad*3600., " arcsec" ######################################################################## AU = 149.6e9 # m a = 0.378*AU e = 0.206 yr = 365.25 * 86400. # s delta_phi = 6.*pi*G*M/(c**2 * a * (1.-e**2)) # a**3/P**2 = GM/4pi^2 P = sqrt(a**3/(G*M)*4*pi**2) delta_phi_100yr = 100.*yr/P * delta_phi print "delta_phi = ", delta_phi, " = ", delta_phi*rad*3600., " arcsec" print "P = ", P, " s = ", P/yr, " yr" print "delta_phi_100yr = ", delta_phi_100yr, " = ", delta_phi_100yr*rad*3600., " arcsec"