求助一个代数方程组的解法

[zhaoyaobang]

小弟刚刚接触matlab不久，遇到如下问题，苦苦钻研了好几天，也没结果，寻求高人指点：
已知方程组：x^2+y^2+z^2=a
x*y*z=b
x+y*z=10
其中当a=[1,2,3,4,5]时，b对应的值为[6,7,8,9,10]。
求出x,y,z的数值解，并在同一个坐标中画出以a为自变量（横坐标),x,y,z的函数曲线。
望各位朋友详细解答一下，非常感谢！！！

[TTT_IOU]

自己试着做了程序如下：
syms x y z a b
f1=x^2+y^2+z^2-a;
f2=x*y*z-b;
f3=x+y*z-10;
[kx,ky,kz]=solve(f1,f2,f3,x,y,z);
b1=[6 7 8 9 10];
for i=1:length(kx)
ffx(i,: )=subs(kx(i),b,b1);
end
for i=1:length(ky)
ffy(i,: )=subs(ky(i),b,b1);
end
for i=1:length(kz)
ffz(i,: )=subs(kz(i),b,b1);
end
然后用ezplot()对其画图；希望有帮助！！！

[zhaoyaobang]

自己试着做了程序如下：
syms x y z a b
f1=x^2+y^2+z^2-a;
f2=x*y*z-b;
f3=x+y*z-10;
[kx,ky,kz]=solve(f1,f2,f3,x,y,z);
b1=[6 7 8 9 10];
for i=1:length(kx)
ffx(i,: )=subs(kx(i),b,b1);
end
for i=1:length(ky)
ffy(i,: )=subs(ky(i),b,b1);
end
for i=1:length(kz)
ffz(i,: )=subs(kz(i),b,b1);
end
然后用ezplot()对其画图；希望有帮助！！！
非常感谢你帮助！！！