已知空间圆弧曲线的两点以及圆心半径，用matlab画出圆弧

已知空间圆弧曲线的两点以及圆心半径，用matlab画出圆弧。

两点坐标：p1 = (p1(1),p1(2),p1(3));p2 = (p2(1),p2(2),p2(3));圆心坐标：c=(c(1),c(2),c(3))

代码：

function pot2plot3(p1,p2,c,r)
% 画空间圆弧
%   输入弧起点p1,终点p2,圆心c
%% 平面 法向量  
p12 = p2 - p1;
pc1 = p1 - c;
pc2 = p2 - c;
pc1 = pc1/norm(pc1);
pc2 = pc2/norm(pc2);
n = cross(p12,pc2);%法向量
n = n/norm(n);
u = pc1/norm(pc1);
v = cross(u,n);%以u为起始向量，满足右手定则坐标系，则夹角thetam是从u向v的
v = v/norm(v);
% v = -1.*v/norm(v); 
thetam = acos(dot(pc1,pc2)/norm(pc1)*norm(pc2));
theta = (0:thetam/100:thetam)';
c1=c(1)*ones(size(theta,1),1);
c2=c(2)*ones(size(theta,1),1);
c3=c(3)*ones(size(theta,1),1);
x=c1+r*u(1)*cos(theta)+r*v(1)*sin(theta);%圆上各点的x坐标
y=c2+r*u(2)*cos(theta)+r*v(2)*sin(theta);%圆上各点的y坐标
z=c3+r*u(3)*cos(theta)+r*v(3)*sin(theta);%圆上各点的z坐标
plot3(x,y,z,'r');
% hold on
% plot3([p1(1),c(1)],[p1(2),c(2)],[p1(3),c(3)],'g');%半径1
% hold on
% plot3([p2(1),c(1)],[p2(2),c(2)],[p2(3),c(3)],'b');%半径2
end