Defining the Model

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

> rate_eqn2 := diff(F(t),t)=(0.00004)*R(t)*F(t)-(0.04)*F(t);

> MODEL := { rate_eqn1, rate_eqn2 };

> vars := { R(t), F(t) };

> init1 := [R(0)=2000, F(0)= 600];

> init2 := [R(0)=2000, F(0)=1200];

> init3 := [R(0)=2000, F(0)=3000];

> domain := t = 0 .. 320;