“…When solving time dependent partial differential equations (PDEs), it is common to first discretize the spatial derivatives to form a system of ordinary differential equations (ODEs), and then solve the system of ODEs by time-stepping methods, such as Runge-Kutta methods and linear multistep methods (LMMs). This is the so-called method of lines (MOL) and is efficient for many time dependent PDEs [4,5,6,7,8,9,10,11,12,13,14]. When the PDEs to be solved are restricted in frequently encountered bounded domains, it is necessary to properly treat boundary conditions so that the overall accuracy and stability are maintained.…”