bibliography about implementation of IIR filters
Okay, according to the wikipedia article you cited (with an appropriate substitution in notation):

$$ \Big|H(j2\pi f)\Big|^2= 10^{\frac{2.00}{10}} \times \frac{12194^4\cdot f^8}{\big(f^2+20.6^2 \big)^2(f^2+107.7^2)(f^2+737.9^2)\big(f^2+12194^2\big)^2}$$

and the dB amplitude curve is this:

$$A(f)=10\log_{10}\left(\Big|H(j2\pi f)\Big|^2\right) $$

I hope you can see how I modified the language and that it is equivalent.  The purpose is to express this transfer function as an $H(s)$ and get where (and how many) the poles as zeros are.

Now you gotta factor some more shit out.  Oh, this is laborious and painful.  Sorry, unless someone wants to finish this answer, I will have to return to it.