代码如下,但是似乎你的题目有问题:
根号X乘lnX 0到1,根号0乘ln0等于0,根号1乘ln1等于1,怎么求?
#include
#include
#include
double fsimpf(double x) /*要进行计算的被积函数*/
{
x05double y;
y=log(x)*sqrt(x);
return(y);
}
double fsimp(double a,double b,double eps,int n) /*辛普森算法:a为积分下限,b为积分上限,eps是希望达到的精度*/
{
int k;
double h,t1,t2,s1,s2,ep,p,x;
h=(float)(b-a)/n;
t1=h*(fsimpf(a)+fsimpf(b))/2.0; /*用梯形公式求出一个大概的估值*/
s1=t1;
ep=eps+1.0;
while (ep>=eps)
x05{
x05x05/*用梯形法则计算*/
x05x05p=0.0;
x05x05for (k=0;k