for n_pts=1:N_pts
v = exp(-j*k(n_pts)*[-(M-1)/2:(M-1)/2]');
- y(n_pts) = abs(v'*R_num*v)/(abs(v'*R_den*v).^2)*exp(1i*k(n_pts)*l);
+ %y(n_pts) = abs(v'*R_num*v)/(abs(v'*R_den*v).^2)*exp(1i*k(n_pts)*l);
+ y(n_pts) = (v'*R_num*v)/(abs(v'*R_den*v).^2)*exp(1i*k(n_pts)*l);
end
\ No newline at end of file