MATLAB program for the design and implementation of Butterworth low pass filter

Program code

clc;
clear all;
close all;
fp=input('Please enter the first pass band frequency: ');
fs=input('Please enter the first stop band frequency: ');
rp=input('Please enter the pass band attenuation: ');
rs=input('Please enter the stop band attenuation: ');
sf=input('Please enter the sampling frequency: ');
wp=2*fp/sf;
ws=2*fs/sf;
[n,wn]=buttord(wp,ws,rp,rs);
[b,a]=butter(n,wn,'low');
freqz(b,a);
title('Butterworth Low Pass Filter');

Example of output

Please enter the first pass band frequency: 2000

Please enter the first stop band frequency: 1600

Please enter the pass band attenuation: .5

Please enter the stop band attenuation: 53

Please enter the sampling frequency: 9500

__________________________________