TBILISI MATHEMATICAL JOURNAL, cilt.14, sa.2, ss.83-91, 2021 (ESCI)
In this research, we discuss the construction of analytic solution of non-homogenous initial boundary value problem including PDEs of fractional order. By means of separation of variables method, the solution is constructed in the form of a Fourier series with respect to the eigen-functions of a corresponding Sturm-Liouville eigenvalue problem including fractional derivative in Liouville-Caputo sense.