Question

Problem 1 The midpoint rule for approximating an integral can be expressed as ["flwydcm) x n s(a – žn+ in where h = (b -a)/n.

Answer 1

https://trinket.io/python3/9020deab0c

Please note this is the function only as asked in question and in order to run it you have to pass f,a,b and n as arguments in it

import numpy as np
def midpointnsum(f,a,b,n=1000):
h=float(b-a)/float(n)
i=np.arange(1,n+1);
x=a-0.5*h+i*h;
y=h*np.sum(f(x));
return y;