[求助]一个二元函数的画图问题

[bit2007]

关于书上一个二元函数ep=0.5*(1-a*a)/(1+a*a+2*a*cos(phy))的画图实现问题，a和phy为变量，我用两个FOR循环来实现，应该是一簇类似平滑的正弦曲线，结果发现图与书上的大象径庭（如附件）。请各位帮我纠错。下面是我的代码。
clc;
clear all;
clf;
for i=1:11
a=(i-1)*0.1;
for j=1:181
phy=(j-91)*2;
ep(i,j)=0.5*(1-a*a)/(1+a*a+2*a*cos(phy));
end
plot(-180:2:180,ep(i,:)); xlabel('初始相位差');ylabel('偏移中心点的距离');
hold on;
end

版主及各位大大请伸出援手啊。

[bit2007]

各位帮下忙:lol:

[slgu]

ep=@(a,phy)0.5*(1-a.*a)./(1+a.*a+2*a.*cos(phy));
[a phy]=meshgrid(0:.05:1,0:pi/40:2*pi);
y=ep(a,phy);
surf(a,phy,y)
shading interp
axis tight

[bit2007]

换了一种思路，谢了啊。不过好像还是有一点问题，想观察a从0到1逐次变化时的曲线的变化情况。