This is the first part of the algorithm; it computes the values of the points (virtual sources) that form a line of the image and returns a Nx3 matrix with the ordered points. The three dimensions represent our Euclidean space, as our beamformer is constructed to compute 3D values of space. This first part of the algorithm depends heavily on the starting point of investigation, investigation depth, the sampling density and thus the step between virtual sources we are keen to accept. The bigger the matrix of points which represent a line of our image in fact is, the higher the number of operations are going to be required from the following steps except for focusing.