Computation of limit cycle via higher order harmonic balance approximation and its application to a 3-bus power system