4.4 Numerical implementations