Исправить ошибку в коде к задаче

Доброго времени суток, форумчане! Написал код для решения рекуррентной формулы, но он не работает. Выдает ошибку ZeroDivisionError: float division by zero в 48 строке программы. Код:
import math as b
Nmax = 10
MAX_COUNT = 2000
y = [0.0]*Nmax
i = 0
k = 0
f = 0
a = 1.0
d = 1.0
c = 1.0
s_i = 0.0
p = 0
l = 0.0
a=0.0
c=0.0
sums = 0.0
Func = 0.0

n = int(input('Введите количество значений Х '))
g = int(input('Введите количество значений точности '))
while (g > Nmax) or (n < 1) or (n>Nmax) or (g < 1):
n = int(input('Введите количество значений Х заново '))
g = int(input('Введите количество значений точности заново '))
x = [0.0]*n
for i in range(n):
x[i] = float(input('Введите Х '))
for i in range(g):
y[i] = float(input('Введите точность'))


for f in range(g):
p = b.ceil(b.fabs(b.log(y[f])/b.log(10)))+1
ps = str(p)
format_1 = 'Точность {0:17.'+ps+'f}'
print(format_1.format(y[f]))
print('\n')
print(' №| x[i] | сумма ряда S(X) |итерации|контр.формула F(X)|контроль точности')
print('\n')
for k in range(n):
Func=(x[k]-b.sin(x[k]))/(x[k]**2)
i = 0
a=x[k]
c=1/(3*2)
s_i = a*c
while not (b.fabs(s_i) < y[f]) and (i != MAX_COUNT):
i+=1
a*=(-1)*(x[k])
c*=1/((2*i+1)*(x[k]**(2*i+1))) (здесь ошибка появляется)
s_i =a*c
sums+=a*c
l = b.fabs(sums - Func)
format_2 = '{0: 2d}|{1:17.'+ps+'f}|{2:17.'+ps+'f}|{3: 8d}|{4:18.'+ps+'f}|{5:17.'+ps+'f}'
print(format_2.format(k+1, x[k], sums, i, Func, l ))
print('\n')
print('\n')

0