genera un codigo que resuelva ecuaciones de park para maquinas sincronas

def park(θ, d, q, b, λ, ω, t):
  return np.array([[np.cos(θ), np.sin(θ), 0],
  				  [-np.sin(θ), np.cos(θ), 0],
  				  [0, 0, 1]]) @ np.array([d, q, b]) + λ * np.array([[np.cos(ω * t), np.sin(ω * t), 0],
  				  [-np.sin(ω * t), np.cos(ω * t), 0],
  				  [0, 0, 1]]) @ np.array([0, 0, 1])