th0 = np.arctan(val0)/alpha