Full code review

Convergence should be tested with an epsilon, not with a max number of iterations. Or maybe, even better, you can add both options. Include comments for each constant that you define. And include all these constants in a separate input file please.