C语言——一些处理脚本(16)——(待完善)

# include <stdio.h>
# include <ctype.h>
# include <math.h>
# include <string.h>
# include <stdlib.h>

# define N   1000
# define kt    5
# define t0    0.873    
# define pi    3.14

# define sig  4

int main(void)
{
  int i, j, k;
  double  t;
  double f1[N], f2[N], rij;
  double x[N], y[N], y1[N];
  double a1, a2, a6, a10, a12;
  double  eps;

  FILE *fp;

  eps=2*pi/N;


  for(i=0;i<N;i++)
  {
      x[i]=(-pi)+i*eps;
      y[i]=kt*(1-cos(x[i]-t0));//+kt*(1-cos(3*(x[i]-t0)));
  }



  fp=fopen("cos5","w");
  for(i=0;i<N;i++)
  {
     fprintf(fp,"%g  %g  %g 0\n",x[i],((x[i]*180)/pi),y[i]);
  }
  fclose(fp);
  
  return 0;
}