Hi there,
I have some basic code for plotting the gravitional time dilation as one moves closer to the event horizon of a black hole of a given mass. See code:
ClearAll["Global`*"]'; tp = 1; G =
6.67408*10^(-11); M = (6.5*10^9)*(1.989*10^30); c =
2.99*10^8; rs = (2*G*
M)/(c^2); Print[]; Print[" Black hole mass:", M, "kg", " \
Schhwarzchild radius:", rs, "m"]; Print[]; Plot[
tp*Sqrt[1 - (rs/r)], {r, 0, 2*10^14}, AxesLabel -> {" ", "Time (s)"},
AxesLabel -> {Style[" (m)", Bold, 26],
Style["Time (s)", Bold, 16]}, LabelStyle -> Directive[Black, 16],
AxesOrigin -> {0, 0}, GridLines -> {{rs}, {}},
GridLinesStyle -> {Directive[{Dashed, Thick}, Red],
Directive[Thick, Red]}, ColorFunction -> "NeonColors"]
I want to create a 3D plot for different masses of black holes ranging from the one given in the code to one a thousand times more massive, with a 3D sheet. I am not great at this, and while my basic code works fine, when I go to do a 3D plot, nothing seems to work.
Attachments:
