The motivation of this study is to determine the analytic solution of initial
boundary value problem including time fractional differential equation with Neumann
boundary conditions in one dimension. By making use of seperation of variables, the
solution is constructed in the form of a Fourier series with respect to the eigenfunctions
of a corresponding Sturm-Liouville eigenvalue problem.