matlab 中点法画椭圆

以下是使用中点算法在MATLAB中画椭圆的代码：

% 定义椭圆的长轴和短轴
a = 5;
b = 3;

% 初始化参数
x = 0;
y = b;
d1 = b^2 - a^2*b + a^2/4;
dx = 2*b^2*x;
dy = 2*a^2*y;

% 画第一象限的1/8椭圆
while dx < dy
    plot(x, y, '.b', 'markersize', 10);
    hold on;
    plot(x, -y, '.b', 'markersize', 10);
    plot(-x, -y, '.b', 'markersize', 10);
    plot(-x, y, '.b', 'markersize', 10);
    axis equal;
    set(gca, 'XMinorGrid', 'on');
    grid on;
    pause(0.1);
    if d1 < 0
        x = x + 1;
        dx = dx + 2*b^2;
        d1 = d1 + dx + b^2;
    else
        x = x + 1;        y = y - 1;
        dx = dx + 2*b^2;
        dy = dy - 2*a^2;
        d1 = d1 + dx - dy + b^2;
    end
end

% 画第一象限的1/4椭圆
d2 = b^2*(x+1/2)^2 + a^2*(y-1)^2 - a^2*b^2;
while y >= 0
    plot(x, y, '.b', 'markersize', 10);
    hold on;
    plot(x, -y, '.b', 'markersize', 10);
    plot(-x, -y, '.b', 'markersize', 10);
    plot(-x, y, '.b', 'markersize', 10);
    axis equal;
    set(gca, 'XMinorGrid', 'on');
    grid on;
    pause(0.1);
    if d2 > 0
        y = y - 1;
        dy = dy - 2*a^2;
        d2 = d2 + a^2 - dy;
    else
        x = x + 1;
        y = y - 1;
        dx = dx + 2*b^2;
        dy = dy - 2*a^2;
        d2 = d2 + dx - dy + a^2;
    end
end