function yex=yexact(t)
yex=exp(-t); % yex=exact solution