find close points
ind = find(S1G,x,epsilon) % find all points in distance epsilon ind = find(S1G,x) % find closest point