[x,y]=clean(x,y)

function [x,y]=clean(x,y)
% this M-file is for manually removing few stray points from a set of x y
% data
plot(x,y)
N=input('How many stray spikes need to be cleaned??    ')
[a,b]=ginput(2*N);
for i=1:2:2*N-1
    k=find(x>a(i)&x<a(i+1));
    x(k)=[];
    y(k)=[];
end
plot(x,y)