梯形积分是这么算的吗？我记得是把每个区间当做梯形来算的，那么你的s0是怎么回事？一上来就算一个大的梯形？

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