I dont know how to code the following ODE
, i have:
Some background: j and y are spherical Bessel functions which I know is coded using scipy.spcecial somehow...
I have been given k=0.05 and k=0.5...
Im really in need of a solution to this code as it is for a project that I am working on with a group and my coding skills are poor to say the least. Any help would be great.
also sorry if this is hard to follow I don't know how the input works in this text box...
Output:dU/dR = -k*R^2*V(R)[j(kR)-y(kR)U(R)]^2
When i find U®, i have:
Output:nl(r)=arctan(U(r))
In order to find nl i need to take limit of $r→∞$ so I need to continue integration until a big number for R...Some background: j and y are spherical Bessel functions which I know is coded using scipy.spcecial somehow...
I have been given k=0.05 and k=0.5...
Output:V(R)=A^2/R^4 where A =35.3 and R=7.0
I think U(0) should be zero as r tends to 0?Im really in need of a solution to this code as it is for a project that I am working on with a group and my coding skills are poor to say the least. Any help would be great.
also sorry if this is hard to follow I don't know how the input works in this text box...