Question

所以我一直在尝试使用过时的turbo c ++编译器（使用过时的turbo c ++编译器，使用simpsons 1/3规则）来实现一个小程序来计算双积分（是的，我们仍然使用它！）但遗憾的是我＆＃39我似乎遇到了扼杀故障。每当我尝试在下面的程序中添加数组ax [i]的元素给我一个零。这些元素似乎完美地出现在自身中，但添加它们会因某些奇怪的原因而导致零。这是一段代码。

#include<iostream.h>
#include<conio.h>
#include<math.h>

double f(double, double);

void main()
{
    clrscr();
    double a, b, c, d, h, k, ans, z[10][10], ax[10];
    int i, j, nx, ny;
    clrscr();
    a=0.0; b=1.0; c=0.0; d=1.8; h=0.25; k=0.30;
    nx = (b-a)/h;
    ny = (d-c)/k;

    for(i=0; i<10; i++)
    {
        for(j=0; j<10;j++)
            z[i][j] = 0;
        ax[i] = 0;
    }

    //Generating the table
    for(i=0; i<=nx; i++)
    {
        for(j=0; j<=ny; j++)
        {
            z[i][j] = f(a+i*h, c+j*k);
        }
    }

    for(i=0; i<=nx; i++)
    {
        ax[i] = 0.0;

        for(j=0; j<=ny; j++)
        {
            if(j==0 || j==ny)
                ax[i] += z[i][j];
            else if(j%2==1)
                ax[i] += 4*z[i][j];
            else
                ax[i] += 2*z[i][j];
        }
        ax[i] *= k/3.0;
        cout<<ax[i]<<endl;
    }

    ans=  1.0;
    //for(int q=0; q<=nx; q++)
        ans= (ax[0])+(ax[1])+(ax[2]);


    cout<<"Value of integral is: "<<ans;
    getch();
}

double f(double x, double y)
{
    float r;
    r = 2*x*y/sqrt(x*x+y*y);
    return r;
}

想法？

Answer 1

double f(double x, double y)
{
    float r;
    r = 2 * x*y / sqrt(x*x + y*y);
    return r;
}

在x=0和y=0时除以零，因此ax[0]等于nan。当您添加数字时，您将获得nan作为输出。您需要使用L'Hopital的规则来评估该限制，并以不同方式处理第一个元素。

尝试添加数组元素给出零

1 个答案: