Курсовая работа: Особые свойства Гамма-функции Эйлера

#define CN 8

static double cof[CN]={

2.5066282746310005,

1.0000000000190015,

76.18009172947146,

-86.50532032941677,

24.01409824083091,

-1.231739572450155,

0.1208650973866179e-2,

-0.5395239384953e-5,

};

double GammLn(double x) {

double lg,lg1;

lg1=log(cof[0]*(cof[1]+cof[2]/(x+1)+cof[3]/(x+2)+cof[4]/(x+3)+cof[5]/(x+4)+cof[6]/(x+5)+cof[7]/(x+6))/x);

lg=(x+0.5)*log(x+5.5)-(x+5.5)+lg1;

return lg;

}

double Gamma(double x) {

return(exp(GammLn(x)));

}

void main()

{

double x[8],g[8];

int i,j;

clrscr();

cout<<"vvedite x[1]";

cin>>x[1];

printf("\n\t\t\t_________________________________________");

printf("\n\t\t\t| x |Gamma(x) |");