clear all

sis1=tf([1],[10^-8 5*10^-5 1])
sis2=tf([1],[10^-8 10^-4 1])
sis3=tf([1],[10^-8 2*10^-4 1])
sis4=tf([1],[10^-8 3*10^-4 1])
sis5=tf([1],[10^-8 5*10^-4 1])

figure(1)
bode(sis1,sis2,sis3,sis4,sis5)
legend('50 Ohm','100 Ohm','200 Ohm','300 Ohm','500 Ohm') 

figure(2)
pzmap(sis1,sis2,sis3,sis4,sis5)
legend('50 Ohm','100 Ohm','200 Ohm','300 Ohm','500 Ohm') 
grid on