On a novel approach for the investigation and approximation of solutions to the systems of higher order nonlinear PDEs