An algorithm based an variable Metric method has been applied to the determination of molecular thermodynamic model parameters from experimental data when any model need several constraints. Implementation of the method is illustrated in the reduction of multi-properxy equation of state data for the systems of argon, nitrogen and methane.