MATLAB求解四面体球心

%x0代表球心，其余代表四个顶点
syms x0 y0 z0
syms x1 y1 z1
syms x2 y2 z2
syms x3 y3 z3
syms x4 y4 z4

%solve((x0-x1)^2 + (y0-y1)^2 + (z0-z1)^2 -((x0-x2)^2 + (y0-y2)^2+(z0-z2)^2),x0)%求x0
%得出x0 = ((y0 - y1)^2 - (y0 - y2)^2 + (z0 - z1)^2 - (z0 - z2)^2 + x1^2 - x2^2)/(2*x1 - 2*x2)

%solve((((y0 - y1)^2 - (y0 - y2)^2 + (z0 - z1)^2 - (z0 - z2)^2 + x1^2 - x2^2)/(2*x1 - 2*x2)-x1)^2 + (y0-y1)^2 + (z0-z1)^2 -((((y0 - y1)^2 - (y0 - y2)^2 + (z0 - z1)^2 - (z0 - z2)^2 + x1^2 - x2^2)/(2*x1 - 2*x2)-x3)^2 +(y0-y3)^2 + (z0-z3)^2),y0)%带入x0求y0
%y0 = ((- x1^2*x2 + x1^2*x3 + x1*x2^2 - x1*x3^2 + x1*y2^2 - x1*y3^2 + x1*z2^2 - 2*z0*x1*z2 - x1*z3^2 + 2*z0*x1*z3 - x2^2*x3 + x2*x3^2 - x2*y1^2 + x2*y3^2 - x2*z1^2 + 2*z0*x2*z1 + x2*z3^2 - 2*z0*x2*z3 + x3*y1^2 - x3*y2^2 + x3*z1^2 - 2*z0*x3*z1 - x3*z2^2 + 2*z0*x3*z2)/(2*x1*y2 - 2*x2*y1 - 2*x1*y3 + 2*x3*y1 + 2*x2*y3 - 2*x3*