Почему не выполняется код в цикле while (Метод Симпсона)

from math import sqrt


def f(x):
    return 1/sqrt(1+x**4)

a=0
b=1
n=100

h = (b - a) / n
s = f(a) + f(b)
for i in range(1, n, 2):
    s += 4 * f(a + i * h)
for i in range(2, n-1, 2):
    s += 2 * f(a + i * h)
print(s * h / 3)

eps = 0.0001
h1 = h
s1 = s
s2 = 0

while abs(s1 - s2) > eps:
    h /= 2
    s2 = 0
    for i in range(0, n, 2):
        s2 += f(a + i * h)
    for i in range(1, n, 2):
        s2 += 2 * f(a + i * h)
    s2 = (s + s2) * h / 3    


print("The result is:", s2)