In[]:=
ClearAll["Global`*"]
In[]:=
c=Quantity[3*10^8,"Meters"/"Seconds"];u=Quantity[2.50*10^5,"Meters"/"Seconds"];γ=SetPrecision[1/Sqrt[1-(u/c)^2],9]
Out[]=
1.00000035
In[]:=
mSun=1;
In[]:=
newtMomentum=mSun*u;
In[]:=
relMomentum=γ*mSun*u;
In[]:=
Round[Simplify[(relMomentum-newtMomentum)/relMomentum],.1*10^-7]
Out[]=
3.5×
-7
10
In[]:=
γMinus1=γ-1
Out[]=
3.5×
-7
10
In[]:=
ans1=Simplify
(3.5*10^-7)*c^2-.5u^2
(3.5*10^-7)*c^2
Out[]=
0.00793651
In[]:=
ans2=Simplify
(γMinus1)*c^2-.5u^2
(γMinus1)*c^2
Out[]=
5.2079×
-7
10