f=logspace(0,5,10^3); 	 

w=2*pi*f; 		
s=j*w;

Rup=480e3;
Rz=3.3e3;
Cz=0.033e-6;
Rp=330;
Rfb=110e3;

Z1=Rup+1./(1/Rz+s*Cz)+Rp;
Z2=Rfb;

KEA=-Z2./Z1;

magKEA2dB=20*log10(abs(KEA)); 
phKEA2deg= angle(KEA)*180/pi;		

hold on;
figure(1);
semilogx(f,magKEA2dB,'linewidth',3);
ylabel('Gain, dB');
xlabel('Frequency, Hz');
title('KEA');
grid on;
figure(2);
semilogx(f,phKEA2deg,'linewidth',3);
ylabel('Phase, deg');
xlabel('Frequency, Hz');
title('KEA');
grid on;
%-----------------------------------------------------------------------