PYTHON: We want to solve the differential equation df/dx = 2 + exp(x), f(0) = 1in the interval [0, 1].

Divide the interval into N divisions, x₀ = 0, x₁ = h = 1/N, x₂ = 2h = 2/N, , x_N = Nh/N = 1, and use df/dx = [f(x+h) f(x)]/h to determine the value of f(x) at x₁, x₂, , x_N.

Calculate the error: assume the exact solution is F(x), and define the error as sum over i from 0 to N of |f(x_i) F(x_i)|/N.

Do the calculation for N = 10, 20, 30, , 200 and plot the error as a function of N.