code with parameters from the numerical example