Based on the parameterisation in Cheng (2008) Ind. Eng. Chem. Res. 47 3285-3288, with some adjustments (see below)
Modified on 28 June 2016. Density calculation has been changed: equation 25 in Cheng's paper to compute the density of the mixture should use the glycerine fraction by VOLUME and not by mass. Thanks to Paul Debue for pointing this out.
Modified on 12 July 2017. Andreas Volk pointed out that the density calculation can be made more accurate by (i) accounting for the volume contraction of the mixture; (ii) adjusting the fit for the density of pure glycerine as a function of temperature. In both cases the fit was chosen to match data from Gregory (table 3 and 7)
A matlab function for this is also available