This is only for first order differential equations. So x0 and x1 are not specified, f(0) is assumed to be C0