clc
clear
close all
s = tf('s')
G = 10/((1+0.1*s)^2*(1+0.001*s))
G = zpk(G)
G.DisplayFormat = 'Frequency'
%bode(G)

%tipo 0
%C = 1
%C = (5*(s+1))/s
%C = 0.5/s
C = (s+1)/s
L = C*G

%BW > 100, wc > 65

margin(L)

W = feedback(L,1)