Factorial estimate for Picard iteration ($|y_{n+1}(x)-y_n(x)|\le M{L}^n|x-a{|}^{n+1}/(n+1)!$)

For Picard iteration $y_{n+1}(x)=b+{\int }_a^{x}f(s,y_n(s))ds$, $y_0=b$, assume $|f|\le M$ and a Lipschitz constant $L$ on a rectangle containing every iterate. The first difference is at most $M|x-a|$. Integrating the preceding bound multiplied by $L$ proves the displayed estimate by mathematical induction. Its uniform bounds sum to a convergent exponential series for every finite $Lh$, so the Weierstrass M-test gives uniform convergence without requiring $Lh<1$.

For two arbitrary initial functions $y,z$ in a function set preserved by the Picard iteration operator $T$, the same induction givesThus some iterate is a contraction mapping on any fixed finite interval for which every iterate stays in that set. The iterated contraction argument gives a unique fixed point of $T$ on a nonempty complete metric space, without requiring $Lh<1$. For a rectangle with field bound $M$ and vertical radius $k$, $Mh\le k$ supplies the required self-mapping property.