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