#define h 1.0/n

这是为什么了？？？

define 也就是 直接替换的

f(i*h) = f(i*1/n) = f(i/n)

俩整数，结果就是 0

# include <stdio.h>
# define n 100000000.0
# define h (1/n)
double f(double x)
{
    double y;
    y=1/(1+x*x);
    return y;
}

double area(double len1,double len2,double high)
{
    return (len1 + len2) * high / 2;
    // （上底+下底）*高/2
}



void main()
{
    int i;
    double f(double x);
    double s=0;
    double h1 = 0;
    double h2 = 0;
    for(i=0;i<n;i++)
    {
        h1 = f(i * h);
        h2 = f(i * h);

        s += area(h1, h2, h);
        // s+=f(i*h);
    }
    // s=(3/4+s)*h;
    printf("%lf\n",s);
}

答案是0.785~~？

[ 本帖最后由 yuccn 于 2013-3-28 19:31 编辑 ]

定积分……还梯形法……[0,1]是上下限吗？好遥远啊……

如果是h=0应该输出0.000000啊，为什么是1.000000

嗯嗯，大侠

你的f函数 f(0) = 0？

# define h 1/n 应该是#include 1.0/n 因为整数除整数结果是整数
s=(3/4+s)*h; 应该是s=(3.0/4+s)*h; 原因与1/n相同