Hi Magdy,
For a change of pace, here is a function that calculates the coefficents of Hn by recursion. I didn't want to get into variable precision arithmetic but with double precision it's good to somewhere around n = 25. Usual Matlab notation so e.g.
ans = 64 0 -480 0 720 0 -120
means 64x^6 - 480x^4 + 720x^2 - 120
H = 2*[A 0] - [0 0 (n-1:-1:1).*A(1:end-1)];