整理后积分如下:

以此表达式编程——

代码文本:
#include "stdio.h"
#include "math.h"
int main(int argc,char *argv[]){
double s,x,t,y,dx=0.000001,dy=0.001;
for(s=x=0;x<=1.0;x+=dx){
for(t=y=0;y<=x;y+=dy)
t+=y;
s+=t*dy/sqrt(1+x*x*x);
}
printf("%g\n",s*dx);
return 0;
}
dx和dy的值取得越小结果越精确,但太小执行时间太长。这里取了个精度和执行时间兼顾的值,大约不到4秒出结果。供参考。