Defining the Model

> rate_eqn1 := diff(h(t),t)=(0.1)*h(t)-(0.005)*h(t)*(1/60)*u(t);

> rate_eqn2 := diff(u(t),t)=(0.00004)*h(t)*u(t)-(0.04)*u(t);

> MODEL := { rate_eqn1, rate_eqn2 };

> vars := { h(t), u(t) };

> init1 := [h(0)=2000,u(0)=600];

> init2 := [h(0)=2000,u(0)=1200];

> init3 := [h(0)=2000, u(0)=3000];

> domain := t = 0 .. 320;