Zahl ^ 0.5 = Wurzel aus Zahl

float Potenziere(float fZahl, int iPotenz)
{
    int ix;
    float ferg = 1;
    for(ix = 0; ix < iPotenz; ix++)
    {
        ferg = ferg * fZahl;
    }
    return ferg;
}

y = x² | - x²
y - x² = 0

x_(i + 1) = x_i - f(x_i) / f'(x_i)

x_(i + 1) = x_i - (y - x_i²) / (-2x_i)

Local y# = 25
Local xi#, xi1# = y

Const EPS# = 0.000001


Repeat
   xi = xi1
   Print xi
   
   xi1 = xi - (y - xi*xi) / (-2 * xi)
Until Abs(xi - xi1) <= EPS

WaitKey

float Potenziere(float fZahl, int iPotenz) //Potenzieren Funktion
{
    int ix; //Zähler
    float ferg = 1; //Ergebnis
    for(ix = 0; ix < iPotenz; ix++) //Solange Zähler kleiner als Potenz
    {
        ferg = ferg * fZahl; //Rechne Ergebnis mal Zahl
    }
    return ferg; //Ergebnis returnen
}

float Wurzel(float fZahl) //Wurzel Funktion
{
    float fStart = 1; //Startwert
    float fErg = 0; //Ergebnis
    const int iGenauigkeit = 5; //je größer desto genauer die berechnung
    int ix; //Zähler
    fStart = (Potenziere(fStart, 2) + fZahl) / (2 * Potenziere(fStart, 1)); //1. Berechnung.
    for(ix = 0; ix < iGenauigkeit; ix++) //Nach 5 durchläufen (6 Rechnungen) Genau genug ;)
    {
        fStart = (Potenziere(fStart, 2) + fZahl) / (2 * Potenziere(fStart, 1)); //Berechnen
    }
    fErg = fStart; //Zuweisung
    return fErg; //return
}

xi = (xi^2 + xi) / (2 * xi^1)

xi = (xi^2 + radicand) / (2 * xi^1)

xi = (xi + radicand/xi) /2

float Wurzel(float fZahl) //Wurzel Funktion
{
    float fStart = 1; //Startwert
    const int iGenauigkeit = 5; //je größer desto genauer die berechnung
    int ix; //Zähler
    for(ix = 0; ix < iGenauigkeit; ix++)
    {
        fStart =  (fStart + fZahl/fStart) / 2
    }
    return fStart; //return Ergebnis
}

#include <stdio.h>

float fsqr( float n ) {
   float result;
   __asm__("fld %0;"
         "fsqrt;"
         "fstp %1;"
          :: "f" (n), "m" (result));
   return result;
}

int main() {
   printf("%f", fsqr(2.0));
   return 0;
}

double fastpow(double x, double y)
{
   double a = 1.0 - y;
   double b = x - 1.0;
   double c = 0.5 * b *b;
   double d = b - c
   return x /((1.0 + (a * d)) + (0.5* (d * d * a*a));   
}