I spent two weeks trying to fix this error the program show my this:: NIntegrate::inumr: The integrand (0.0131326 r Sin[q r])/(1+(E^(1.75439 Plus[<<2>>])+E^(1.75439 Plus[<<2>>])) (0.5 +0.0882389 Power[<<2>>])^1.5) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,50}}.

`Subscript[F, pp1][q_] = ((4*Pi)/q)*NIntegrate[Subscript[\[Rho], p][r]*Sin[q*r]*r, {r, 0, 50}] `

those some definition I used

` Subscript[\[Rho], p][r_] = Subscript[\[Rho]0, p]/(1 + (0.5 + 0.5*(r/R)^2)^1.5*(Exp[(r - R)/Subscript[a, p]] + Exp[-(r + R)/Subscript[a, p]])); 4*Pi*NIntegrate[Subscript[\[Rho], p][r]*r^2, {r, 0, Infinity}]; Nn = 6; R = 1.31*Nn^(1/3); Subscript[a, p] = 0.570; Subscript[\[Alpha], p] = 1.5; `