贝塞尔函数表达式



以下为程序代码
#include<stdio.h>
#include<math.h>
void main()
{
  
  long jiecheng(int i);
  int n,k;
  double x,jn;
  
  FILE *fp;
  fp=fopen("D:\\my project\\c\\0031.dat","w");
  jn=0.0;
  printf("in put n");
  scanf("\n%d",&n);
  for(x=0.0;x<15.0;x=x+0.1)
  {
    for(k=0;k<15;k++)
    {
      jn=jn+pow(-1,k)*pow(x/2,(2*k+n))/jiecheng(k)/jiecheng((n+k));  
    }
    fprintf(fp,"%f,%f\n",x,jn);
    printf("\n%f,%f\n",x,jn);
    jn=0.0;
  }
   fclose(fp);
}


   long jiecheng(int i)
   {
     long int p=1;
     int n;
     if(i==0)
         p=1;
     else
     {
     for(n=1;n<=i;n++)
     {
       p=p*n;
     }
     }
     return(p);
   }

scanf("\n%d",&n);改成
scanf("%d",&n);

#include<stdio.h>
#include<math.h>
void main()
{

 

  long jiecheng(int i);
  int n,k;
  double x,jn;

 

  //FILE *fp;//你这里为什么要打开文件呀？你把这些注释调能行啊。
  //fp=fopen("D:\\my project\\c\\0031.dat","w");
  jn=0.0;
  printf("in put n");
  scanf("\n%d",&n);
  for(x=0.0;x<15.0;x=x+0.1)
  {
    for(k=0;k<15;k++)
    {
      jn=jn+pow(-1,k)*pow(x/2,(2*k+n))/jiecheng(k)/jiecheng((n+k)); 

    }
    //fprintf(fp,"%f,%f\n",x,jn);
    printf("\n%f,%f\n",x,jn);
    jn=0.0;
  }
   //fclose(fp);
}


   long jiecheng(int i)
   {
     long int p=1;
     int n;
     if(i==0)
         p=1;
     else
     {
     for(n=1;n<=i;n++)
     {
       p=p*n;
     }
     }
     return(p);
   } 

我是想把数据导入一个.dat的文件

大神们，刚才我都试了一下，数据跟原来还是一样的

阶乘计算溢出了。在32位编译器里 long int 和 int 是一样的。所以你的 jiecheng 只能正确计算0-12，13及以后的结果都是错误的。

另外，直接用公式计算受计算机计算精度的影响太大效率也很低下，调整一下计算顺序可以很有效的降低误差提高运算效率。

5分，就提示这么多了。

我就是来看看贝塞尔函数的公式。

这个我坚决不看   坚决抵制  恶心数学方面的  什么 什么 什么的  太费脑子