Modified projection method for solving a system of monotone equations with convex constraints