(i, n, h)
| 1405 | eps = np.sqrt(10**(ripple/10)-1) |
| 1406 | h = np.tanh((1/n)*np.arcsinh(1/eps)) |
| 1407 | def a_i(i, n, h): return np.sqrt( |
| 1408 | 1/(1-h**2) - (np.sin((2*i-1)/n*np.pi/2))**2) |
| 1409 | |
| 1410 | def b_i(i, n, h): return np.sqrt(1 + 1/(h*np.tan((2*i-1)/n*np.pi/2))**2)/2 |
| 1411 | if n % 2: |