
function zs = func(t)
%-------------------
a = .25 * t;
b = -3/16.0;
c = 19.0 * exp(4*t)/16.0;
zs = a + b + c;

