
where \(f^{k+1/2}=(f^{k+1}+f^k )/2\) . However, the CN method may introduce spurious oscillations to the numerical solution for non-smooth data unless the algorithm parameters satisfy the maximum principle [1]}.
