Le he estado dando vueltas y para cubrir todo ese rango con buena precisión habría que construir un polinomio interpolante considerablemente grande,y esto conllevaría bastante más cantidad de instrucciones y de carga para el micro de lo que ya te genera la función sqrt,por lo que no veo viable intentarlo por este método.
A ver si a alguien se le ocurre otra manera