Question

我已经工作了几个小时试图弄清楚我做错了什么。我需要做的就是使用牛顿方法从数组表示的多项式中找到一个根。这两个函数（polyval和polyder）似乎给了我正确的答案，我觉得主代码正确地做了Newton的方法。我希望有经验的人可以给我一些建议。

#include <iostream>
#include <cmath>

using namespace std;

float polyval(float*, int, float);
void polyder(float*, int, float*);

int main(void)  {
    int n;
    float x=-1,f;
    float tol=pow(10,-5);

    cout << "Enter polynomial order:" << endl;
    cin >> n;

    float* p=new float[n+1];
    float* dp=new float[n];

    cout << "Enter coefficients, starting with the highest power:" << endl;
    for (int k=0;k<n+1;k++) {
        cin >> p[k];
    }

    polyder(p,n,dp);
    f=polyval(p,n,x);

    while (fabs(f)>tol) {
        x=x-f/polyval(dp,n,x);
        f=polyval(p,n,x);
        cout << x << endl;
        cout << f << endl;
    }

    cout << "A real root is at x= " << x << endl;
    cout << "To verify: p(" << x << ") = " << polyval(p,n,x) << endl;

    return 0;
}

float polyval(float* p, int n, float x) {
    float px;
    px=0;

    for (int k=0;k<n+1;k++) {
        px=px+p[k]*pow(x,n-k);
    }

    return px;
}

void polyder(float* p, int n, float* dp) {
    for(int k=0;k<n;k++) {
        dp[k] = p[k+1] * (k+1);
    }
}

Answer 1

您对polyval(dp,n,x)的来电将超出所分配的dp空间，该空间包含n条，而不是n+1。

我是初学者，试图编写一个找到数组函数根的函数

1 个答案: