#include?
#include?
#define?N?8???//N可自定義，越大對應最后的結果精度越高，這里最大只能取8，可精確到小數點后面8位以上，可以參考matlab的計算結果（matlab積分精度特別高）
double?gass_integral（double?（*）（double）doubledoubleint）;//高斯求積分函數的聲明

void?main（）?
{??
	double?function_name（double）;//自定義函數的聲明
	double?a=1b=3;//a，b分別表示積分函數的下、上限，可自行更改
	double?ans;//最終的結果
	ans=gass_integral（function_nameabN）;//調用高斯求積分函數
	printf（“ans=%.10f\n“ans）;//輸出結果，顯示小數點后10位
?}?

double?gass_integral（double?（*func）（double?x）?double?a?double?b?int?n?）//高斯求積分函數，勿更改（如果自定義函數的輸入參數有多個，double?x這塊還是要改的。。。）
?{
	//高斯點及其求積系數列表，話說我也是蠻拼的，眼睛都看疼了，才輸完------------------------------------------------------------------
	double??x1[1]={0.0};
	double??A1[1]={2};
	double??x2[2]={-0.55735026920.5573502692};
	double??A2[2]={11};?
	double??x3[3]={-0.77459666920.00.7745966692};
	double??A3[3]={0.55555555560.88888888890.5555555556};
	double??x4[4]={0.3399810436-0.33998104360.8611363116-0.8611363116};
	double??A4[4]={0.65214515490.65214515490.34785484510.3478548451};
	double??x5[5]={0.00.5384693101-0.53846931010.9061798459-0.9061798459};
	double??A5[5]={0.56888888890.47862867050.47862867050.23692688510.2369268851};
	double??x6[6]={0.2386191861-0.23861918610.6612