computation.cpp: rhs, time step (AH)