I have this problem and I want to solve this and the conditions are:
when p[r] = 0 we have the radius of the star;
and after this lim r-> \infinity R[r] = 0.
How can I solve this problem?
I mean how can I write a looping for solve the equations inside and outside the star simultaneously.
tryStarobinsky.nb (22.2 KB)