Iterasi Integrasi Numerik Metode Trapesium

fx = lambda x : 6*x**2 - 2*x

def integrasiTrapesium(a, b, n) :
  h = (b - a) / n
  s = fx(a) + fx(b)
  i = 1
  while(i < n) :
    s += 2 * fx(a + i * h)
    i += 1
  return (h/2) * s


eksak = 3311
tol =  0.1
err = 1000
xa = 1
xb = 12
n = 1

while(err>tol) :
  hasil = integrasiTrapesium(xa, xb, n)
  err = abs(hasil-eksak)
  n=n+1
  print("hasil : ",hasil," , error : ", err)
  if(n==500) :
    break

print(n)