function[z0]=subz0_1(zs0,zr0,zl0,zlM0,m,p)
z0=zeros(6);
z0(1,1)=zs0;
z0(2,2)=zr0;
z0(3,3)=zl0*m;
z0(4,4)=zl0*p;
z0(5,5)=zl0*(1-m);
z0(6,6)=zl0*(1-p);
if m<=p
z0(3,4)=zlM0*m;
z0(4,5)=-zlM0*(p-m);
z0(5,6)=zlM0*(1-p);
z0(4,3)=z0(3,4);
z0(5,4)=z0(4,5);
z0(6,5)=z0(5,6);
else
z0(3,4)=zlM0*p;
z0(3,6)=-zlM0*(m-p);
z0(5,6)=zlM0*(1-m);
z0(4,3)=z0(3,4);
z0(6,3)=z0(3,6);
z0(6,5)=z0(5,6);
end
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
