求:y=√x
1.二分法求零点
由于函数在(0,+oo)为单调函数,取值为正
y*y - x = 0 设 g(y) = y*y - x ,函数值在自变量y属于(0,+oo)也单调。g(0)<=0 ,g(x+1)>0 所以可在区间[0,x+1]上采用二分法
#includeconst double eps=1e-10;//精度小数点后10位
//求y=根号x
//则y^2=x
double Abs(double val) {//求val的绝对值
if (val < 0)
val = -val;
return val;
}
int main(void) {
double x;
scanf("%lf", &x);
double l = 0, r = x + 1, mid = l;
while (Abs(r - l) > eps) {//区间长度大于eps就继续二分区间
mid = (l + r) / 2;
if (mid * mid > x)
r = mid;
else
l = mid;
}
printf("%.20f", mid);
return 0;
}
二分法的时间复杂度为 log2(x/精度)
eg.
输入:2
输出:1.41421356229693630000
2.三