function [Gw ] = Bode2( G, w ) Gw = 0*w; for i=1:length(w) Gw(i) = evalfr(G, j*w(i)); end end