椭圆方程（相对于中心的极坐标形式）：
参考代码：
t=linspace(0,2*pi,100);a=3;b=5;r=a*b./sqrt(a^2*sin(t).^2+b^2*cos(t).^2);polar(t,r)
