clear
a=8; b=2; c=5; d=3; dt=.001;
x=input('x di partenza? ');
y=input('y di partenza? ');
nmax=input('quante iterate? ');
n=input('quale figura? ');
figure(n)
set(n,'position',[407 12 612 684])
fl=input('cancella figura?  ','s');
if fl=='y'
  clf
end
axis([0 5 0 10])
grid on;
hold on;
xlabel('Prede','fontsize',15,'color','r')
ylabel('Predatori','fontsize',15,'color','r')
fl=input('auto? ','s');
passo=input('passo?  ');
set(n,'defaultlinelinewidth',2)
plot(x,y,'r*')
for i=1:nmax
  title(['Iterazioni=' num2str(i)],'fontsize',15,'color','r')
  if fl=='y'
    if ceil(i/passo)==i/passo
      pause(.001)
    end
  else
    if ceil(i/passo)==i/passo
      pause
    end
  end
  xx=x+dt*(a*x-b*x*y);
  yy=y+dt*(-c*y+d*x*y);
  plot(x,y,'b*')
  plot(xx,yy,'r*')
  x=xx;
  y=yy;
end
