How to create a matrix with for-loop when the columns are unknown?

I want to find the intersections of multiple points (known) to a function when the number of intersections for each point is unknown. For one point I have:

intersections=find(abs(y_function-y_point)<=(0.001));

how do I create a matrix, that includes the intersections for n points?

Thank you in advance

Answer 1

[Andrei Bobrov]

Let y_point - vector [n x 1]

[ii,jj] = find(abs(bsxfun(@minus,y_function(:),y_point(:)')) < 1e-3);

[Theodosios Stergiou]

Thanks Andrei. How can I gather the results into different n arrays?