# Plotting the derivative of the result of NDsolve [duplicate]

I have the following set of equations:

X[r_] := (-4 a*A[r]/(r^2 ϕ[r] Sqrt[A[r] B[r]])
D[r^2/Sqrt[A[r] B[r]] D[ϕ[r], r], r] + 2 a*m^2)/(6 α^2);

Scalar1[r_] := β*D[r^2*Sqrt[A[r]/B[r]]*D[ψ[r], r], r]- 3*k*α^2*X[r]*(ϕ[r])^2*r^2*Sqrt[B[r]/A[r]];

TensorTT1[r_] := (A[r] (-B[r] + B[r]^2 + r Derivative[1][B][r]))/(r^2 B[r]^2) + (3 α^2 X[r] - a*m^2) (ϕ[r])^2/2 -2 a*(D[ϕ[r], r])^2/(
2 B[r]) - β*(D[ψ[r], r])^2 A[r]/(2 B[r]);

TensorRR1[r_] := (A[r] - A[r] B[r] + r Derivative[1][A][r])/(r^2 A[r]) + (3 α^2 X[r] \[Minus] a*m^2) B[r] (ϕ[r])^2/(
2 A[r]) + 2 a*(D[ϕ[r], r])^2/(2 A[r]) - β*(D[ψ[r], r])^2/2;

Tensorθθ1[r_] := (r (-r B[r] Derivative[1][A][r]^2 - 2 A[r]^2 Derivative[1][B][r] +
A[r] (-r Derivative[1][A][r] Derivative[1][B][r] +
2 B[r] (Derivative[1][A][r] + r A''[r]))))/(4 A[r]^2 B[
r]^2) + (3 α^2 X[r] - a*m^2) r^2 (ϕ[r])^2/(
2 A[r]) - 2 a*r^2 (D[ϕ[r], r])^2/(2 A[r] B[r]) + β*
r^2*(D[ψ[r], r])^2/(2 B[r]);

And I am solving it using NDSolve:

sol1 = NDSolve[{Scalar1[r] == 0, TensorTT1[r] == 0, TensorRR1[r] == 0,
Tensorθθ1[r] == 0, A[100] == 1, B[100] == 1, ϕ[100] == 0.01, ψ[100] ==
0.01, ϕ'[100] == 0.1, ψ'[100] == 0.001, A'[100] == 0.001}, {A[r], B[r], ϕ[r], ψ[r]}, {r, 0.1, 100}]

I get a nice clean plot when I try:

Plot[Evaluate[A[r] /. sol1], {r, 0.1, 0.2}, PlotRange -> All]

But I am more interested in looking at the plot of the derivative, but doing

Plot[Evaluate[A'[r] /. sol1], {r, 0.1, 50}, PlotRange -> All]`

gives me an empty plot.