function [y, yp] = ftosol(x)
% y gives the function value at x, while yp gives the derivative at x:
y = sin(x) - exp(-x);
yp = cos(x) + exp(-x);
end