Wavelets play an essential part in numerical analysis. In this study, a novel numerical technique to solve fractional differential equations (FDEs) corresponding to initial conditions is presented using Haar wavelet approximations. Haar wavelet is first presented with an operational matrix of fractional order integration. Then, illustrative examples are presented to signify the validity and applicability of the proposed method.