A Boundary Integral Method For Modeling Axisymmetric Flow Around A Rising Bubble In A Vertical Tube And Accurate Numerical Evaluation Of Orthogonal Polynomials