double *ab = expansion(x0, d, alfaTab[j], Nmax);

				exp << x0 << ";" << solution::f_calls << ";" << ab[0] << ";" << ab[1] << endl;
				cout << "x0: " << x0 << "\t ab[0]: " << ab[0] << "\t" << "ab[1]: " << ab[1] << " N: " << solution::f_calls << endl;
				solution::clear_calls();

				solution f = fib(ab[0], ab[1], epsilon);
				fibb << x0 << ";" << f.x << ";" << f.y << ";" << solution::f_calls << endl;
				cout << x0 << "," << alfaTab[j] << "," << f.x << f.y << solution::f_calls << endl;
				solution::clear_calls();

				solution l = lag(ab[0], ab[1], epsilon, Nmax);
				lagg << x0 << ";" << l.x << ";" << l.y << ";" << solution::f_calls << endl;
				cout << x0 << "," << alfaTab[j] << "," << l.x << l.y << solution::f_calls << endl;
				solution::clear_calls();