This paper presents the results of computer simulations carried out to determine coordination numbers for a system of parallel cylindrical fibres distributed at random in a circular matrix according to twodimensional pattern created by random sequential addition scheme. Two different methods to calculate coordination number were utilized and compared. The first method was based on integration of pair distribution function. The second method was the modified sequential analysis. The calculations following from ensemble average approach revealed that these two methods give very close results for the same neighbourhood area irrespective of the wide range of radii used for calculation.