This work is licensed under Creative Commons Attribution 4.0 License GJES.MS.ID.000647.of the value between 0.3315 and 0.33175 and later the value was determined with sufficient precision as f′′(0) =0.332057336. using, among other algorithms, the shooting method and its well recommended variants for BVPs [5][6][7][8][9][10][11][12] as well as the finite differences and collocation methods. Parallel to the continuous appearence of new schemes and families of numerical methods for solving ODEs [13][14][15][16], ubiquitous stiff problems are emerging namely in electronics, solid mechanics, plasma physics, chemical