clc

pkg load control
pkg load symbolic

wo_ideal=sqrt(1.567*10^7)
fo_ideal=wo_ideal/(2*pi)
Q_ideal=wo_ideal/1847

s=tf('s');
h_ideal=(1847*s)/(s^2 + 1847 * s + 1.567*10^7);

[y_step_ideal, t_step_ideal] = step(h_ideal);

[mod_bode_ideal, fase_bode_ideal, w_bode_ideal] = bode(h_ideal);

return;