Как определить специальную матричную функцию в цикле for в MATLAB?

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))

Мне было интересно, не могли бы вы помочь мне с моей проблемой. Заранее спасибо.