#include
#include
#define maxsize 100
#define eps 0.000001
int main(void)
{
    double a,b,t,x[maxsize],step,s0,s1;
    double f(double m);
    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++;
        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;
}
    这是数值分析或者计算方法里面：数值积分与数值微分里面的题目，求积分的近似值 。把不规则的图形换成梯形，然后将梯形分割成若干个小梯形，使小梯形相加，最后近似的求出不规则图形的面积

   

你这没有注释说明很不好理解，只能给你提点思路：
你既然把不规则图形分成n个小梯形，你就应求出x0，x1对应的y0，y1，y0 y1就是两条边的长，x1-x0就是高，梯形面积=（顶边+底边）*高/2，把n个小梯形面积加起来就是近似值

直接用积分求解不行吗？

楼主代码有点小问题
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是没有赋值的，为随机值，所以会运行出错。