The 7-dimensional Ising model with nearest-neighbor pair interactions is simulated on the Creutz cellular automaton by using the finite-size lattices with the linear dimensions L = 4, 6, 8. The temperature variations and the finite-size scaling plots of the specific heat and Binder parameter verify the theoretically predicted expression near the infinite lattice critical temperature. The approximate values for the critical temperature of the infinite lattice T-c = 12.866(2), T-c = 12.871(2) and T-c = 12.871(49) are obtained from the intersection points of specific heat curves, Binder parameter curves and the straight line fit of specific heat maxima, respectively. These results are in agreement with the more precise values of the Creutz cellular automaton results of T-c = 12.8700(42), the 1/d-expansion result of T-c = 12.8712, the series expansion result of T-c = 12.86902(33), the dynamic Monte Carlo result of T-c = 12.8667(50). The values obtained for the critical exponent of the specific heat, i.e., alpha = 0.011(76), alpha = -0.002, alpha = 0.011(5) and alpha = 0.082(32) corresponding to the above T-c values, respectively, are in agreement with alpha = 0 predicted by the theory. Moreover the values for the Binder parameter are calculated as gL(T-c) = - 1.022(28), gL(T-c) = -1.09, gL(T-c) = -1.01(12) and gL(T-c) = -1.34(55) corresponding to the above T-c values, respectively.