MATLAB 平面三点生成圆弧

MATLAB新手，最近要写一个程序，参考了多个百度的结果以及https://blog.csdn.net/libing403/article/details/78234640?utm_source=blogxgwz6这篇文章提供的方法，自己想了一个办法，目前测试过的点都成功了，所以想给大家分享一下。欢迎大家指出错误和不足(๑╹◡╹)ﾉ"""。

%判断是圆还是圆弧
if get(handles.rad_round_1,'value')
    round_flag=1;
elseif get(handles.rad_round_2,'value')
    round_flag=2;
else
    warndlg('请选择圆或圆弧^-^','提示','modal')
end

%输入三个点的坐标

%%获取三点所在圆圆心及半径
x0=[x1 x2 x3]';
y0=[y1 y2 y3]';
TR=triangulation([1,2,3],x0,y0); %表示成三角网格
[D,r]=circumcenter(TR);         %利用内置函数求外接圆心和半径
a=D(1);b=D(2);                  %圆心（a,b）  半径 r

%开始生成
switch(round_flag)
    case 1                      %根据圆心半径生成圆
        deg=0:0.01:360;
        rx=a+r*cos(deg);
        ry=b+r*sin(deg);
        plot(a,b,'*r');hold on;
        plot(x0,y0,'ro');hold on;
        plot(rx,ry,'g');hold on;
    case 2                      %根据圆心半径生成圆弧
        A=[x1 y1 0];B=[x2 y2 0];C=[x3 y3 0];
        AB=B-A;BC=C-B;
        dis=cross(AB,BC);
        if dis(3)==0
            warndlg('三点共线，请重新输入坐标^-^','提示','modal')
        else
            e=c_theta(x1,y1,a,b)；
            f=c_theta(x2,y2,a,b)；
            g=c_theta(x3,y3,a,b)；
            if max(max(e,f),g)==e
                if f>g
                    deg=e:-0.1:g;
                else
                    deg=e:0.1:(360+g);
                end
            elseif max(max(e,f),g)==f
                if e>g
                    deg=e:0.1:(g+360);
                else
                    deg=e:-0.1:(g-360);
                end
            elseif max(max(e,f),g)==g
                if e>f
                    deg=e:-0.1:(g-360);
                else
                    deg=e:0.1:g;
                end
            end
            rx=a+r*cosd(deg);
            ry=b+r*sind(deg);
            plot(a,b,'*r');hold on;
            plot(x0,y0,'ro');hold on;
            plot(rx,ry,'g');hold on;
        end
end

set(gca,'XAxisLocation','origin');
set(gca,'YAxisLocation','origin');
axis equal
hold on;grid on;

%判断输入点的角度
function theta=c_theta(x,y,a,b)
theta=atand((y-b)/(x-a));
if x >= a
    if y >= b
        theta=theta;     %第一象限
    elseif y= b
        theta=theta+180;     %第二象限
    elseif y