Python-сообщество

wn = []                         # empty array for coefficients
M = 25                            # just a number
N = 2*M + 1                      # number of coefficients
r = 10**(-3)                      # sidelobe attenuation
x0 = cosh(1./2*M)*pow(cosh(1./r), -1)    # stop band frequency
for n in range(-M, M+1):
   S = 0; theta_n = 2*pi*n*N
   for m in range(1, M+1):
      theta_m = 2*pi*m/N
      S += cheb(2*M, x0*cos(theta_m/2.))*cos(m*theta_n)
   wn.append((1./N)*(1+2*r*S))

def cheb(q, x): 
	if q == 0:
		return 1
	if q == 1:
		return x
	if q > 1:
		return 2*x*cheb(q-1, x) - cheb(q-2, x)

r = 10**(-3)

sys.float_info(max=1.7976931348623157e+308...