* np.cos(2. * np.pi * y / sampling * np.cos(f_angle)