function dy = lotka_volterra(t,y,hist,dy,alpha,x_buffered,weight)

	x_hist0 = hist(1:2, 1);

	x = y(1:2);
	x_buffered(1) = x_hist0(1);
	x_buffered(2) = x(2);
	x_in = matmul(weight, x_buffered);
	
	dy(1:2) = x.*(alpha + x_in);

end


function z = matmul(x,y)
    z = mtimes(y,x');
end