function q = tex2(x,Ts)
T = ppval(Ts,x);
A = pi*(-9.375e-4 * x + 2.5e-5);
q = T.*A;