function u=DIPPR105(U,T); A=U(1);B=U(2);C=U(3);D=U(4); u=A./(B.^(1+(1-T/C).^D));