function y=u(x)
% Exact sol (from Maple session)
% u.m
y=x.*exp(x).*(1-x);                 % ex2
%y=exp(-5*x.^2+5*x-5/4)-exp(-5/4);   % ex1
