exb += ((a * np.cos((i+1) * t * np.pi / 12))