
#include
#include
#include
main()
{ void I1();
void I2();
void I3();
I1();
I2();
I3();
}
void I1()
{ double f(double x);
int N,i;
double a,b,deltax,I1=0;
printf("a,b,N\n");
scanf("%lf%lf%d",&a,&b,&N);
deltax=(b-a)/N;
for(i=0;i I1=I1+deltax*f(a+deltax*i);
printf("\n输出I1值\n");
printf("%lf\n",I1);
}
void I2()
{ double f(double x);
int N,i;
double a,b,deltax,I2=0;
printf("a,b,N\n");
scanf("%lf%lf%d",&a,&b,&N);
deltax=(b-a)/N;
I2=I2+0.5*deltax*(f(a)+f(b));
for(i=1;i I2=I2+deltax*f(a+deltax*i);
printf("\n输出I2值\n");
printf("%lf\n",I2);
}
void I3()
{ double f(double x);
int N,i;
double a,b,deltax,I3=0;
printf("a,b,N\n");
scanf("%lf%lf%d",&a,&b,&N);
deltax=(b-a)/N;
I3=I3+1.0/3*deltax*(f(a)+f(b)); //有改动
for(i=1;(2*i-1) I3=I3+4.0/3*deltax*f(a+deltax*(2*i-1)); //有改动
for(i=1;(2*i) I3=I3+2.0/3*deltax*f(a+deltax*(2*i)); //有改动
printf("\n输出I3值\n");
printf("%lf\n",I3);
}
double f(double x)
{ double y;
y=cos(1/(1+x*x));
return (y);
}