In this work, we consider a system of coupled Cahn-Hilliard equations describing the phase separation of a copolymer and a homopolymer blend. We propose some numerical methods to approximate the solution of the system which are based on suitable combinations of existing schemes for the single Cahn-Hilliard equation. As a verification for our experimental approach, we present some tests and a detailed description of the numerical solutions’ behaviour obtained by varying the values of the system’s characteristic parameters.