In using dynamic programming, by taking only accessible states for the x-grid and using an iterative procedure employing region contraction, only a small number of grid points are required at each iteration to yield very good accuracy even if the dimension of the system is high. The effect of the number of grid points and the choice of the contraction factor are analysed by considering a non-linear system consisting of eight ordinary differential equations and four control variables.