function dydt=lv(t,y)
a=2;b=1;c=0.01;d=0.01;
dydt=[a*y(1) - c*y(1)*y(2); -b*y(2) + d*y(1)*y(2)];