Miscellany
public text v1 · immutables=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;
ends=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