In this paper, the time-fractional heat equation with the Caputo derivative of order α where 0 < α ≤ 1 is considered. The parametric Crank-Nicholson type method for direct problems is used. But for ...