Вычисление квадратного корня без использования стандартных функций

        public static double SqrtByGeron(double value)
        {
            return GeronFunction(value, 3.0, 1, 1000);
        }
 
        public static double GeronFunction(double number, double limit, int step, int count)
        {
            double answer = 0.0;
 
            for (int iteration = step; iteration <= count; ++iteration)
            {
                answer = (limit + (number / limit)) / 2.0;
 
                double diff = (answer - limit) - 1.0;
                if (diff == -1.0 || diff == 1.0)
                {
                    break;
                }
 
                limit = answer;
            }
 
            return answer;
        }

        // вызываем так: double result = SqrtByGeron(25.0, 1E-4);
 
        public static double SqrtByGeron(double value, double delta)
        {
            return GeronFunction(value, 3.0, delta, 1, 1000);
        }
 
        public static double GeronFunction(double number, double limit, double delta, int step, int count)
        {
            double answer = 0.0;
 
            for (int iteration = step; iteration <= count; ++iteration)
            {
                answer = (limit + (number / limit)) / 2;
 
                double diff = (answer - limit);
                if (
                    ((diff > 0.0) && (diff <= delta)) ||
                    ((diff < 0.0) && (-diff <= delta))
                   )
                {
                    break;
                }
 
                limit = answer;
            }
 
            return answer;
        }

float Sqrt(float x)
{
    float xhalf = 0.5f * x;
    int i = *(int*)&x;            // evil floating point bit level hacking  
    i = 0x5f375a86 - (i >> 1);    // what the fuck?  
    x = *(float*)&i;                
    x = 1 / (x * (1.5f - xhalf * x * x)); // 1st iteration
    return x;
}