function yp = greenhouse(t,y)
% Greenhouse DE (rhs)
yt1 = -(13/8)*y(1) +(3/4)*y(2) + 14;
yt2 = (1/4)*y(1) - (1/4)*y(2);
yp = [yt1,yt2]';

end

