The combustion problem by the Newton method
int main(
int argc,
char**argv) {
string approx = (argc > 2) ? argv[2] : "P1";
Float tol = (argc > 4) ? atof(argv[4]) : eps;
size_t max_iter = (argc > 5) ? atoi(argv[5]) : 100;
size_t n_iter = max_iter;
return (
residue <= sqrt(tol)) ? 0 : 1;
}