Numerical solution of the two-point boundary value problems arising in modeling viscoelastic flows

Shooting methods will, be used to compute solutions to the primary two-point boundary value problem arising in modeling viscoelastic flow. The resulting initial value problems will be reguralized and then solved using fifth order Runge-Kutta Fehlberg scheme. Due to the presence of the singularity, such scheme will have at least h(2)-order of convergence. Comparison with the work of others will also be presented through some numerical results and examples.