Lagrange's error

perturbation theory in celestial mechanics involved very long and technical calculations, with many successive canonical transforms: imagine that there’s an error with a probability \(p\), and the calculation involves \(n\) steps: when would the probability for an error be larger than \(0.5\) as a function of \(p\) and \(n\)?