求定积分的程序，跟大家分享一下，有错的话请指出，谢谢！

/*题 目：编写利用矩形法计算定积分的通用函数。然后利用它分
别计算以下3个数学函数的定积分：
(1)f(x)=x*x-5*x+1
(2)f(x)=x*x*x+2*x*x-2*x+3
(3)f(x)=x/(2+x*x)
*/
//**********************************************************

#include < stdio.h >

//******求第一个公式的积分值******

double funtion1( double x )
{
double fx ;
fx = x * x - 5 * x + 1;
return fx ;

}

//******求第二个公式的积分值******

double funtion2( double x )
{
double fx ;

fx = x * x * x + 2 * x * x - 2 * x + 3 ;
return fx ;
}

//******求第三个公式的积分值******

double funtion3( double x )
{
double fx ;

fx = x / ( 2 + x * x ) ;
return fx ;
}

// 调用函数求通用积分公式

void def_integral( double x1,double x2,int n , double (*p)( double x ))
{
double y1,y2 ; //积分区间X1、X2对应的函数值
double per_len ,height ;
double area , fx = 0.0 ;
int i ;

per_len = ( x2 - x1 ) / n ; //等分成n分

y1 = (*p)( x1 ) ;

for( i = 1 ; i <= n ; i++ )
{
x2 = x1 + per_len ; //单位区间长度
y2 = (*p)(x2) ;
height = (y1 + y2 ) / 2.0 ; //高
area = per_len * height ; //每一等份长度对应的面积
x1 = x2 ; //下限右移
y1 = y2 ; //右移
fx = fx + area ; //面积之和为积分的值
}

printf("所求积分的值为：/n");
printf("%.3f/n" , fx ) ;

}

//*****主函数**************

void main()
{
double a,b; //a、b分别为积分区间的上下限
int n ; //n为要输入的等份 即 区间分成的份数

int selection ;

printf("请输入区间的等分数：/n");
scanf("%d" , &n );

printf("请输入积分区间的上下限(中间用空格隔开)：/n");
scanf("%lf%lf", &a, &b );

printf("请输入你的选择（1、2、3）：/n");
scanf("%d" , &selection );

if ( selection == 1 )
def_integral( a,b, n , funtion1 );
else if( selection == 2 )
def_integral( a,b, n , funtion2 );
else
def_integral( a,b, n , funtion3 );

}