x = linspace(-2*pi, 2*pi, 60);
y = sin(x).^2./(x+eps);
plot(x,y)
title('Plot of sin(x)^2/x')
[a,b] = ginput(8)  % get up to 8 points
hold on
plot(a,b,'co')     % plot the collected points
hold off