theta=pi*(-m:m)/m; z=r*exp(i*theta)-0.5; rho=abs(z); theta=angle(z);
w=rho.^(1/n).*exp(i*theta./n);
surf(real(z),imag(z),real(w),imag(w)); hold on
w=rho.^(1/n).*exp(i*pi+i*theta./n); surf(real(z),imag(z),real(w),imag(w)); % colorbar('vert') title('(z-0.5)^1/2')
因篇幅问题不能全部显示，请点此查看更多更全内容