I have this differential equation:
eq=D[\[Psi][r],r]+A*\[Psi][r]-(B/(A*r))\[Psi][r]==0;
And DSolve sucessfully solves it:
sol=DSolve[eq,\[Psi][r],r];
But when I check by replacing all, in this case it doesn't replace the derivative of the function:
Simplify[eq/.sol]
Derivative still must be replaced too. If I do:
Simplify[eq/.sol/.D[sol,r]]
Then this gives True.
In other cases Mathematica replaced the function and its derivatives in one shot. What's happening in this case?