The order of phase transition of a two-dimensional Ising model with m-body interactions in the horizontal direction and n-body interactions in the vertical direction in the m = 3, n = 3 case is of the first-order which is determined by our Monte Carlo algorithm, while in the m= 3,n=2 case the transition is of the second order.