function dydt=odefun3(t,y)

a=1;
b=.3;
d=2;

r=y(1);
f=y(2);

drdt=a*r-b*r*f;
dfdt=-d*f+b*r*f;

dydt=[drdt; dfdt];