function T0 = Mu0(a,b) frequency; mu3=0; ll0=0; for i=1:length(w) mu3=mu3+Mu01(ll(i),a,b)*w(i); ll0=ll0+wl(i)*ll(i); end T0=mu3; end