等势线的画法问题

xyp88fly

如标题，一个等势线分布的方程F，是极坐标方程，要在直角坐标下，画出等势线F的分布图。
先把我的程序列出来：
theta=linspace(0,pi/2);rho=linspace(0+eps,2)
[theta,rho]=meshgrid(theta,rho);
[x,y]=pol2cart(theta,rho);
for m=1:5
    n=2*m-1;
F1=(rho/2).^(n-1).*sin((n-1)*theta-n*pi/2);
F2=(1/2)^(2*n)*(2./rho).^(n+1).*sin((n+1)*theta-n*pi/2);
F3=1-(1/2)^(2*n);
F4=(F1+F2)/F3;
F=2*F4/(pi*2);
end
[C,h]=contour(x,y,F);
clabel(C,h);

其中，F1,F2,F3,F4只是为了简化方程而定义的，最终组成的方程F表示的是等势线的分布，是极坐标方程，theta和rho是极坐标，x和y是直角坐标（不知道这样变换对不对？），下面是我编写的程序运行后得到的图形


但是，理想中得到的图形应该是类似下面这样的（因为是对称，所以只是给出了1/4的图形）


请问，我的程序该如何修改才能得到类似的图形呢？
请大家帮帮忙，先谢过！