Many researchers apply Legendre and Chebyshev wavelet operational matrices(of derivatives or of integration) to solve singular boundary value problems. At x=0(Independent variable) also the solutions are given for y (Dependent variable) for Lane-Emden type equations in many research articles.How the singularity is overcome?