版块导航
正在加载中...
客户端APP下载
登录
注册
帖子
帖子
用户
本版
应《网络安全法》要求,自2017年10月1日起,未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用,请尽快对帐号进行手机号验证,感谢您的理解与支持!
24小时热门版块排行榜
>
论坛更新日志
(483)
>
虫友互识
(34)
>
休闲灌水
(34)
>
基金申请
(14)
>
仿真模拟
(9)
>
考博
(6)
>
导师招生
(5)
>
晶体
(5)
>
海归之家
(4)
>
论文道贺祈福
(4)
>
公派出国
(4)
>
考研
(4)
>
论文投稿
(4)
>
博后之家
(3)
>
教师之家
(3)
>
硕博家园
(2)
小木虫论坛-学术科研互动平台
»
计算模拟区
»
程序语言
»
C/C++
»
【求助】请高手指出一下此程序有多少错误?
5
1/1
返回列表
查看: 727 | 回复: 12
只看楼主
@他人
存档
新回复提醒
(忽略)
收藏
在APP中查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖
zyj8119
木虫
(著名写手)
应助: 65
(初中生)
贵宾: 0.003
金币: 915.1
帖子: 2936
在线: 1329.4小时
虫号: 664177
[交流]
【求助】请高手指出一下此程序有多少错误?
CODE:
#include
#include
#define M 16.043
#define omiga 0.012
#define Tc 190.6
#define Pc 45.99
#define R 83.14
#define epsilon 1-pow(2,0.5)
#define sigma 1+pow(2,0.5)
#define alpha (1+(0.37464+1.54226*omiga-0.26992*pow(omiga,2)*pow(1-pow((T/Tc),(0.5))),2)
#define a ((pow(0.45724*R,2)*pow(Tc,2)/Pc)*alpha
#define b 0.07779*R*Tc/Pc
#define beta b*P/(R*T)
#define q a/(b*R*T)
int fugacity(float P, int T)
{ int s,I,N;
double phi,f,P1;
s=sizeof(P);
for(int k=0;k<=s;k++){
double z0(k)=M(s,1)=0;
double z1(k)=N(s,1)=1;
if((z0(k)-z1(k)).le.1e-6){
z0(k)=z1(k);
z1(k)=1+beta-q*beta*(z0(k)-beta)/((z0(k)+epsilon*beta)*(z0(k)+sigma*beta));
}
}
I=(1/(sigma-epsilon))*log((z1(k)+sigma*beta)/(z1(k)+epsilon*beta));
phi=exp(z1(k)-1-log(z1(k)-beta)-q*I);
f=phi*P1;
return f;
}
void main()
{
fugacity(1000,298);
}
回复此楼
» 猜你喜欢
药学硕士,第一、第二作者已发表6 篇 SCI,药理方向及相关方向2026年/2027年博士申请
已经有6人回复
一篇MDPI论文改变了学习工作和生活
已经有5人回复
26年博士申请自荐-电催化
已经有3人回复
中国地质大学(北京)博士招生补录,数理学院材料科学与工程专业和材料与化工专业
已经有6人回复
收到国自然专家邀请后几年才会有本子送过来评
已经有4人回复
考博
已经有5人回复
26年申博自荐-计算机视觉
已经有4人回复
药化及相关博士的申请
已经有3人回复
高级回复
» 本主题相关价值贴推荐,对您同样有帮助:
【求助】write()输出疑问?
已经有8人回复
【求助】增塑剂除色
已经有2人回复
【求助】急求高手帮忙!只有cif如何晶体画图(绝对没造假!)
已经有10人回复
【整理】《提问的智慧》——献给那些不会提问的人
已经有78人回复
【求助】关于用液相进行甲醇测试的问题
已经有10人回复
【求助/交流】一封老外教授给我的回信~(8.4更新)!
已经有57人回复
【求助】色谱A投稿
已经有3人回复
【求助】(石墨烯)紧束缚模型下的几个问题
已经有3人回复
【求助】一题一金币
已经有8人回复
【求助】Arcgis与GS+在地统计分析中的区别
已经有39人回复
» 抢金币啦!回帖就可以得到:
查看全部散金贴
湖南大学化工院招收2026年入学博士
+
1
/190
2026年江西师范大学药学院陈芬儿院士课题组招收智能药学博士生
+
2
/126
2026年江西师范大学人文地理与低碳经济方向拟招收1名“申请入学制”学术博士研究生
+
2
/124
中科大 冯伟 课题组招生-机器人界面控制、电学机械界面设计
+
1
/84
2026年江西师范大学流动化学与工程研究院陈芬儿院士团队化学博士招生
+
3
/83
电子科技大学(深圳)高等研究院陈明豪课题组招收推免硕士研究生
+
1
/82
温州医科大学吴平课题组招聘博士后,年薪≥50W/年
+
1
/80
山东征女友,坐标济南
+
1
/73
生命学部的面地青本子开始分发了么?
+
1
/64
双一流高校-南京林业大学-化学工程学院-国家海外优青团队招2026级博士(5月15号截止)
+
1
/56
西北工业大学物理科学与技术学院-诚聘【AI与光学交叉方向】博士后
+
1
/31
坐标北京,83年男征女友
+
1
/28
坐标南京
+
1
/22
北京大学新材料学院招聘柔性电子方向博士后
+
2
/8
南京大学医学院附属鼓楼医院院士团队招聘博士后、科研助理
+
1
/6
浙江大学长江学者特聘教授黄飞鹤招聘博士后(年薪35万+)
+
1
/5
招收2026年秋季入学博士生1名(河北工业大学/北京科技大学联合 增材制造/生物材料)
+
1
/5
英国伦敦布鲁内尔大学高薪招聘欧盟玛丽居里全奖博士
+
1
/5
[b] 南方科技大学 张作泰国家杰青 2026级全日制博士招生(固废资源化+深度减碳)[/b]
+
1
/3
Postdoc Position in SLAM & Spatial AI – University of Edinburgh
+
1
/2
1楼
2010-11-24 16:15:40
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
zyj8119
木虫
(著名写手)
应助: 65
(初中生)
贵宾: 0.003
金币: 915.1
帖子: 2936
在线: 1329.4小时
虫号: 664177
引用回帖:
Originally posted by
imyourkobe
at 2010-11-24 21:23:31:
你 上面
#define beta b*P/(R*T)
#define q a/(b*R*T)
用到了T 和P,这说明T和P是常量。
与主程序中要区分开。
那该怎么改呢?
回复此楼
高级回复
9楼
2010-11-24 21:29:57
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
查看全部 13 个回答
imyourkobe
铁杆木虫
(著名写手)
应助: 26
(小学生)
贵宾: 0.126
金币: 7477
帖子: 1688
在线: 459.7小时
虫号: 183723
★
ajian04(金币+1):谢谢参与交流~ 2010-11-24 20:48:45
呵呵,错误很多的啊。你到底要问i什么,是需要帮助修改错误?
赞
一下
(1人)
回复此楼
2楼
2010-11-24 18:35:29
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
imyourkobe
铁杆木虫
(著名写手)
应助: 26
(小学生)
贵宾: 0.126
金币: 7477
帖子: 1688
在线: 459.7小时
虫号: 183723
zyj8119(金币+5):那么程序应该怎么改呢?高手? 2010-11-24 19:09:15
ajian04:谢谢参与交流~ 2010-11-24 20:48:53
其中:z0(k),z1(k)要定义为数组
double z0[100],z1[100];
double z0(k)=M(s,1)=0;初始化要重新改写。
例如:
for(int i=0;i
z0
=0
赞
一下
(1人)
回复此楼
3楼
2010-11-24 18:39:23
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
zyj8119
木虫
(著名写手)
应助: 65
(初中生)
贵宾: 0.003
金币: 915.1
帖子: 2936
在线: 1329.4小时
虫号: 664177
ajian04:谢谢参与交流~ 2010-11-24 20:48:59
引用回帖:
Originally posted by
imyourkobe
at 2010-11-24 18:35:29:
呵呵,错误很多的啊。你到底要问i什么,是需要帮助修改错误?
希望您帮忙修改下,我是新手。。。。。
赞
一下
(1人)
回复此楼
4楼
2010-11-24 19:10:39
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
查看全部 13 个回答
如果回帖内容含有宣传信息,请如实选中。否则帐号将被全论坛禁言
普通表情
龙
兔
虎
猫
高级回复
(可上传附件)
百度网盘
|
360云盘
|
千易网盘
|
华为网盘
在新窗口页面中打开自己喜欢的网盘网站,将文件上传后,然后将下载链接复制到帖子内容中就可以了。
信息提示
关闭
请填处理意见
关闭
确定