function fazna_slika_duffing(a,z_0)
  %nariše fazno sliko Duffingovega oscilatorja

  F=@(z,t)[z(2); z(1)-z(1)^3-a*z(2)]; %desne strani DE
  n=size(z_0,2);
  t=linspace(0,10);
  for i=1:n
    z=lsode(F, z_0(i,:),t); %rešimo sistem z danimi začetnimi pogoji
    plot(z(1,:),z(2,:));
    hold on;
  endfor

  hold off;

endfunction