i have 2 uncorrlated sources in 62 and 67 degree . addative withe gussian noise snr 10 db i want to calcaulate covariance matrix of X. number of element 8 number of snapshot 200 why my code is not working?

N = 8;
fc = 3e8;                              % Operating frequency
fs = 8192; c = 3e8;                               % Sampling frequency
lambda = c/fc;      % Wavelength
ha = phased.ULA('NumElements',N,'ElementSpacing',lambda/2);
ang1 =62; ang2 =67;           % Direction of the signals
angs = [ang1 ang2];
Nsamp = 200;                             % Number of snapshots
noisePwr = 0.01;                         % Noise power
pos = getElementPosition(ha)/lambda;      % Element position in wavelengths                        
x = sensorsig(pos,Nsamp,angs,noisePwr);
y = awgn(x,10) ;
Y=transpose(y) %transpose
r=zeros(8)
for q =1:1:200;
p=Y(:,q)*Y(:,q)';
r=r+p;
end
r=1/200*r;

[Geoff Hayes]

farid - please clarify what you mean by why my code is not working. Are you observing an error, and if so, what is it? Or, is your code just producing an unexpected or incorrect result?