The method of the fourth-order cumulant of Challa, Landau, and Binder is used together with the Monte Carlo histogram technique of Ferrenberg and Swendsen to study the order of the phase transitions of two-dimensional Ising systems with multispin interactions in the horizontal direction and two-body interactions in the vertical direction.