Orthogonal decomposition of derivatives and antiderivatives for easy evaluation of extended Gram matrix

Simple and efficient algorithms for orthogonal decomposition of derivatives and antiderivatives of a function with rational Laplace transform are proposed. Based on a new theorem related to Routh αα–ββ expansion, they enable direct evaluation of the extended Gram matrix which has proved to be very useful in model-reduction applications.