function yprime = lv(t,y)
% lotka volterra vector field
r=2;  %coefficient
yprime =[0;0];   %creates a column vector with two elements
yprime(1) =  r*y(1)*(1  -  y(2));
yprime(2) =  y(2)*(y(1) - 1);
