贝塞尔函数表达式



以下为程序代码
#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的文件

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