24小时热门版块排行榜

>论坛更新日志 (3315)
>虫友互识 (349)
>导师招生 (201)
>文献求助 (197)
>硕博家园 (107)
>公派出国 (90)
>考博 (85)
>找工作 (71)
>论文投稿 (69)
>博后之家 (68)
>基金申请 (64)
>论文道贺祈福 (60)
>考研 (60)
>休闲灌水 (60)
>教师之家 (51)
>招聘信息布告栏 (38)

返回列表

xmq135828

铁杆木虫 (著名写手)

应助: 96 (初中生)
金币: 323.4
散金: 5804
红花: 48
帖子: 1535
在线: 1694.1小时
虫号: 1267622
注册: 2011-04-16
专业: 无机非金属材料

[求助] 有没有人运行过这样的程序，C语言的

#include
#include
#include
#include
#include
double a,b,c,al,be,ga,sal,sbe,sga,cal,cbe,cga,mm,nn,ss,f1,f2,f3;
/*以下是变量的说明
a，b，c：三个基矢的长度
al，be，ga：三个基矢之间的夹角（角度值）
sal，sbe，sga：三个基矢之间的夹角的正弦值
cal，cbe，cga：三个基矢之间的夹角的余弦值
mm，nn，ss：倒易点阵三个基矢之间的夹角的余弦值
f1，f2，f3：倒易点阵的三个基矢的长度
fvect：矢量类型

以下为函数说明
angle：计算倒易点阵中两个矢量间的夹角，即正空间中两个面的夹角
input：输入参数
*/
struct fvect
{
float x,y,z;
};
struct fvect cef,rif;
float angle(struct fvect aa,struct fvect bb)
{
float se1=aa.x*bb.x*f1*f1+aa.y*bb.y*f2*f2+aa.z*bb.z*f3*f3+(aa.x*bb.y+bb.x*aa.y)*ss+(aa.x*bb.z+bb.x*aa.z)*nn+(aa.y*bb.z+bb.y*aa.z)*mm;
float se2=sqrt(aa.x*aa.x*f1*f1+aa.y*aa.y*f2*f2+aa.z*aa.z*f3*f3+(aa.x*aa.y+aa.x+aa.y)*ss+(aa.x*aa.z+aa.x*aa.z)*nn+(aa.y*aa.z+aa.y*aa.z)*mm);
float se3=sqrt(bb.x*bb.x*f1*f1+bb.y*bb.y*f2*f2+bb.z*bb.z*f3*f3+(bb.x*bb.y+bb.x*bb.y)*ss+(bb.x*bb.z+bb.x*bb.z)*nn+(bb.y*bb.z+bb.y*bb.z)*mm);

float ee=se1/se2/se3;
if(fabs(ee)<0.01) ee=0;
return (ee);
}

void input(void)
{
double pi,angle1;int t;
struct fvect p1,p2,p3;
printf("please input a,b,c\n"

;
scanf("%lf %lf %lf",&a,&b,&c);
printf("please input al,be,ga\n"

;
scanf("%lf %lf %lf",&al,&be,&ga);
pi=3.1415926/180;t=1;
al=al*pi;be=be*pi;ga=ga*pi;
sal=sin(al);sbe=sin(be);sga=sin(ga);
cal=cos(al);cbe=cos(be);cga=cos(ga);
mm=(cal*cbe-cga)/sal/sbe;
nn=(cal*cga-cbe)/sal/sga;
ss=(cbe*cga-cal)/sbe/sga;
f1=a*sal/(sqrt(t-cal*cal-cbe*cbe-cga*cga+cal*cbe*cga+cal*cbe*cga));
f2=b*sbe/(sqrt(t-cal*cal-cbe*cbe-cga*cga+cal*cbe*cga+cal*cbe*cga));
f3=c*sga/(sqrt(t-cal*cal-cbe*cbe-cga*cga+cal*cbe*cga+cal*cbe*cga));
printf("please input the first plane\n"

;
scanf("%lf %lf %lf",&cef.x,&cef.y,&cef.z);
printf("please input the second plane\n"

;
scanf("%lf %lf %lf",&rif.x,&rif.y,&rif.z);
angle1=angle(cef,rif);
angle1=acos(angle1);
printf("The angle between two plane is %f",angle1*180/3.1415926);
}

void main()
{
input();
system("pause"

;
}

各位帮忙看看，这段代码有什么错误啵
运行时总是提示main函数返回值不能为空，改为int后，输入参数后的运行结果如附图

运行结果

回复此楼