function ax=PlotPoints(X)
% PLOTPOINTS plots the points X and returns the axis handle

clf; hold off;
plot(X(:,1),X(:,2),'o');
ax=axis; hold on;