polynomial system solver via eigen-computations