Approximate Crank–Nicolson Algorithm with Higher-Order PML Implementation for Plasma Simulation in Open Region Problems. (10th February 2021)