function yp = lv_ode(t,y)
% LV model for ODE
yp = [10*y(1) - 0.01*y(1)*y(2);
     -10*y(2) + 0.01*y(1)*y(2)];
end

