#include // // ODE definition // double f(double x, double y) { // return y+x-1; return x-y; // return exp(x); }