求救!!!!错误提示:exp overflow error,

我的程序如下,不知道哪出问题,郁闷

 

#include “glaguerr.c”
#define Cq 0.1377
#define Q 1464
#define j 1281.75
#define k 1
#define q 240
#define L 8
extern double glaguerre(int n);
double f(double x);
main()
{
int i,n;
double integral;
n=8;
system(“cls”);
printf(“Gauss-Laguerre method:\n”);
for(i=2;i<n;i++)
{
integral=glaguerre(i);
printf(“n=%2d,integral=%16.8f\n”,i,integral);
}
getch();
}
double f(double x)
{
double y,E,U;
E=1/(sqrt(2*3.1415)*Cq*Q);
U=2*Cq*Cq*Q;
y=E*exp(-(x-Q)*(x-Q)/U)*(1-pow(1-j*exp((q-x)/k),L))*(1/(k*j))*exp((x-q)/j)/exp(-x);
return(y);
}