  The polynomial spline collocation method is proposed for solution of Volterra
integral equations of the first kind with special piecewise continuous kernels.
The Gauss-type quadrature formula is used to approximate integrals during the
discretisation of the proposed projection method. The estimate of accuracy of
approximate solution is obtained. Stochastic arithmetics is also used based on
the Contr\^{o}le et Estimation Stochastique des Arrondis de Calculs (CESTAC)
method and the Control of Accuracy and Debugging for Numerical Applications
(CADNA) library. Applying this approach it is possible to find optimal
parameters of the projective method. The numerical examples are included to
illustrate the efficiency of proposed novel collocation method.
