The bosonic IIB matrix model, which contains the bosonic part of the IIB matrix model conjectured to be a non-perturbative definition of the type IIB superstring theory, is studied using a numerical method. The large N scaling behavior of the model is shown performing a Monte Carlo simulation. The expectation value of the Wilson loop operator is measured and the string tension is estimated. From the numerical results, the prescription of the double scaling limit is obtained.