C语言 梯形积分法

#include<iostream>
#include<math.h>
using namespace std;
int main(){
 float fun(float x);
 float a,b;
 cout<<"请输入函数X^2的定积分的下限a和上限b：";
 cin>>a>>b;
 int n=50;//将区间划分成50份 
 float h=(b-a)/n;//h是每个区间分大小 
 float s=0;//s是矩形的面积的和
 float i=0;
 for(i=a;i<b;i+=h){
  s=s+fun(i)*h;
 } 
 cout<<"\n结果是:"<<s<<endl;
 cout<<endl; 
}
float fun(float x){
 return pow(x,2);
}