Matlab code for ternary model for unbleached molecules seems to work properly now. Problem before: volume fraction not defined properly and thus low mobility inside.

Deleting simple Cahn Hilliard example. First trial Flory Huggins with many terms (simple substitution of y=laplace(phi_tot) instead of the entire chemical potential) still doesn't work properply. For now continuing with amended Cahn Hilliard which accounts for Flory Huggins free energy. The other Flory Huggins code probably has a coding/fenics issue I don't see.

Adding new file with Cahn Hilliard extended to Flory Huggins binary mixture. This was used in email to Chris.

Amended Cahn Hilliard code in a minimal way to account for Flory Huggins binary mixture with the corresponding mobility scaling. Seems to work, but not sure about the volume fraction to concentration conversion.

2D/1D Flory Huggins works in principle. Currently output only via paraview, need to figure out plotting from fenics.