Equilibrium concentrations of species are computed before signals from equilibrium experiments can be calculated.
The model definitions and derivations were performed in Equilibrium_thermodynamic_models/
The functions to perform this task are placed in code/+equilibrium_thermodynamic_equations/. These functions may be used as stand-alone functions to investigate reaction mixture composition separately from fitting sessions.
Example of functional interface is:
[concentrations, species_names] = U_model(Rtotal, LRratio, K_A,... model_numeric_solver, model_numeric_options)
where
The solution of equilibrium thermodynamic equations is either analytical or numeric. A switch between the two is through model_numeric_solver='analytical' or ='ode45' (as an example). Some models wouldn't have analytical solution.
You can display progression of a titration using plot_populations.m (see Tutorial 6 for examples).
For numeric solutions the function to solve is provided separately (for example, see U_R_model_numeric.m)