m1=1.6726*10^-27
m2=9.109*10^-31
s=RandomReal[5.29*10^-11,{200}]
c=RandomReal[300000000,{200}]
n=RandomInteger[100000000000,{200}]
n1=RandomInteger[100000000000,{200}]
n2=RandomInteger[100000000000,{200}]
t=RandomReal[1000000000000,{200}]
G=6.67*10^-11
e=m1*c^2
e1=m2*c^2
s1=s^-1
v2=s1/s
v1=(e-e1)*v2/(m2-m1)*c^2
f=(e1/e)/Sqrt[(e1/e)-((v2/c)*e1/e)^2]*m1-m2
g=(v2*(n-1))/(c*Sqrt[1-v1^2/c^2])
t1=t*c*v2*(n-1)/(1-v1^2)
gama=1/Sqrt[n1^2/Pi^2+(n1^2*v1^2)/(n2^2*c^2)]
gama1=1/Sqrt[n1^2/Pi^2-(n1^2*v1^2)/(n2^2*c^2)]
gama2=gama^-1
gama3= gama+gama1
t2=f/t1
t3=t2*g
ListPlot[v1]
ListPlot[v2]
ListPlot[t1]
ListPolarPlot[f]
ListPolarPlot[g]
ListPolarPlot[t3]
ListLinePlot[gama3]
data=Table[SinIntegral[f*g],{g,0,-Pi,Pi},{f,0,-Pi,Pi}]
ListPlot[data]
ContourPlot3D[gama*gama1==Cos [f] ,{gama,-Pi,Pi},{f,-Pi,Pi},{gama1,-Pi,Pi}]