function [M] = sol(t)
%
%   this function returns M = solution y(t) = e^t
%
%
%      M =  exp(t);

% or the function form example 3 on page 107:

% the exact solution is y(t)=-7/4+0.5*t+11*exp(2*t)/4

% M = -7/4+0.5*t+11*exp(2*t)/4;

% or the function form example 1 on page 455:

% the exact solution is y(t)= (19/16)*exp(4*t)+0.25*t-3/16

 M = (19/16)*exp(4*t)+0.25*t-3/16;

end