变步长梯形公式.cpp

来自「数值计算程序 包括：变步长梯形公式算法、二叉树算法、二分法、高斯列主消去法、曲线」· C++ 代码 · 共 23 行

#include<stdio.h>
#include<math.h>
void main( )
{float a,b,c,s,h,T1,T2,x;

scanf("%f,%f,%f",&a,&b,&c);
h=b-a;
T1=(h/2.0)*(1+sin(b)/b);
   s=0;x=a+(h/2.0);
     while(x<b)
        {s+=sin(x)/x;
           x+=h;}
    T2=1.0/2*T1+(h/2.0)*s;
    if(fabs(T2-T1)<c)
       printf("%f\n",T2);
    else
        {while(fabs(T2-T1)>=c)
     {T1=T2;h=h/2;
     s=0;x=a+(h/2.0);
               while(x<b)
                   {s+=sin(x)/x; x+=h;}
   T2=1.0/2*T1+(h/2.0)*s;}
   printf("%f\n",T2);}}