function u=DIPPR101(U,T); A=U(1);B=U(2);C=U(3);D=U(4);E=U(5); u=exp(A+B./T+C*log(T)+D*T.^E);