Application of the Pade approximation via Lanczos (PVL) algorithm to electromagnetic systems with expansion at infinity