b = np.pi/2