【问题标题】:What is wrong with my logic for finding square root by Newton's Method?我用牛顿法求平方根的逻辑有什么问题?
【发布时间】: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.下面是泰勒对代码的回答。您遇到的是整数除法(不允许余数),因此您的某些函数没有返回正确的浮点值。

标签: c algorithm math


【解决方案1】:
double guess(int n)
{
    return  n / (double) 2;
}

【讨论】:

  • 或除以 2.0 而不是 2。
【解决方案2】:

你忘了一个return

double rootByNewtonApprox(int n)
{
if(goodEnough(guess(n)))
    return guess(n);
else
    return rootByNewtonApprox(improve(guess(n)));
    ^
}

【讨论】:

  • 先生,在进行了上述所有更改后,仍然没有输出,只是 segmentation Fault 盯着我的脸!
  • 我认为没有理由重复 J.S.泰勒指出。所以我将原因添加到您的第二个问题 - n=2 返回 nan。
猜你喜欢
  • 1970-01-01
  • 2012-10-02
  • 1970-01-01
  • 2019-06-27
  • 2021-12-09
  • 2014-03-07
  • 1970-01-01
  • 2012-10-24
  • 2013-09-25
相关资源
最近更新 更多