The filter was tested by computing the convolution of the input x(n) with h(n).
Using a sampled sine wave at frequency=f_0=2000 Hz:
Using a sampled sine wave at frequency = f_0+B/2 = 2250 Hz:
The results compare exactly with the results obtained by the method of difference equations in question 8.
The Matlab Code from 10 can be seen here.