Automatic computation of polyharmonic balance equations for non-linear differential systems