[求助] 为何运行出的反正切函数不能连续 怎么才能让它连续的 我的程序有什么问题？

figure(3) h=ezplot('y-0.02/(sqrt((1-x^2+0.15*y^2)^2+(0.06*x)^2))',[0 2],[0 pi]); C=get(h,'contourMatrix'); a=C(1,; b=C(2,; if a<1.01     y=atan(0.06.*a./(1-a.^2+0.15.*b.^2)) ; else     y=atan(0.06.*a./(1-a.^2+0.15.*b.^2))+pi; end plot(a,y) hold on

回复此楼