Miscellany

s=zpk('s');
Fs = SERD + s*SERE;
for m=1:Norder,
	prodVar = 2 * SERC(m) * SERA(m);
	for n = 1:Norder
		if (n~=m)
			prodVar = prodVar * (s*s - SERA(n)^2);
		end
	end 
	Fs = Fs + prodVar;
end