Generating Uncorrelated fading channel coefficients using complex gaussian distribution ~CN(0,1)

I donot have communication toolbox for this. I want to generate channel state information or channel coefficients randomly from complex gaussian distribution. The values of my channel coefficients must be uncorrelated to each other. How shall I implement this?(like the rayleigh fading model)

I have implemented a function which gives me complex numbers as channel coefficient matrix. But I want my channels to be uncorrelated to each other. How can i fix this?

function result =channel(R,T)
         rp=1/sqrt(2)*randn(R,T); 
         imgp=1/sqrt(2)*randn(R,T);
         result=rp+1j*imgp;
end
%function to generate complex gaussian distributed channel coefficients
%~CN(0,1)
%the channel coefficients must be uncorrelated 
%it is a fading channel model