%pylab inline
Populating the interactive namespace from numpy and matplotlib
Rewrite it as a system of first order ODEs. Find the region of absolute stability. (Hint: find the eigenvalue decomposition of the coefficient matrix, and the region of absolute stability is the domain where the step size is chosen so that each eigen value satisfies the condition for the 1D ODE.)