In this work, we give an adaptive grid generation method which allows a single point to be added in the regions of large variation. This method uses a quadrature rule as a weight function. Our weight function measures the variation of the solution function on each subinterval of the solution domain. The method is applied to obtain the numerical solutions of some differential equations. A comparison of the numerical solution obtained by this method and other methods is given.