prefactor | Tool for writing Python refactorings | Parser library

I am trying to define a function that requires me to make use of something called an "Associated Laguerre Polynomial". It's listed here under the library. In visual studio code, intellisense predicts "assoc_laguerre()" as a function so it clearly exists!

Yet, when building the code, it highlights the assoc_laguerre() function with the message: "Identifier not found".

Any help would be very appreciated! Thanks!

Code:

Edit: this question is essentially identical to this one and I encourage viewers to go there for more information. To add to that answer: because my function lacks narrow bumps, singularities, etc; has limited integration bounds; and is known to have a small integer solution (I only need 1-2 sigfigs), it is a near-ideal choice for a Monte Carlo algorithm. In the end 1e6 points were more than enough to solve it.

I have a pretty complicated trig function of five variables x1, x2, x3, y1, z1 I'm trying to integrate over (-pi,pi) (once for every variable). I know the solution is an integer, so I'm only asking for one digit of precision during my numerical integration with scipy.integrate.nquad(); but it still ran for 74.5 hours on my laptop with no result before I stopped it. I also had no idea what my progress was by that point.

My actual code (and the function):

I have three 2D-arrays of same shape, lets call them theta, phi and A. Let theta and phi be angles to the normal vector seen from different distances on a surface:

Edit I believe there is a problem with the normalization of the histogram, since one must divide with the radius of each element.

I am trying trying to calculate the fluctuations of particle number and the radial distribution function of a 2d LennardJones(LJ) system using python3. Although I believe the particle fluctuations come out right, the pair correlation g(r) come right for small distances but then blow up ( the calculation uses numpy's histogram method).

The thing is, I can' t find out why such a behavior emerges- perhaps of some misunderstanding of a method? As it is, I am posting the relevant code right below, and if needed, I could also upload other parts of the code or the entire script.

Note first, that since we are working with the Grand-Canonical Ensemble, as the number of particles changes, so is the array that stores the particles- and perhaps that's another point where a mistake in implementation could exist.

Particle removal or insertion

The power spectral density St of a signal u may be computed as the product of the FFT of the signal, u_fft with its complex conjugate u_fft_c. In Python, this would be written as:

I'm running into an error with RSpec and I cannot figure out why it's givin the error. I use FactoryBot to create a sweep with sets a value for fsp_name. This value is somehow not recognized in the validations when trying to validate the first_order.

The error that I'm getting is

I am trying to plot a log-likelihood function in python for the Berthe-Blocke equation and I keep getting a nonsensical plot. Currently I am taking doing

I'm trying to create normal distribution out of evenly distributed values in Fortran with the rejection method. It actually works more or less well, but I don't get exactly the result that I want.

I generate the normal distribution with this segment of code