Numerical determination of the critical isotherm of the ferromagnetic Ising model