The paper describes a new numerical method for the calculation of noncrossing probabilities for arbitrary boundaries by a Poisson process. We find the method to be simple in implementation, quick and efficient - it works reliably for Poisson processes of very high intensity n, up to several thousand. Hence, it can be used to detect unusual features in the finite-sample behaviour of empirical process and trace it down to very high sample sizes. It also can be used as a good approximation for noncrossing probabilities for Brownian motion and Brownian bridge, in particular when the boundaries are not regular. As a numerical example we demonstrate the divergence of normalized Kolmogorov-Smirnov statistics from their prescribed limiting distributions (Eicker (1979), Jaeshke (1979)) for quite large n in contrast to very regular behaviour of statistics of Mason (1983). For the Brownian motion case we considered square-root, Daniels' (1969) and Grooneboom's (1989) boundaries.