Introduction - If you have any usage issues, please Google them yourself
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