close all;
clear;

f = 1;
w = 2 * pi * f;
t = linspace(0, 3 / f, 100);
Ui = sin(w * t);
L = 1 * 10 ^ (-4);
C = 22 * 10 ^ (-6);

filter = [];

for f=1:10000
    w = 2 * pi * f;
    Ui = sin(w * t);
    Uo = (Ui * w ^ 2 * L * C) / (w ^ 2 * L * C - 1);

%plot(t, Ui);
%hold on
%plot(t, Uo);
%legend('Ui', 'Uo');

    filter = [filter, max(Uo) / max(Ui)];
end

plot(dB(filter))