function[]=gershbands(G)
switch class(G)
　　case 'tf'　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　% Transfer Function Form
　　　　G_TF=G; [num,den]=tfdata(G);　　　　　　　　　　　　　　　　　　　　　　　
　　　　 for i=1:size(num,1)　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　　　　　for j=1:size(num,2)　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　
　　　　　　　　num_symbolic(i,j)=poly2sym(num{i,j},'s'); den_symbolic(i,j)=poly2sym(den{i,j},'s');
　　　　　　end
　　　　end