function step=step_pd(A,O,PHI0,draw_plot)

pd = [0:180];
step = zeros(1,length(pd));

N=5;

for i=[1:length(pd)]
  [x,y]=trajectory(A,O,pd(i),PHI0,N,false);
  step(i)=max(x)/N;
endfor;

plot(pd,step);
grid on;
axis([0,180]);
title("Step vs phase difference");
ylabel("Step (mm)");
xlabel("Phase difference (degrees)");




