function [StaticLoopSensitivity,DCgain]=tfanalysis(sys)
% function [StaticLoopSensitivity,DCgain]=tfanalysis(sys)
% sys is any siso lti system
%

tfsys=tf(sys);

num=getfield(tfsys,'num');
den=getfield(tfsys,'den');

j=find(num{:}~=0);
StaticLoopSensitivity=num{:}(1,j(1));

[drow,dcol]=size(den{:});
[nrow,ncol]=size(num{:});
DCgain=num{:}(1,ncol)/den{:}(1,dcol);
