function f = predpreyfunc(t,x), f = [x(1)*(1-0.5*x(2)); x(2)*(-0.75+0.25*x(1))];