sysdiag — block diagonal system connection
r=sysdiag(a1,a2,...,an)
Returns the block-diagonal system made with subsystems put in the main diagonal
subsystems (i.e. gains, or linear systems in state-space or transfer form)
Used in particular for system interconnections.