% PG62n function [z1]=subz1T(zs1,zr1,zsr1,zl1,m); z1=zeros(5); z1(1,1)=zs1; z1(2,2)=zr1; z1(3,3)=m*zl1; z1(4,4)=(1-m)*zl1; z1(5,5)=zsr1;