add numerical_integration codes