function LVexpeuler2(x,y,T,N)
% x, y: nyul -ill. rokapopulacio merete kezdetben
% T - ido
% N - idolepesek szama a numerikus megoldas soran
h=T/N;
xvektor=x; 
yvektor=y; 
for n=1:N 
    u=f(x,y); 
    v=g(x,y); 
    x=x+h*u; 
    y=y+h*v; 
    xvektor=[xvektor,x]; 
    yvektor=[yvektor,y]; 
end 
t=0:h:T; 
plot(t,xvektor,'r', t, yvektor,'b') 
xlabel('idõ')
ylabel('préda (piros), ragadozó (kék)') 
function U=f(x,y)
U=2*x-x.*y;     % ha a = 2 és b = 1 
function V=g(x,y)
V=-1*y+1*x.*y; 