An accurate pentadiagonal matrix solution for the time-dependent Schr\"{o}dinger equation