I have a project to fit data to an activity coefficient model to determine if the experimental data given is reliable and to determine if a solution diverges from ideal behavior.
The activity coefficient model is Margules and an objective function to be used is Barker's method which uses only experimental temperature, pressure and liquid phase mole fractions.
Partial data was given for the experimental pressure and liquid and vapor mole fractions for a binary system of acetone/cyclohexane at 50°C.
I am supposed to fit the data using the activity coefficient but need to solve for the A12 and A21 values first to determine the activity coefficient for each different mole fraction.
I am stuck as to how I am to solve for these two parameters using Matlab??
Any help would be appreciated, I could be going down the wrong path and would like to know if I am. Thanks.