function[z0]=subz0(zs0,zr0,zt0,zl0u,zlM0u,L1,L2,L3,a,b,c,d) % pu
z0=zeros(11);
z0(1,1)=zs0;
z0(2,2)=zr0;
z0(3,3)=zt0;
z0(4,4)=zl0u*a*L1;
z0(5,5)=zl0u*c*L1;
z0(6,6)=zl0u*(1-b)*L2;
z0(7,7)=zl0u*(1-d)*L3;
z0(8,8)=zl0u*(1-a)*L1;
z0(9,9)=zl0u*b*L2;
z0(10,10)=zl0u*(1-c)*L1;
z0(11,11)=zl0u*d*L3;
if a<=c
z0(4,5)=zlM0u*a*L1;
z0(5,8)=-zlM0u*(c-a)*L1;
z0(8,10)=zlM0u*(1-c)*L1;
z0(5,4)=z0(4,5);
z0(8,5)=z0(5,8);
z0(10,8)=z0(8,10);
else
z0(4,5)=zlM0u*c*L1;
z0(4,10)=-zlM0u*(a-c)*L1;
z0(8,10)=zlM0u*(1-a)*L1;
z0(5,4)=z0(4,5);
z0(10,4)=z0(4,10);
z0(10,8)=z0(8,10);
end
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
