function rabbits=f(t,y,a)
rabbits=[2*y(1)-a(1)*y(1)*y(2);-y(2)+a(1)*y(1)*y(2)];
