#include
#include
float?y（float?x）
{
return?sqrt（1+pow（3*pow（x2）+6*x+12））;
}
float?fsimp（float?afloat?bfloat?（*p）（float））
{
float?cs;
c=（a+b）/2;
s=（b-a）/6*（p（a）+4*p（c