Within the framework of a three-dimensional Ising model with nearest magnetic neighbor interactions, and using a Monte Carlo–Metropolis algorithm for equilibrium and energy minimization, we compute the canonical ensemble averages for magnetization, magnetic susceptibility and specific heat of stoichiometric magnetite Fe3O4. In the model, atoms are distributed on an inverse spinel structure for which interactions arising from antiferromagnetic Fe3+A–Fe3+A, Fe3+A–Fe3+B, Fe3+A–Fe2+B couplings and ferromagnetic Fe3+B–Fe3+B, Fe3+B–Fe2+B, Fe2+B–Fe2+B bonds are taken into account. Labels A and B refer to tetrahedral and octahedral sites, respectively. On the basis of finite-size scaling theory, our best estimates of critical exponents for the correlation length, specific heat, magnetization and susceptibility are ν=0.65(2), α=0.23(2), β=0.35(1) and γ=1.13(2), respectively. Simulation results as a function of temperature also allow obtaining information of the contributions per magnetic ion to the total magnetization and to the susceptibility in a differentiated fashion.