H = tf([0,1/(R*C),1/(R*C)],[1,1/(R*C),1/(L*C)]);

for n = 0:4

    %Insignal
    A = 1; % amplitud
    freq = 100; % frekvens
    omega = 2*pi*freq*(2*n+1); % omega

    t = linspace(0, 10e-2, 10000); % tidsvektor
    z = 1j*omega; % skicka in komplext tal
    eval = evalfr(H,z);  % evaluerar tf för valt z

    %Utsignal
    utsignal = abs(eval)*A*sin(omega*t + angle(eval));
    plot(t, utsignal)
    hold on

end