光学编程

这是我用Ｃ语言编的二元积分，但是运行不出来，有哪位高手能指点一下．　好让我对人生不太失望

#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);} 这是我用Ｃ语言编的二元积分，但是运行不出来，有哪位高手能指点一下．　好让我对人生不太失望