How to modify the following code for nakagami fading channel without considering a MIMO system as shown below?

function H=nak_m(m,Nr,Nt)
m=2;
Nr=5;
Nt=5;
n=zeros(Nr,Nt);
for i=1:2*m
    n=n+randn(Nr,Nt).^2;
end
n=n/(2*m);
phi=2*pi*rand(Nr,Nt);
H=(n^0.5).*cos(phi)+j*(n^0.5).*sin(phi);
v=abs(H)
end