Loop infinito

um programa para achar resultado aproximado de integral pelo método trapézio está em loop infinito e não consigo achar o erro.

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) 

Enquanto c >= 1, f(c) sempre será menor q b pq é p inverso de uma exponencial, ou então diminua o valor de b