function res=funmin(x,m,PW,Tr)
a=x(1);
b=x(2);
res=sum((Tr-exp(-a*(m.*PW).^b)).^2);