matlab 画平面圆,求助：Matlab如何绘制空间圆

a=input('圆心坐标：');

r=input('输入半径：');

b=input('输入平面参数：');%按平面的标准方程Ax+By+Cz+D=0输入数组[A,B,C,D]

x=a(1)-r:0.0001:a(1)+r;

y1=1/2/(b(3)^2+b(2)^2).*(2*b(3)^2*a(2)-2*x*b(1)*b(2)-2*b(2)*b(3)*a(3)-2*b(2)*b(4)+2*(-b(1)^2*b(3)^2*x.^2-2*b(1)*b(2)*a(2)*x*b(3)^2-b(3)^4*x.^2-x.^2*b(2)^2*b(3)^2-2*b(3)^3*x*a(3)*b(1)-b(3)^4*a(1)^2+b(3)^4*r^2-b(3)^4*a(3)^2-b(3)^2*b(4)^2-2*a(2)*b(3)^3*b(2)*a(3)-2*a(2)*b(3)^2*b(2)*b(4)-2*b(3)^2*x*b(1)*b(4)+2*b(2)^2*b(3)^2*x*a(1)+2*b(3)^4*x*a(1)-2*b(3)^3*b(4)*a(3)-b(2)^2*b(3)^2*a(1)^2+b(2)^2*b(3)^2*r^2-b(2)^2*b(3)^2*a(2)^2).^(1/2));

%解出的y关于x的表达式(利用solve解出的表达式，未给出过程)

for ii=length(x):-1:1

f(ii)=isreal(y1(ii));

end

x=x(f);

y1=y1(f);

z=-b(3)*(x*b(2)+b(1)*y1-b(1)*b(2))/b(1)/b(2);

plot3(x,y1,z);

hold on;

y2=1/2/(b(3)^2+b(2)^2).*(2*b(3)^2*a(2)-2*x*b(1)*b(2)-2*b(2)*b(3)*a(3)-2*b(2)*b(4)-2*(-b(1)^2*b(3)^2*x.^2-2*b(1)*b(2)*a(2)*x*b(3)^2-b(3)^4*x.^2-x.^2*b(2)^2*b(3)^2-2*b(3)^3*x*a(3)*b(1)-b(3)^4*a(1)^2+b(3)^4*r^2-b(3)^4*a(3)^2-b(3)^2*b(4)^2-2*a(2)*b(3)^3*b(2)*a(3)-2*a(2)*b(3)^2*b(2)*b(4)-2*b(3)^2*x*b(1)*b(4)+2*b(2)^2*b(3)^2*x*a(1)+2*b(3)^4*x*a(1)-2*b(3)^3*b(4)*a(3)-b(2)^2*b(3)^2*a(1)^2+b(2)^2*b(3)^2*r^2-b(2)^2*b(3)^2*a(2)^2).^(1/2));

z=-b(3)*(x*b(2)+b(1)*y2-b(1)*b(2))/b(1)/b(2);

plot3(x,y2,z);

view(3);

grid on;