Quit[]

s[f_]:=FullSimplify[f,r>2];

gtt=1-2/r; grr=-1/gtt;

ṫ=(( (Sqrt[r ṙ^2+(r-2)]))/(Sqrt[4/r+r-4]));

pt=s[-gtt ṫ];

ṙ=s[Solve[pt==-1,ṙ]][[1,1,2]]

ŕ=s[ṙ /ṫ]

plot[f_]:=Plot[f,{r,rg,2rg}, Frame->True, ImagePadding->{{28,2},{14,1}},
PlotStyle->Thick, ImageSize->320, AspectRatio->1/3, PlotRange->{{rg,2rg},All}]

rg=9001/4000;

Grid[{{"    v"->plot[ŕ (Sqrt[-(grr/gtt)])], "dτ/dt"->plot[1/ṫ]},
{"dR/dr"->plot[(Sqrt[-grr])], "dт/dt"->plot[(Sqrt[gtt])]},
{"    ṙ"->plot[ṙ], "dr/dt"->plot[ŕ]},
{"    ρ"->plot[0], "    p"->plot[0]}}]