R1=3.3e+3;
C1=0.01e-6;
C2=0.01e-6;


for R2=[1000 6000 100e+3]
H=tf([1 0 0], [1 (R1*(C1+C2))/(R1*R2*C1*C2) 1/(R1*R2*C1*C2)])
%w=logspace(0, 3)
%bode(H)
h=bodeplot(H);
setoptions(h,'FreqUnits','Hz','PhaseVisible', 'off');

grid on
hold on
end