For y=F(t,y)y'=F(t,y), y(t0)=y0y(t_0)=y_0, Picard iteration starts from a trial curve, often y(0)(t)=y0y^{(0)}(t)=y_0, and sets

y(n+1)(t)=y0+t0tF(s,y(n)(s))ds.y^{(n+1)}(t)=y_0+\int_{t_0}^tF(s,y^{(n)}(s))\,ds.

This is fixed-point iteration for the .

Convergence mechanism

On a time interval of length at most hh from t0t_0, a state Lipschitz bound LL gives TyTzsuphLyzsup\|Ty-Tz\|_{\sup}\le hL\|y-z\|_{\sup}. One must also keep the iterates in a ball where FF is defined and bounded. The makes both requirements precise. Repeated time integration also gives factorial bounds, which can prove convergence without a one-step contraction on a longer interval.