Matlab code to generate Nakagami fading channel

I am doing some wireless communications simulations using Matlab. There are two famous channel models 1) Rayleigh fading and 2) Nakagami fading. I can easily generate Rayleigh fading which is a complex Gaussian random variable:

rayRv=(randn(1,1)+1i*randn(1,1))/sqrt(2)

Similarly, I want to generate Nakagami random variable but It makes me confuse as it is not a complex Gaussian but may be a sum of such variables (but I may have to carefully deal with the phase???). Many technical papers use Nakagami fading (say nakRv) but they mostly focus on abs(nakRv) which is easy to generate as we have clear distribution, or power version which is abs(nakRv)^2, this can be generated with gamma distribution.

If anyone has experience of generating nakRv similar to rayRv, please share the code with me. Thanks !!!