發送短信 : Similarity solution and Runge Kutta method to a thermal boundary layer model at the entrance region of a circular tube: the Leveque approximation)