A Schur vector method to solve higher order Lyapunov equations

A general Schur vector method for solving a partial differential equation which appears in designing nonlinear optimal control laws for linear or nonlinear dynamical systems is proposed. The method reduces the dimensionality problem of the direct method. Numerical experiments show that it also reduces the numerical instability of the eigenvector method. >