Question

我正忙着参加我正在上课的课程。我必须在C ++中使用割线方法找到函数的根。我把它全部写完了，我发现第二个零点为0.146;但是，我的教授正在寻找的根是.064。我已经尝试了所有我能做到的事情，但我无法弄清楚如何使其输出0.064。他告诉我们放入一个包含该函数的头文件，这里是头文件：

    #define FX0 pow(x0, 3) - 0.165*pow(x0, 2) + 3.993E-4
    #define FX1 pow(x1, 3) - 0.165*pow(x1, 2) + 3.993E-4

以及代码：

    #include <iostream>
    #include <cmath>
    #include <iomanip>
    #include "Secant.h"
    #include "Keyboard.h"
    using namespace std;


    int main()
    {
    //Define Variables
    float x0,x1,x2,tolerance,maxIterations,count,FX;
    count = 0;
    x0 = 0.02;
    x1 = 0.05;
    tolerance = .000001;
    maxIterations = 100;
    FX = pow(x0, 3) - 0.165*pow(x0, 2) + 3.993E-4;

    //Loop statement that runs until a Zero is found
    while(fabs(FX0-FX1)>tolerance && count < maxIterations && fabs(FX)>tolerance)
    {
    x2=x1-(FX1*((x0-x1)/(FX0-FX1)));
         FX = pow(x2, 3) - 0.165*pow(x2, 2) + 3.993E-4;

    x0=x1;
    x1=x2;

    count++;
    }


    //Display the zero
    if (fabs(FX)<tolerance)      
    cout << "The zero is at x = " << setprecision(4) << x2;
    //Or Report that no zero was found
    else
    cout << "No zeroes were found within the given function.";

    return 0;
    }

Answer 1

当您使用#define时，编译器只会将您的FX1宏替换为其文本值（在本例中）。所以

 FX1*((x0-x1)/(FX0-FX1))

变为

 pow(x1, 3) - 0.165*pow(x1, 2) + 3.993E-4*((x0-x1)/(pow(x0, 3) - 0.165*pow(x0, 2) + 3.993E-4-pow(x1, 3) - 0.165*pow(x1, 2) + 3.993E-4))

导致错误位置的括号出现问题，而是将3.9993E-4相乘。

尝试在您的定义周围添加括号，或者将它们更改为函数。

#define FX0 (pow(x0, 3) - 0.165*pow(x0, 2) + 3.993E-4)

使用C ++中的割线函数查找根

1 个答案: