Phase transitions and the thermodynamic properties of the two-dimensional ferromagnetic Izing model with spin state number q = 4 on a kagome lattice are studied by the Monte Carlo method on the base of the Wang–Landau algorithm. The phase transition character is analyzed using the fourth-order Binder cumulant method and the histogram analysis of the data. It is found that this model undergoes a first-order phase transition.