Decription of the file MSD_paper.m

In the file the code is provided to reduce a systems consiting of 200
mass-spring-dampars (in series).

By calculating the Hankel singular values it is clear to which order the 
system could be approximated. In this case 20.

mPoissonMatrix.m makes the standard Poisson matrix.

sSMOR.m contains the construction of the orthosymplectic basis, 
see Algorithm 1 in the manuscript. Furthermeore, it construct the 
reduced model on basis of this. The reduced order is 20.

sSMORH2.m does the Symplectic Model Order Redution method of 
the manuscript but, combined it with the H2-optimal method.

mIRKA.m contains the IRKA-PH algorithm as described in Algorithm 1. 
see p. 1967 of Gugercin et al (2012)

The plots comparing the different approximation methods are constructed 
at the end of the file. 
