root2 = np.sqrt((1.0 - kx) * (1.0 - kx) + 4.0 * tan2delta)