The four-dimensional ferromagnetic Ising model in external magnetic field is simulated on the Creutz cellular automaton algorithm using finite-size lattices with linear dimension 4 ≤ L ≤ 8. The study aims to investigate the finite-size scaling behavior of the model, which is crucial for understanding the critical properties of the system. The results of the simulation will provide valuable insights into the behavior of the model at different sizes, allowing for a better understanding of the phase transitions and critical phenomena in the system.