Loop infito

[Danado_Rj]

o programa é para achar o valor aproximado de uma integral pelo método  trapézio , porém  o mesmo está em loop infinito.

intervalo = input("Entre com o  número de intervalos : ")
a = 1
b = 3
n = float (intervalo)
c = float (a)
soma = 0
# FUNÇÃO : 1/x^2
def f(x):
  return 1/(x)**2
h = math.fabs (b-a)/n
c += h
def tpz(f,a,b,c,n):
  while f(c) < b):
    soma = (f(c) + f(c + h))
    c += h
  soma *= 2
  soma = ((f(a) + f(b)) + soma) * h/2
d = tpz(f,a,b,c,n)
print("O resultado da integral é aproximadamente :",d)