“…Following the ideas in [4], re-formulate the problem (1) either by introducing a physical parameter or identify a physical parameter that appears either in the differential equation or in the boundary conditions as independent variable. Let us introduce λ as physical parameter in differential equation 1, so we have obtained, We wish to determine the approximate numerical solution u(x) of the problem (2) for different values of λ and λ 0 ≤ λ ≤ λ M . We define M − 1 numbers of nodal points in [λ 0 , λ M ], as λ 0 < λ < λ 2 < ...... < λ M using a uniform step length H such that λ i = λ 0 + jH, j = 0, 1, .., M. Also, we define N − 1 numbers of nodal points in [a,b], the domain in which the solution of the problem (1) is desired, as a ≤ x 0 < x 1 < x 2 < ...... < x N = b using a uniform step length h such that x i = x 0 + ih, i = 0, 1, .., N. To simplify the expression we denote the numerical approximation of u(x) at the node x = x i as u i .…”