#ifndef SOLVE_H #define SOLVE_H double t1_solve (const double x_0, const int n, const double * restrict X, const double * restrict Y); #endif