function[z0]=subz0sr(zs0,zr0,zsr0,zl0,zlM0,m,p,q)
z0=zeros(9);
z0(1,1)=zs0;
z0(2,2)=zr0;
z0(3,3)=zl0*m;
z0(4,4)=zl0*p;
z0(5,5)=zl0*q;
z0(6,6)=zl0*(1-m);
z0(7,7)=zl0*(1-p);
z0(8,8)=zl0*(1-q);
z0(9,9)=zsr0;
if m<=p & p<=q
z0(3,4)=zlM0*m;
z0(3,5)=zlM0*m;
z0(4,5)=zlM0*p;
z0(4,6)=-zlM0*(p-m);
z0(5,6)=-zlM0*(q-m);
z0(5,7)=-zlM0*(q-p);
z0(6,7)=zlM0*(1-p);
z0(6,8)=zlM0*(1-q);
z0(7,8)=zlM0*(1-q);
z0(4,3)=z0(3,4);
z0(5,3)=z0(3,5);
z0(5,4)=z0(4,5);
z0(6,4)=z0(4,6);
z0(6,5)=z0(5,6);
z0(7,5)=z0(5,7);
z0(7,6)=z0(6,7);
z0(8,6)=z0(6,8);
z0(8,7)=z0(7,8);
elseif m<=q & q<=p
z0(3,4)=zlM0*m;
z0(3,5)=zlM0*m;
z0(4,5)=zlM0*q;
z0(4,6)=-zlM0*(p-m);
z0(4,8)=-zlM0*(p-q);
z0(5,6)=-zlM0*(q-m);
z0(6,7)=zlM0*(1-p);
z0(6,8)=zlM0*(1-q);
z0(7,8)=zlM0*(1-p);
z0(4,3)=z0(3,4);
z0(5,3)=z0(3,5);
z0(5,4)=z0(4,5);
z0(6,4)=z0(4,6);
z0(8,4)=z0(4,8);
z0(6,5)=z0(5,6);
z0(7,6)=z0(6,7);
z0(8,6)=z0(6,8);
z0(8,7)=z0(7,8);
elseif p<=m & m<=q
z0(3,4)=zlM0*p;
z0(3,5)=zlM0*m;
z0(3,7)=-zlM0*(m-p);
z0(4,5)=zlM0*p;
z0(5,6)=-zlM0*(q-m);
z0(5,7)=-zlM0*(q-p);
z0(6,7)=zlM0*(1-m);
z0(6,8)=zlM0*(1-q);
z0(7,8)=zlM0*(1-q);
z0(4,3)=z0(3,4);
z0(5,3)=z0(3,5);
z0(7,3)=z0(3,7);
z0(5,4)=z0(4,5);
z0(6,5)=z0(5,6);
z0(7,5)=z0(5,7);
z0(7,6)=z0(6,7);
z0(8,6)=z0(6,8);
z0(8,7)=z0(7,8);
elseif p<=q & q<=m
z0(3,4)=zlM0*p;
z0(3,5)=zlM0*q;
z0(3,7)=-zlM0*(m-p);
z0(3,8)=-zlM0*(m-q);
z0(4,5)=zlM0*p;
z0(5,7)=-zlM0*(q-p);
z0(6,7)=zlM0*(1-m);
z0(6,8)=zlM0*(1-m);
z0(7,8)=zlM0*(1-q);
z0(4,3)=z0(3,4);
z0(5,3)=z0(3,5);
z0(7,3)=z0(3,7);
z0(8,3)=z0(3,8);
z0(5,4)=z0(4,5);
z0(7,5)=z0(5,7);
z0(7,6)=z0(6,7);
z0(8,6)=z0(6,8);
z0(8,7)=z0(7,8);
elseif q<=m & m<=p
z0(3,4)=zlM0*m;
z0(3,5)=zlM0*q;
z0(3,8)=-zlM0*(m-q);
z0(4,5)=zlM0*q;
z0(4,6)=-zlM0*(p-m);
z0(4,8)=-zlM0*(p-q);
z0(6,7)=zlM0*(1-p);
z0(6,8)=zlM0*(1-m);
z0(7,8)=zlM0*(1-p);
z0(4,3)=z0(3,4);
z0(5,3)=z0(3,5);
z0(8,3)=z0(3,8);
z0(5,4)=z0(4,5);
z0(6,4)=z0(4,6);
z0(8,4)=z0(4,8);
z0(7,6)=z0(6,7);
z0(8,6)=z0(6,8);
z0(8,7)=z0(7,8);
elseif q<=p & p<=m
z0(3,4)=zlM0*p;
z0(3,5)=zlM0*q;
z0(3,7)=-zlM0*(m-p);
z0(3,8)=-zlM0*(m-q);
z0(4,5)=zlM0*q;
z0(4,8)=-zlM0*(p-q);
z0(6,7)=zlM0*(1-m);
z0(6,8)=zlM0*(1-m);
z0(7,8)=zlM0*(1-p);
z0(4,3)=z0(3,4);
z0(5,3)=z0(3,5);
z0(7,3)=z0(3,7);
z0(8,3)=z0(3,8);
z0(5,4)=z0(4,5);
z0(8,4)=z0(4,8);
z0(7,6)=z0(6,7);
z0(8,6)=z0(6,8);
z0(8,7)=z0(7,8);
end
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
