HOUGH transform detection circle written in MATLAB