25/08/09 10:50:07.17 hJupetNI.net
relativeVelocity[u_, v_, theta_: Pi, c_: 1] := Module[{
numerator,
denominator
},
numerator = Sqrt[u^2 + v^2 - 2*u*v*Cos[theta] - (u^2*v^2*Sin[theta]^2)/c^2];
denominator = 1 - (u*v*Cos[theta])/c^2;
numerator/denominator
]
u = Range[0, 1, 1/99];
v = Range[0, 1, 1/99];
z = Outer[relativeVelocity, u, v];
ContourPlot[relativeVelocity[x, y], {x, 0, 1}, {y, 0, 1},
PlotLegends -> Automatic,
Contours -> 50,
FrameLabel -> {"u", "v"}
]
V = Range[0, 1, 1/1000];
Plot[relativeVelocity[v, v], {v, 0, 1},
AspectRatio -> 1
]