function z = F_sys(x,y,t)

% This is coded up for the system:
%     dx/dt = x + 4 y
%     dy/dt = x + y
%
z(1) = x+4*y;
z(2) = x+y;