% [It0]=subIt0(Ism0,Irm0,zl0,zlM0,m); function [It0]=subIt0(Ism0,Irm0,zl0,zlM0,m) It0=(zlM0/zl0)*(-m*Ism0+(1-m)*Irm0);