underoverflow euler-cromer SIUM(3.5) A_03 RUNGE KUTTA METHOD O(4) euler-cromer Practice 1(IV) A_01 CHEBYSHEV_METHOD gghhhhh Check identity matrix or not? Gauss backward elimination