On page 446 of the C++ edition, near the end of the routine simp2(), there is a line that says:

q0=-a[i][k+1]/a[i][kp];

It should read:

q0=-a[i+1][k+1]/a[i+1][kp];