function [f] = f_ivp(y,t)
   f = y;
   f(1) = y(2);
   f(2) = -exp(y(1));
end