Integral of Abel de-projection equation

Question

I am Trying to solve one integral for my problem which is based on the de-projection of density profile using the Abel inversion equation from 2d to 3d(https://en.wikipedia.org/wiki/Abel_transform) and i am writing my function first then solving the integral from (r to maximum value 4)

   
 # my function is:

def nu(R,r,dn/dR):
    return (-1/np.pi)*((dn/dR)/((R**2 - r**2)**0.5))

# where dndr is the derivative 

nu=lambda r :integrate.quad(nu, r,4,args=(r))[0]
# Here i am stuck i am getting an error whie solving like this .

This is the equation .for upper bound i am using 4 as a maximum value instead of infinity

Any kind of help is warmly welcome.