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