 reduced order model (Ar,Br,Cr,Dr) for an original
  state-space representation (A,B,C,D) by using the stochastic
  balancing approach in conjunction with the square-root or
  the balancing-free square-root Balance & Truncate (B&T)
  or Singular Perturbation Approximation (SPA) model reduction
  methods for the ALPHA-stable part of the system.

</PRE>
<A name="Specification"><B><FONT SIZE="+1">Specification</FONT></B></A>
<PRE>
      SUBROUTINE AB09HD( DICO, JOB, EQUIL, ORDSEL, N, M, P, NR, ALPHA,
     $   