| 查看: 883 | 回复: 4 | ||
[求助]
有没有人运行过这样的程序，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后，输入参数后的运行结果如附图  运行结果 |
回复此楼» 猜你喜欢
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 一个区间搜索的c语言程序 出问题了
已经有8人回复
- 求C语言基本运算的运行时间分析
已经有10人回复
- 【求助】本人对C语言有一点基础，打算学习一VC，请各位推荐几本经典入门教材
已经有17人回复
- 【求助】求将c语言改写为C++，求教高人啊！
已经有6人回复
- 【分享】谭浩强C语言教程讲义，强烈推荐
已经有19人回复
2楼2012-03-12 12:15:02
wydmhy
至尊木虫 (职业作家)
- 应助: 161 (高中生)
- 金币: 10554.3
- 散金: 320
- 红花: 41
- 帖子: 4055
- 在线: 719.6小时
- 虫号: 903420
- 注册: 2009-11-15
- 性别: GG
- 专业: 无机非金属类光电信息与功
3楼2012-03-12 12:29:48
4楼2012-03-12 17:24:32
5楼2012-03-12 17:58:50
投票:
30