Today, in many different sciences, including engineering and physics, it is important to obtain the solutions of differential equations, among which the most important of these equations is partial differential equations. In most cases, there is no analytical answer to these equations or it is very difficult to obtain them. That is why obtaining a suitable numerical method is an important issue to find the approximate solution of these equations. Parabolic equations are a kind of equations of partial derivatives. In recent years, the numerical solution of time dependent separation parabolic equations has been of interest to researchers by using high-order splitting methods with complex coefficients.
