function [z0]=subz0(zs0,zr0,zSAm0,zRAm0,zSBp0,zRBp0,zlM0,...
zpTS,zsTS,ztTS,zpTR,zsTR,ztTR,m,p)
z0=zeros(12);
z0(1,1)=zs0;
z0(2,2)=zr0;
z0(3,3)=zSAm0;
z0(4,4)=zRAm0;
z0(5,5)=ztTS;
z0(6,6)=ztTR;
z0(7,7)=zSBp0;
z0(8,8)=zRBp0;
z0(9,9)=zpTS;
z0(10,10)=zsTS;
z0(11,11)=zpTR;
z0(12,12)=zsTR;
if m<=p
z0(3,7)=zlM0*m;
z0(4,8)=zlM0*(1-p);
z0(4,7)=-zlM0*(p-m);
z0(7,3)=z0(3,7);
z0(8,4)=z0(4,8);
z0(7,4)=z0(4,7);
else
z0(3,7)=zlM0*p;
z0(4,8)=zlM0*(1-m);
z0(3,8)=-zlM0*(m-p);
z0(7,3)=z0(3,7);
z0(8,4)=z0(4,8);
z0(8,3)=z0(3,8);
end
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
