First-order necessary conditions for quadratically optimal reduced-order modeling of linear time-invariant systems are derived in the form of a pair of modified Lyapunov equations coupled by an oblique projection which determines the optimal reduced-order model. This form of the necessary conditions considerably simplifies previous results of Wilson [1] and clearly demonstrates the quadratic extremality and nonoptimality of the balancing method of Moore [2]. The possible existence of multiple solutions of the optimal projection equations is demonstrated and a relaxation-type algorithm is proposed for computing these local extrema. A component-cost analysis of the model-error criterion similar to the approach of Skelton [3] is utilized at each iteration to direct the algorithm to the global minimum.