# include <stdio.h>
# include <math.h>
# include <string.h>
# include <stdlib.h>
# define N 10
double th0=2.182;
double ph0=1.658;
double l=4.0;
int main(void)
{
int i,j;
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++)
{
for(j=0;j<N;j++)
{
x[i]=i;
y[i]=j;
z[i]=zm*y[i]/(ym-2*l*sin(th0/2));
}
}*/
fp=fopen("snp1","w");
for(i=0;i<N;i++)
{
for(j=0;j<N;j++)
{
x[i]=i*4;
y[i]=j*4;
z[i]=zm*y[i]/(ym-2*l*sin(th0/2));
fprintf(fp,"%.3f %.3f %.3f\n",x[i],y[i],z[i]);
}
}
fclose(fp);
return 0;
}
C语言——一些处理脚本(12)——(待完善)
最新推荐文章于 2024-07-18 11:11:00 发布