Numerical Estimation of the Current Large Deviation Function in the Asymmetric Simple Exclusion Process with Open Boundary Conditions