  In this work, we study a system of coupled Cahn-Hilliard equations describing
the phase separation of a copolymer and a homopolymer blend. The numerical
methods we propose are based on suitable combinations of existing schemes for
the single Cahn-Hilliard equation. As a verification for our approach, we
present some tests and a detailed description of the numerical solutions'
behaviour obtained by varying the values of the parameters.
