function f=impt(TT)
global coeff T xb xt xz JT
 P = 1;
 yb(JT) = feval('Pb',TT)/P*xb(JT);
 yt(JT) = feval('Pt',TT)/P*xt(JT);
 yz(JT) = feval('Pz',TT)/P*xz(JT);
 sum_y = yb(JT)+yt(JT)+yz(JT);
 f = sum_y - 1.0;
