A numerical method for solving differential equations in which the dependent variables are expanded as series of orthogonal basis functions and the original equations are reduced to a set of algebraic equations for the modal amplitudes.