First-order probabilistic harmonic power flow