梯形法积分.cpp

数值计算的源代码

#include<stdio.h>
#include<math.h>
float A,B;

float y(float x)
{return A*x*x+B;}

void main()
{float a,b,x,h,sum,trapz;
int i,n;
printf("请输入抛物线y=Ax(2)+B的系数A,B:");
scanf("%f%f",&A,&B);
printf("抛物线为y=%fx(2)+%f\n",A,B);
printf("请输入积分的上下限a,b的值:");
scanf("%f%f",&a,&b);
printf("请输入积分区间的等份n的值:");
scanf("%d",&n);

h=(b-a)/n;
sum=(y(a)+y(b))/2;
for(i=1;i<=n-1;i++)
{x=a+i*h;
sum=sum+y(x);
}


trapz=sum*h;
printf("该抛物线与x轴所围成的面积:%f",fabs(trapz));
}