function [Ydot] = Gompertz(t,Y)
Ydot(1) = Y(1)*Y(2);
Ydot(2) = -Y(2);
Ydot=Ydot';