מצאתי את זה ברשת:

#include <stdio.h>
#include <math.h>


const double epsilon=0.0000001;


double sqr2(double x)
{
	double a0=0;
	double a1=1;
	while (fabs(a0-a1)>=epsilon)			// |ai - ai+1| < epsilon
	{
		a0=a1;
		a1=(x/a0+a0)/2;				// ai+1= (x/ai+ai) /2
	}
	return a1;


}

הבעיה היחידה היא שנעשה שימוש בmath.h על אף שרק בערך מוחלט..

אין ממש בעיה לכתוב פונקציה ערך מוחלט אם זאת הבעיה.

פתרון דומה היה כאן כבר בתראד אחר.