function D36 = D36eq(T)
Tinv = 1./T;
D36 = -428367800000000.*Tinv.^6+9035851000000.*Tinv.^5-71879880000.*Tinv.^4+246554900.*Tinv.^3-184362.6.*Tinv.^2+45.40216.*Tinv+0.0000003836289;