Four numerical techniques based on the linear B-spline functions are presented for the numerical solution of the Lane-Emden equation. Some properties of the B-spline functions are presented and are utilized to reduce the solution of the Lane-Emden equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the new techniques.Several authors have investigated this equation. The interested reader can see [1,[6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25]. However, the approach in the current paper is different. This paper deals with the numerical solution of the Lane-Emden equation [1] by means of the B-spline techniques where either the first or second derivatives of the solution are approximated by B-splines. In the current paper, we give four different numerical methods to solve the Lane-Emden equation (1.1). In the present methods, we reduce the problem to a set of algebraic equations by expanding the unknown function as the linear B-spline functions specially constructed on bounded intervals with unknown coefficients. Some properties of the linear B-spline function like the operational matrices of the derivative and integration are given. These matrices together with the B-spline functions are then utilized to evaluate the unknown coefficients.Also, there have been numerous applications of B-spline to the numerical solution of singular initial and/or boundary value problems. But we will not review them here. Also, we refer the interested reader to [26][27][28] for alternative approaches in the numerical solution of initial or boundary value problems [29][30][31][32][33][34][35][36][37][38].This article is organized as follows: In Section 2, we describe the formulation of the B-spline functions on OE0, 1, the operational matrices of the derivative and integration and some properties required for our subsequent development. In Section 3, the present techniques are proposed to approximate the solution of the given problem. As a result, a set of algebraic equations is formed, and a solution of the considered problem is introduced. In Section 4, we report our computational results and demonstrate the accuracy of the proposed numerical schemes by presenting several test problems. Section 5 completes this paper with a brief conclusion.