A multicenter numerical integration scheme for polyatomic molecules