Numerical implementation