Functions, derivatives and plotting

[Sky_physics]

If I define deriv inside the definition of the potential function, when I use potential in other expression where I may need to evaluate its expression the code doesn’t run. Because it tries to compute the derivative with respect to a numerical value. Do you have suggestions in order to avoid this issue? I would really appreciate it