plot_integration_error function

Plot numerical integration error