function y=exact_sol(x)

y=2*exp(x)+x-1;  % solution of y'=x+y, y(0)=1