如何在C ++中找到20位精度的平方根?

时间:2016-12-06 16:13:22

标签: c++ math precision

嗨,伙计们,我要坚持这个功课,我需要找到方程的根源,使用 Bisection方法,精度为10 ^ -20又名0.00000000000000000001所以起初我虽然这是因为我不是'使用long double和L在数字的末尾,但是即使我使用它时我的最后3位数字也不正确,因为下面给出的代码要求你给出我的情况下的数字是5,所以我得到2.3227751229355622087

虽然正确的答案应该是2.3227751229355622988,我真的找不到我的错误,如果有人帮我解决这个问题会很高兴。

供您参考,以下是Bisection method 的说明和说明。

这是我的代码:

#include<iostream>
#include<cmath>
#include<math.h>
#include<iomanip>

using namespace std;

long double f(long double   x, long double a);
long double F = 123456L % 100L;

long double f(long double  x, long double a)
{
    long double  sum = pow(x, 5) - a*x - F;
    return sum;
}

int main()
{
    cout.setf(ios::fixed);
    long double a, b, c, fa, fb, fc;
    long double e;
    long double aa;
    bool flag = true;
    while (cin >> aa)
    {
        cout.precision(19);
        flag = true;
        a = 0L;
        b = 10L;
        e = 0.00000000000000000001L;

        if (f(a, aa)*f(b, aa)>0)
        {
            flag = false;
        }

        while(fabs(a-b)>=e){
            c = (a + b) / 2.0L;
            fa = f(a, aa);
            fb = f(b, aa);
            fc = f(c, aa);

            if (fc == 0)
            {       
                break;
            }

            if (fa*fc>0)
            {
                a = c;
            }
            else if (fa*fc<0)
            {
                b = c;
            }
        } 

        if (flag == true)
        {
            cout  << c << endl;
        }
        else 
        {
            cout << "NO SOLUTION" << endl;
        }
    }
    return 0;
 }

1 个答案:

答案 0 :(得分:-2)

问题是我使用了错误的编译器(Visual Studio)。

一旦安装了另一个编译器,20位精度就没有任何问题。在我进行更多调查时,似乎有些编译器在小数点后有一个限制为14或15。 如果c ++拒绝将数字舍入到更高的精度,则更改编译器:)