function u=dDIPPR101(U,T);

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