for j=1:length(wt1) lamdaH=1./b*10^(a0+(j-1)*s0); H1=f(lamdaH)*wt1(j); H2=H2+H1; end HH2(i)=H2/b;