【发布时间】:2011-04-30 06:11:47
【问题描述】:
我编写了以下代码,用于通过牛顿法通过逐次逼近求平方根,但它没有给我正确的答案。有人可以解释一下吗?
#include<stdio.h>
#include<stdlib.h>
#define square(x) x*x
double rootByNewtonApprox(int n);
double improve(double n);
double average(double a,double b);
int goodEnough(double guess);
double guess(int n);
int number;
int main(void)
{
double root;
printf("\nEnter the number you want square root of: ");
scanf("%d",&number);
if(number<0)
number = -1* number;
root = rootByNewtonApprox(number);
printf("\nThe square root of %d is %lf\n",number,root);
return 0;
}
double guess(int n)
{
return n/2;
}
double rootByNewtonApprox(int n)
{
if(goodEnough(guess(n)))
return guess(n);
else
rootByNewtonApprox(improve(guess(n)));
}
double improve(double guess)
{
return average(guess,(number/guess));
}
double average(double a,double b)
{
return ((a+b)/2);
}
int goodEnough(double guess)
{
if(abs(square(guess) - number) <= 0.001)
return 1;
else
return 0;
}
现在,当我给n = 2 时,它给出了输出nan,当我给n = 9 时,它告诉segmentation Fault。
【问题讨论】:
-
见@J.S.下面是泰勒对代码的回答。您遇到的是整数除法(不允许余数),因此您的某些函数没有返回正确的浮点值。