The function computes the balanced realization for the stable system .

As output you obtain the balanced realization

with .

The balancing proceeds as follows:

- First, the function computes the -norm, of the stable system .
- With
- These solutions can now be used to perform a standard balancing.