| Date 08/09/2010 |
|
||||||
|
| Script S6_2_14.m | |||
|
|||
|
%==============================================
%bradley %============================================== % %============================================== %============================================== %Corrigenda %in the sixth row of pag. 267 %change %(one sixtieth of a minute) with %(" means one sixtieth of a minute) %============================================== %============================================== % %light speed in vacuum in m/sec c=2.9979e+08; %mean orbital speed of the Earth in m/sec v=29770; %their ratio vr=v/c % vr2=2*vr %beam emitted along the y axis teta=pi/2 % %tetap at time t num1=cos(teta)-vr; fatt=vr*cos(teta); den1=1-fatt; arg1=num1/den1; tetap=acos(arg1); tetapg=tetap*180/pi % %tetap at time t+6 months num2=cos(teta)+vr; fatt=vr*cos(teta); den2=1+fatt; arg2=num2/den2; tetas=acos(arg2); tetasg=tetas*180/pi % %the first difference %in degrees deltateta1g=(tetapg-90) %in second (a second is one 3600th of a degree) deltateta1m=(tetapg-90)*3600 %in radians deltateta1r=(tetapg-90)*pi/180 % % %the second difference %in degrees deltateta2g=(tetasg-90) %in seconds (a second is one 3600th of a degree) deltateta2m=(tetasg-90)*3600 %in radians deltateta2r=(tetasg-90)*pi/180 %============================================== % |
|||
| Top | |||
