This MATLAB function designs an order n lowpass digital Butterworth filter with normalized cutoff frequency Wn.
www.mathworks.com