这是我用Ｃ语言编的二元积分，但是运行不出来，有哪位高手能指点一下． 好让我对人生不太失望
| #include stdio.h
#include math.h
#define N 100
double f(x,y)
double x,y;
{
double f;
f=exp(-x*x-y*y);
return(f);
}
double integral(fun,a,b,c,d)
double (*fun)(),a,b,c,d;
{
double h1,h2,v;
int i,j;
v=0.0;
h1=(b-a)/N;
h2=(d-c)/N;
for(i=1;i=N;i++)
for(j=1;j=N;j++)
v=v+(*fun)(a+j*h1,b+i*h2)*h1*h2;
return(v);
}
main()
{
double f(),integral();
double y;
y=integral(f,1,20,3,40);
printf(" The result is: %f\n",y);}
这是我用Ｃ语言编的二元积分，但是运行不出来，有哪位高手能指点一下． 好让我对人生不太失望 |
| 欢迎光临 光电工程师社区 (http://bbs.oecr.com/) | Powered by Discuz! X3.2 |