matlab调试错误种类,程序调试出错

function f=fx2(x)

k=10;

for i=1:k

x(i)=sym (['x',num2str(i)]);

end

%  %13即非方阵不能用mulnewton解

f(1)= 8/3 - (16*x(1)*x(6))/(4*x(1)^2 - 4*x(6)^2);

f(2)=26 - (16*x(2)*x(6) + x(1)*(24*x(6) + 16*x(7)))/(4*x(1)^2 + 8*x(2)*x(1) - 8*x(6)^2 - 8*x(7)*x(6));

f(3)=(24*x(2)*x(6) + 16*x(2)*x(7) + 16*x(3)*x(6) + x(1)*(57*x(6) + 24*x(7) + 16*x(8)))/(16*x(6)*x(7) + 8*x(6)*x(8) - 9*x(1)^2 - 4*x(2)^2 + 20*x(6)^2 + 4*x(7)^2 - x(1)*(8*x(2) + 8*x(3))) - 273/20;

f(4)=(57*x(2)*x(6) + 24*x(2)*x(7) + 24*x(3)*x(6) + 16*x(2)*x(8) + 16*x(3)*x(7) + 16*x(4)*x(6) + x(1)*(36*x(6) + 57*x(7) + 24*x(8) + 16*x(9)))/(40*x(6)*x(7) - 8*x(2)*x(3) + 16*x(6)*x(8) + 8*x(6)*x(9) + 8*x(7)*x(8) - x(1)*(18*x(2) + 8*x(3) + 8*x(4)) - 4*x(1)^2 - 4*x(2)^2 + 16*x(6)^2 + 8*x(7)^2) - 593/84;

f(5)=(36*x(2)*x(6) + 57*x(2)*x(7) + 57*x(3)*x(6) + 24*x(2)*x(8) + 24*x(3)*x(7) + 24*x(4)*x(6) + 16*x(2)*x(9) + 16*x(3)*x(8) + 16*x(4)*x(7) + 16*x(5)*x(6) + x(1)*(36*x(6) + 36*x(7) + 57*x(8) + 24*x(9) + 16*x(10)))/(32*x(6)*x(7) - 8*x(2)*x(4) - 8*x(2)*x(3) + 40*x(6)*x(8) + 16*x(6)*x(9) + 16*x(7)*x(8) + 8*x(6)*x(10) + 8*x(7)*x(9) - x(1)*(8*x(2) + 18*x(3) + 8*x(4) + 8*x(5)) - 4*x(1)^2 - 9*x(2)^2 - 4*x(3)^2 + 16*x(6)^2 + 20*x(7)^2 + 4*x(8)^2) - 2211/464;

f(6)=(36*x(2)*x(6) + 36*x(2)*x(7) + 36*x(3)*x(6) + 57*x(2)*x(8) + 57*x(3)*x(7) + 57*x(4)*x(6) + 24*x(2)*x(9) + 24*x(3)*x(8) + 24*x(4)*x(7) + 24*x(5)*x(6) + 16*x(2)*x(10) + 16*x(3)*x(9) + 16*x(4)*x(8) + 16*x(5)*x(7) + x(1)*(36*x(7) + 36*x(8) + 57*x(9) + 24*x(10)))/(32*x(6)*x(7) - 8*x(2)*x(4) - 8*x(2)*x(5) - 8*x(3)*x(4) - 18*x(2)*x(3) + 32*x(6)*x(8) + 40*x(6)*x(9) + 40*x(7)*x(8) + 16*x(6)*x(10) + 16*x(7)*x(9) + 8*x(7)*x(10) + 8*x(8)*x(9) - x(1)*(8*x(2) + 8*x(3) + 18*x(4) + 8*x(5)) - 4*x(2)^2 - 4*x(3)^2 + 16*x(7)^2 + 8*x(8)^2) - 1465/287;

f(7)=(36*x(2)*x(7) + 36*x(3)*x(6) + 36*x(2)*x(8) + 36*x(3)*x(7) + 36*x(4)*x(6) + 57*x(2)*x(9) + 57*x(3)*x(8) + 57*x(4)*x(7) + 57*x(5)*x(6) + 24*x(2)*x(10) + 24*x(3)*x(9) + 24*x(4)*x(8) + 24*x(5)*x(7) + 16*x(3)*x(10) + 16*x(4)*x(9) + 16*x(5)*x(8) + x(1)*(36*x(8) + 36*x(9) + 57*x(10)))/(32*x(6)*x(8) - 18*x(2)*x(4) - 8*x(2)*x(5) - 8*x(3)*x(4) - 8*x(3)*x(5) - 8*x(2)*x(3) + 32*x(6)*x(9) + 32*x(7)*x(8) + 40*x(6)*x(10) + 40*x(7)*x(9) + 16*x(7)*x(10) + 16*x(8)*x(9) + 8*x(8)*x(10) - x(1)*(8*x(3) + 8*x(4) + 18*x(5)) - 4*x(2)^2 - 9*x(3)^2 - 4*x(4)^2 + 16*x(7)^2 + 20*x(8)^2 + 4*x(9)^2) - 3037/636;

f(8)=(36*x(2)*x(8) + 36*x(3)*x(7) + 36*x(4)*x(6) + 36*x(2)*x(9) + 36*x(3)*x(8) + 36*x(4)*x(7) + 36*x(5)*x(6) + 57*x(2)*x(10) + 57*x(3)*x(9) + 57*x(4)*x(8) + 57*x(5)*x(7) + 24*x(3)*x(10) + 24*x(4)*x(9) + 24*x(5)*x(8) + 16*x(4)*x(10) + 16*x(5)*x(9) + x(1)*(36*x(9) + 36*x(10)))/(32*x(6)*x(9) - 8*x(2)*x(4) - 18*x(2)*x(5) - 18*x(3)*x(4) - 8*x(3)*x(5) - 8*x(4)*x(5) - 8*x(2)*x(3) + 32*x(7)*x(8) + 32*x(6)*x(10) + 32*x(7)*x(9) + 40*x(7)*x(10) + 40*x(8)*x(9) + 16*x(8)*x(10) + 8*x(9)*x(10) - 4*x(3)^2 - 4*x(4)^2 + 16*x(8)^2 + 8*x(9)^2 - x(1)*(8*x(4) + 8*x(5))) - 1887/248;

f(9)=- (36*x(1)*x(10) + 36*x(2)*x(9) + 36*x(3)*x(8) + 36*x(4)*x(7) + 36*x(5)*x(6) + 36*x(2)*x(10) + 36*x(3)*x(9) + 36*x(4)*x(8) + 36*x(5)*x(7) + 57*x(3)*x(10) + 57*x(4)*x(9) + 57*x(5)*x(8) + 24*x(4)*x(10) + 24*x(5)*x(9) + 16*x(5)*x(10))/(8*x(1)*x(5) + 8*x(2)*x(4) + 8*x(2)*x(5) + 8*x(3)*x(4) + 18*x(3)*x(5) + 8*x(4)*x(5) - 32*x(6)*x(10) - 32*x(7)*x(9) - 32*x(7)*x(10) - 32*x(8)*x(9) - 40*x(8)*x(10) - 16*x(9)*x(10) + 4*x(3)^2 + 9*x(4)^2 + 4*x(5)^2 - 16*x(8)^2 - 20*x(9)^2 - 4*x(10)^2) - 20665/2618;

f(10)=- (36*x(2)*x(10) + 36*x(3)*x(9) + 36*x(4)*x(8) + 36*x(5)*x(7) + 36*x(3)*x(10) + 36*x(4)*x(9) + 36*x(5)*x(8) + 57*x(4)*x(10) + 57*x(5)*x(9) + 24*x(5)*x(10))/(8*x(2)*x(5) + 8*x(3)*x(4) + 8*x(3)*x(5) + 18*x(4)*x(5) - 32*x(7)*x(10) - 32*x(8)*x(9) - 32*x(8)*x(10) - 40*x(9)*x(10) + 4*x(4)^2 + 4*x(5)^2 - 16*x(9)^2 - 8*x(10)^2) - 3345/79;

%  f(11)=579/13 - (36*x(3)*x(10) + 36*x(4)*x(9) + 36*x(5)*x(8) + 36*x(4)*x(10) + 36*x(5)*x(9) + 57*x(5)*x(10))/(4*x(4)^2 + 8*x(4)*x(5) + 9*x(5)^2 + 8*x(3)*x(5) - 16*x(9)^2 - 32*x(9)*x(10) - 20*x(10)^2 - 32*x(8)*x(10));

%  f(12)=315/31 - (36*x(4)*x(10) + 36*x(5)*x(9) + 36*x(5)*x(10))/(4*x(5)^2 + 8*x(4)*x(5) - 16*x(10)^2 - 32*x(9)*x(10));

%  f(13)=15/7 - (36*x(5)*x(10))/(4*x(5)^2 - 16*x(10)^2);

end