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

[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

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