#include<stdio.h>
#include<math.h>
#define maxsize 100
#define eps 0.000001
int main(void)
{
    double a,b,t,x[maxsize],step,s0,s1,f();
    int i,n;

    printf("please input a= and b=\n");
    scanf("%lf  %lf",&a,&b);

    n=1;
    s0=(f(a)+f(b))*(b-a)/2;
    do
    {
            n=n+1;
            step=(b-a)/n;
        for(i=0;i<n;i++)
        {
        
            x[i]=a+i*step;
            s1=0;
                for(i=1;i<=n;i++)
               
                    s1+=(f(x[i+1])+f(x[i]))*step/2;
               
                    t=s0;
                    s0=s1;
               
                    
        
        }
        
    }
    while (fabs(s0-t)>eps);
    {
        printf("s0=%lf\n",s0);
        printf("此时步长为:%d",n);
        printf("\n");
    }
    return 0;
}
double f(double m)
{
    double y;
    y=sqrt(1+pow(m,2));
    return y;
}
   

   

直接给我那些分，告诉你你是怎么错的
你在用在 s0=(f(a)+f(b))*(b-a)/2;没有在开头对f(a)和f(b)即函数f(x)进行声明，改进的方法：
 在   “double a,b,t,x[maxsize],step,s0,s1,f();
    int i,n;”的前面写进 “double f(double m);”即可，赶快试试吧，要是我说错了就不要你给我分

for(i=0;i<n;i++)
        {
            x[i]=a+i*step;
            s1=0;
                for(i=1;i<=n;i++)
                    s1+=(f(x[i+1])+f(x[i]))*step/2;
                     t=s0;
                    s0=s1;

这个循环中的赋值要达到什么目的，第二个for后面是否应该有{}？

确实没见过这样做函数声明的，我测试了几次，发现今天我的编译器抽风，不管被调用的函数在什么位置都不用声明，我狂晕，以前不声明函数就报错的代码也不用声明了。
楼主你能说说你的代码是做什么的吗？

一楼，你说的double f(double m)我试过了没用，也是一样的，我这个函数是用来：将积分区间分成N个等份，用梯形公式计算机积分f(1,0)sqrt(1+pow(x,2))的近似值

这是数值分析或者计算方法里面：数值积分与数值微分里面的题目，求积分的近似值

你先试试矩形法，看看结果

#include <stdio.h>
#include <math.h>
#define N 200

double f(double m);
double trapz(double lower,double upper);
int main(void)
{
    double a,b,s;
   

    printf("please input a= and b=\n");
    scanf("%lf  %lf",&a,&b);
    s=trapz(a,b);
    printf("\nthe result is:%lf",s);
    getch();
    return 0;
}
double trapz(double lower,double upper)
{
int i;
double step;
static double s=0;
step = (upper-lower)/N;
for(i=0;i<N;i++)
s+=(f(lower+step*i)+f(lower+step*(i+1)))*step/2;
return s;
}
double f(double m)
{
    double y;
    y = m*m;
    return y;
}

这个梯形法，我用N=1和2，测试过了，正常；而且N越大越接近2.333，跟积分计算的精确值7/3结果一致

[ 本帖最后由 pauljames 于 2012-11-22 20:17 编辑 ]

楼主代码有点小问题
do
    {
        n++;
        step=(b-a)/n;
        for(i=0;i<=n;i++)
        {
            x[i]=a+i*step;
            s1=0;
         
            for(i=1;i<=n;i++)
                s1+=(f(x[i+1])+f(x[i]))*step/2;
            t=s0;
            s0=s1;
         
        }
    }
这里面内层for循环运行到最后的时候x[i+1]超出输入设定的界限了，另外比如n=20对于某一个外层for循环i=5，内层for循环x[i]中i从6到20是没有赋值的，为随机值，所以会运行出错。