§Example:  Consider a transfer matrix
•
•             
•
• with
•
•
• Then the command h2norm(Gs) gives ||Gs||2= 0.6055 and h2norm(cjt(Gu)) gives ||Gu||2= 3.182. Hence            
•                                                            
• >> P = gram(A,B); Q = gram(A´,C´); or P = lyap(A,B*B´);
• >> [Gs,Gu] = sdecomp(G); % decompose into stable and             antistable parts.