The accuracy of a numerical fission gas release algorithm developed by Forsberg and Massih for solving the problem of diffusive flow to a spherical grain boundary is analysed. Estimates of numerical errors are derived for both steady-state and time varying conditions. We also present a method through which the accuracy of the algorithm can be improved or optimised for most applications.