In DSP math, it is common to use Chebyshev polynomial approximation. You can get incredibly precise results within your required bounds.