%[z2]=subz2(zs2,zr2,zl2,m); function [z2]=subz2(zs2,zr2,zl2,m); z2=zeros(4); z2(1,1)=zs2; z2(2,2)=zr2; z2(3,3)=m*zl2; z2(4,4)=(1-m)*zl2;