Разбиение функции на n равных площадей

from sympy import *
x=Symbol('x')
i=[]
c=[]
 
def integral(func,a,b,step):
    global x
    a1=a
    while a<b:
        if 0.9<=((integrate(func,(x,a1,a)))/(integrate(func,(x,a,b))))<=1.1:
            cd=a
        a=a+step
        a=round(a,len(str(step))-2)
    return cd
 
def integ(func,a,b,step,count):
    global x,i,c
    while len(i)<=count:
        for n in range(1,abs(len(i)-1)):
            c.append(integral(func,i[n-1],i[n],step))
        count+=1
        i=c
    i.append(a)
    i.append(b)
    i.sort()
    return i

from sympy import *
x=Symbol('x')
i=[]
c=[]
 
def integral(func,a,b,step):
    global x
    a1=a
    while a<b:
        if 0.9<=((integrate(func,(x,a1,a)))/(integrate(func,(x,a,b))))<=1.1:
            cd=a
        a=a+step
        a=round(a,len(str(step))-2)
    return cd
 
 
def integ(func,a,b,step,count):
    global x,i,c
  
    while len(i)-2<=count:
        if len(i)<=1:
            i.append(a)
            i.append(b)
        for n in range(1,abs(len(i))):
            c.append(integral(func,i[n-1],i[n],step))
            count+=1
        i=c
    if not(a in i):
        i.append(a)
    if not(b in i):
        i.append(b)
    i.sort()
    return i

Traceback (most recent call last):
  File "<pyshell#20>", line 1, in <module>
    integ(((sin(pi*x)**2)/2)+0.25,1,2,0.005,3)
  File "C:/Python32/full", line 25, in integ
    c.append(integral(func,i[n-1],i[n],step))
  File "C:/Python32/full", line 14, in integral
    return cd
UnboundLocalError: local variable 'cd' referenced before assignment