Computation of balanced realisations from transfer functions