Numerical calculation of the incomplete Toronto function