We study a particular case of a convergence model of interest rates. The bond prices are given as solutions of a parabolic partial differential equation and we consider different possibilities of approximating them, using approximate analytical solutions. We consider an approximation already suggested in the literature and compare it with a newly suggested one for which we derive the order of accuracy. Since the two formulae use different approaches and the resulting leading terms of the error depend on different parameter sets of the model, we propose their combination, which has a higher order of accuracy. Finally, we propose one more approach, which leads to higher accuracy of the resulting approximation formula.