Peak/Bell/Parametric filter on Matlab

Hello everybody.

I am looking for a function which creates a peak/bell/parametric filter with configurable gain and frequency.

I am looking for something like this...

[B,A]= iirnotch (Wo,BW); [N,d] = freqz(B,A,2*length(w),'whole');

but with the kind of filters I am looking for instead of notch filters.

----------------------------------------------------------------------------------

iirpeak function is not usefull because it does not create bell filters. http://www.mathworks.es/es/help/dsp/ref/iirpeak.html

-----------------------------------------------------------------------------------

The filter I am looking for should look like this:

https://dl.dropbox.com/u/27700839/peak.png