Finite difference method for solving the Schrödinger equation with band nonparabolicity in mid-infrared quantum cascade lasers