Matlab - Problema con solve dentro de un bucle

clc

clear all

syms Ax Ay N1 N2 x Fm

%Data

x0=5;

n=0.47;

p=2;

e=2;

vm=(n*p*e)/(2*3.14);

angle=0.26;

M=11.55;

m=0.3;

am=0

tr=11:22

H=316;

a=-51;

altura=377.5;

base=162;

g=9.81;

%Geometria

U=1;

for i=1:length(tr)

    t=tr(1,U);

x2=(x0+vm*cos(angle)*t+(am*cos(angle)*t^2))/2;

x2r(U,:)=x2;

U=U+1;

end

U=1;

for i=1:length(tr)

    t=tr(1,U);

h=(vm*t*sin(angle)+(am*cos(angle)*t^2))/2;

hr(U,:)=h;

U=U+1;

end

U=1;

for i=1:length(hr)

    h=hr(U,1);

    x2=x2r(U,1);

tita=asin(h/x2);

titar(U,1)=tita;

U=U+1;

end

%Velocidades de giro

U=1;

for i=1:length(titar)

    tita=titar(U,1);

    x2=x2r(U,1);

w2=(vm*sin(angle)/cos(angle))/(x2*cos(tita));

w2r(U,:)=w2;

U=U+1;

end

%Aceleracion de giro

U=1;

for i=1:length(titar)

    tita=titar(U,1);

    w2=w2r(U,1);

    x2=x2r(U,1);

alfa2=(am*sin(angle)*cos(angle)+x2*sin(angle)*w2.^2-2*w2*vm)/(x2*cos(tita));

alfa2r(U,:)=alfa2;

U=U+1;

end

%ag1

ag1x=am

ag1y=am*sin(angle)/cos(angle)

%ag2

U=1;

for  i=1:length(alfa2r)

    alfa2=alfa2r(U,1);

    w2=w2r(U,1);

ag2x=-alfa2*(H/2)-w2^2*a;

ag2y=alfa2*a-w2^2*(H/2);

ag2xr(U,1)=ag2x;

ag2yr(U,1)=ag2y;

U=U+1;

end

H=H/1000;

a=-a/1000;

altura=altura/1000;

base=base/1000;

s=11/1000;

Ig=(altura^2+base^2)*M/12;

nu1=0.6;

nu2=0.8;

L=10.8/1000;

U=1;

for i=1:11

    ag2x=ag2xr(U,1);

    ag2y=ag2yr(U,1);

    alfa2=alfa2r(U,1);

    x2=x2r(U,1)/1000;

Eq1=Ax+N2*nu1-M*ag2x;

Eq2= Ay+N2-M*g-M*ag2y;

Eq3=N2*x2-M*g*a+M*ag2x*H/2-M*ag2y*a-Ig*alfa2;

Eq4=N2*nu1-Fm+N1*nu2*cos(angle)+N1*sin(angle)-m*ag1x;

Eq5=-N2-m*g-N1*nu2*sin(angle)+N1*cos(angle)-m*ag1y;

Eq6=-Fm*L+N1*nu2*cos(angle)*s+N1*nu2*sin(angle)*x+N1*sin(angle)*s+N1*cos(angle)*x-N2*x2-N2*nu1*L;

[Ax,Ay,N1,N2,x,Fm]=solve(Eq1,Eq2,Eq3,Eq4,Eq5,Eq6,Ax,Ay,N1,N2,x,Fm)

Axr(U,:)=Ax;

Ayr(U,:)=Ay;

N1(U,:)=N1;

N2(U,:)=N2;

x(U,:)=x;

Fm(U,:)=Fm;

U=U+1;

end

[Axx,Ayy,N11,N22,xx,Fmm]=solve(Eq1,Eq2,Eq3,Eq4,Eq5,Eq6,Ax,Ay,N1,N2,x,Fm);

Axrr(U,1)=Axx;

Ayrr(U,1)=Ayy;

N1r(U,1)=N11;

N2r(U,1)=N22;

xr(U,1)=xx;

Fmr(U,1)=Fmm;

>>sol_sistemas_no_lineales1