An efficient technique for solving polynomial systems with a particular structure is presented. This structure is very specific but arises naturally when computing the critical points of a symmetric polynomial energy function. This novel numerical solution method is based on homotopy continuation and is a particular application of results due to Canny and Rojas. An illustrative example from magnetism is presented.