%Sólo hay que introducir los valores de      psi, wn y k

w=logspace(-1,3,1000); % vector de frecuencias

%     Apuntes.CoyaN.es

psi=0.354;        %coeficiente de amortiguamiento
wn=104.977;       %frecuencia natural
k=1.15;           %ganancia


	num=[k*wn^2];
	den=[1 2*psi*wn wn^2];

	figure(1);
	hold on
	bode(num,den,w);
	grid on;
	title('Diagrama de Bode del sistema de 2º orden');
	legend('Apuntes.CoyaN.es');
	hold off

	figure(2);
	hold on
    step(tf(num,den));
    grid on;
	title('Diagrama de Bode del sistema de 2º orden');
	legend('Apuntes.CoyaN.es');
	hold off