Simulate and plot Butterworth filter frequency response using transfer function

I am attempting to the simulate the frequency response of a butterworth fiilter. The transfer function for the butterworth filter is, R/(s^3(L1)(L2)+S^2RC(L1)+S(L1+L2)+R). The filter is designed to be a bandpass filter at 1kHz.

This is the code I have

R = 10;

L1 = 10;

L2 = 10;

C = 1;

numerator = R;

a = L1*L2*C;

b = R*C*L1;

c = L1+L2;

denominator = [a, b, c, R];

sys = tf(numerator, denominator);

stepplot(sys)

Answer 1

[Star Strider]

If you want to see the frequency response, use the bodeplot function —

R = 10;
L1 = 10;
L2 = 10;
C = 1;                                                                   % Changed
numerator = R;
a = L1*L2*C;
b = R*C*L1;
c = L1+L2;
denominator = [a, b, c, R];
sys = tf(numerator, denominator);

hbp = bodeplot(sys);
grid
setoptions(hbp, 'FreqUnits','kHz');

The component values need to be changed to achieve the desired frequency response. I leave that to you.

.