Optimal feedback control system design of a nuclear power plant has usually been conducted in the time domain by solving the Riccati equation. In this paper, the problem will be treated in the s-domain by means of the MacFarlane extension criterion (MEC). The algorithm is based on generalized spectral factorization (GSF), and the MEC is resolved by the Newton-Raphson method featuring rapid convergence and high precision. Whether the plant is stable or unstable, the initial matrix can be chosen effectively and a general computer code for finding the optimal feedback gain matrix of the multi-input system is presented. The control system design of a boiling water reactor (BWR) shows that the results obtained from the proposed method are consistent with those from solving the conventional Riccati equation. A satisfactory closed-loop multivariable system with a prescribed degree of stability can also be designed in the s-domain. The methodology of this paper can also be generalized to treat similar control sys...