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

# include 
# include 
# include 
# include 

# define N   40

double th0=2.182;
double ph0=1.658;
double l=4.0;

int main(void)
{
  int i;
  double x[N], y[N], z[N];
  double xm, ym, zm;
  FILE *fp;

  ym=-cos(ph0)*sin(th0)*cos(th0/2)+sin(th0/2)-cos(th0)*sin(th0/2);
  ym=ym*l;

  xm=(l*cos(ph0)*sin(th0)+ym*cos(th0/2))/sin(th0/2);
  xm=-xm;

  zm=l*l-xm*xm-l*l*cos(th0/2)*cos(th0/2)-2*xm*l*cos(th0/2)-ym*ym-l*l*sin(th0/2)*sin(th0/2)+2*ym*l*sin(th0/2);
  zm=sqrt(zm);

  for(i=0;i<N;i++)
  {
      x[i]=0;
      y[i]=i*4;
      z[i]=-y[i]*(ym-2*l*sin(th0/2))/zm;
  }
  fp=fopen("snp0","w");
  for(i=0;i<N;i++)
  { 
    fprintf(fp,"%.3f    %.3f    %.3f\n",x[i],y[i],z[i]); 
  }
  fclose(fp);
  return 0;
}

# include 
# include 
# include 
# include 

# define N   40

double th0=2.182;
double ph0=1.658;
double l=4.0;

int main(void)
{
  int i;
  double x[N], y[N], z[N];
  double xm, ym, zm;
  FILE *fp;

  ym=-cos(ph0)*sin(th0)*cos(th0/2)+sin(th0/2)-cos(th0)*sin(th0/2);
  ym=ym*l;

  xm=(l*cos(ph0)*sin(th0)+ym*cos(th0/2))/sin(th0/2);
  xm=-xm;

  zm=l*l-xm*xm-l*l*cos(th0/2)*cos(th0/2)-2*xm*l*cos(th0/2)-ym*ym-l*l*sin(th0/2)*sin(th0/2)+2*ym*l*sin(th0/2);
  zm=sqrt(zm);

  for(i=0;i<N;i++)
  {
      x[i]=0;
      y[i]=i*4;
      z[i]=y[i]/(ym-2*l*sin(th0/2))*zm;
  }
  fp=fopen("snp0","w");
  for(i=0;i<N;i++)
  { 
    fprintf(fp,"%.3f    %.3f    %.3f\n",x[i],y[i],z[i]); 
  }
  fclose(fp);
  return 0;
}