%function [TT,T]=temp(T,Kp)
function f=temp(T,Kp)

%TT=2.113+2074.8/T-2.4943*log10(T)-1.256e-4*T+1.8564e-7*T^2-log10(Kp);
f=2.113+2074.8/T-2.4943*log10(T)-1.256e-4*T+1.8564e-7*T^2-log10(Kp);
%f=2.113+2074.8/T-2.4943*log(T)-1.256e-4*T-log(Kp);
