a = d*np.sin(wt+delta) + q*np.cos(wt+delta) + z