rmsGain.m

function g = rmsGain(H,f1,f2,N)
%function g = rmsGain(H,f1,f2,N=100)
%Compute the root mean-square gain of the discrete-time
%tf H in the frequency band (f1,f2)

if nargin<4
    N = 100;
end
w = linspace(2*pi*f1,2*pi*f2,N);
g = norm( evalTF(H,exp(j*w)) ) / sqrt(N);