Watch this Author's files
I think there is an error in the code: besselh(n,1,k*radius) is used instead of the besselh(n,2,k*radius) (Hankel function of 2nd kind)
Contact us at files@mathworks.com