如何在MATLAB中的for循环中定义特殊的矩阵函数？

Question

这是我的代码：

syms  y1 y2 y3 y4 y5;
f1(y1,y2,y3,y4,y5)=(-(y5^3)/(y3^2))*(3*(y2-y1+y3^(-1)-y5/10)^2+(1/5)*(y2-y1+y3^(-1)-y5/10))-y4;
f2(y1,y2,y3,y4,y5)=(1/10)*y5-y4;
f3(y1,y2,y3,y4,y5)=(y5^3)*(3*(y2-y1+y3^(-1)-y5/10)^2+(1/5)*(y2-y1+y3^(-1)-y5/10));
f4(y1,y2,y3,y4,y5)=y1-y3^(-1);
J(y1,y2,y3,y4,y5)=jacobian([f1,f2,f3,f4],[y1,y2,y3,y4]);
y0=[1,1,1,1,1];
for i=1:2
    J=J(y0(1),y0(2),y0(3),y0(4),y0(5))
    y0=y0+[1,1,1,1,1]
end

对于第2步，我看到以下错误：

Error using sub2ind (line 72)
Out of range subscript.

Error in sym/subsref (line 766)
                R_tilde = sub2ind(size(L), Idx.subs{:});

Error in Untit11led (line 9)
    J=J(y0(1),y0(2),y0(3),y0(4),y0(5))

我想知道您是否可以帮助我解决我的问题。预先感谢。