Summary: We develop a numerical method for computing the semistabilizing solution of a generalized algebraic Riccati equation (GARE). The semistabilizing solution of such a GARE has been used to characterize the solvability of the \((J, J^\prime)\)-spectral factorization problem for general rational matrices which have poles and zeros on the extended imaginary axis. The main difficulty for solving such a GARE is that its associated skew-Hamiltonian/Hamiltonian pencil has eigenvalues on the extended imaginary axis; consequently, it is not clear which eigenspace of the associated skew-Hamiltonian/Hamiltonian pencil can characterize the desired semistabilizing solution; i.e., it is not clear which eigenvectors and principal vectors corresponding to the eigenvalues on the extended imaginary axis should be contained in the eigenspace that we wish to compute, and hence the well-known generalized eigenvalue approach for the classical algebraic Riccati equations cannot be directly employed for it. Our proposed method consists of computations of the eigendecomposition of the system pencil corresponding to the eigenvalues on the extended imaginary axis and the stable eigenspace of an augmented matrix pencil; hence, it is a generalization of the generalized eigenvalue approach for the classical algebraic Riccati equations.