function w=paixu(t)
b1=t(1);b=t(2);b2=t(3);
n=10;alfa=0.5;aa=0;
for j=1:n+1
    lamta=(j-1)/n;
    ad(j)=b1+(b-b1)*lamta;
    au(j)=b2-(b2-b)*lamta;
    I(j)=lamta;
    M(j)=(1-alfa)*ad(j)+alfa*au(j);
    d(j)=au(j)-ad(j);
end
Ju=sum(M.*I)/n;
Jd=sum(I)/n;
J=Ju/Jd;
for j=n/2+1:n+1
    aa=aa+I(j)*d(j);
end   
H=aa/30;
beta=0.00001;
w=J/(1+beta*H)^0.5;