Abstract:
Here is a replacement for invhilb.m with the following features: 1. Fast (only one level of for-loop) 2. Robust against overflow: only returns Inf when the entry concerned cannot be represented in IEEE double precision 3. Accurate (produces closest machine number to actual value)