function f = f691(x,y)
f = [y(1) - 2*y(2) - 2*exp(-x) + 2
        2*y(1) - y(2) - 2*exp(-x) + 1];
%f = [y(1)/y(2)
%      y(1)/2];