The system of two coupled maps under in-phase harmonic forcing is investigated numerically. By calculating of the lyapunov exponent spectrum the structure of control parameters space is investigated, the existence regions of the various multistable states are presented, the structure of the basins of attraction is studed.