Loop infinito

eletric_rj
(usa Outra)

Enviado em 15/08/2021 - 14:19h

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) 

leandropscardua
(usa Ubuntu)

Enviado em 15/08/2021 - 14:48h

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