[g,a]=meshgrid(0:0.02:0.2, 0:0.2:2);z=20*log10(sqrt((g.^2+2*g.*cos(a)+1)/(g.^2-2*g.*cos(a)+1))); surf(a,g,z);
