版块导航
正在加载中...
客户端APP下载
论文辅导
申博辅导
登录
注册
帖子
帖子
用户
本版
应《网络安全法》要求,自2017年10月1日起,未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用,请尽快对帐号进行手机号验证,感谢您的理解与支持!
24小时热门版块排行榜
>
论坛更新日志
(3337)
>
虫友互识
(311)
>
文献求助
(286)
>
休闲灌水
(278)
>
导师招生
(195)
>
考博
(80)
>
硕博家园
(71)
>
招聘信息布告栏
(68)
>
论文投稿
(59)
>
基金申请
(52)
>
博后之家
(51)
>
考研
(42)
>
教师之家
(40)
>
公派出国
(30)
>
绿色求助(高悬赏)
(22)
>
外文书籍求助
(22)
小木虫论坛-学术科研互动平台
»
计算模拟区
»
程序语言
»
C/C++
»
【求助】请高手指出一下此程序有多少错误?
5
1/1
返回列表
查看: 673 | 回复: 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);
}
回复此楼
» 猜你喜欢
假如你的研究生提出不合理要求
已经有8人回复
萌生出自己或许不适合搞科研的想法,现在跑or等等看?
已经有4人回复
Materials Today Chemistry审稿周期
已经有4人回复
参与限项
已经有3人回复
实验室接单子
已经有4人回复
全日制(定向)博士
已经有4人回复
对氯苯硼酸纯化
已经有3人回复
求助:我三月中下旬出站,青基依托单位怎么办?
已经有12人回复
所感
已经有4人回复
要不要辞职读博?
已经有7人回复
高级回复
» 本主题相关价值贴推荐,对您同样有帮助:
【求助】write()输出疑问?
已经有8人回复
【求助】增塑剂除色
已经有2人回复
【求助】急求高手帮忙!只有cif如何晶体画图(绝对没造假!)
已经有10人回复
【整理】《提问的智慧》——献给那些不会提问的人
已经有78人回复
【求助】关于用液相进行甲醇测试的问题
已经有10人回复
【求助/交流】一封老外教授给我的回信~(8.4更新)!
已经有57人回复
【求助】色谱A投稿
已经有3人回复
【求助】(石墨烯)紧束缚模型下的几个问题
已经有3人回复
【求助】一题一金币
已经有8人回复
【求助】Arcgis与GS+在地统计分析中的区别
已经有39人回复
» 抢金币啦!回帖就可以得到:
查看全部散金贴
散金
+
5
/650
科瑞赛生物内皮细胞培养基试用装限时大放送,助力你的实验高效进阶!
+
1
/83
供应德国EXAKT艾卡特半导体导热散热材料三辊研磨机50 PLUS
+
1
/82
上海大学昝鹏教授、军事医学研究院伯晓晨研究员/倪铭副研究员 课题组招聘博士生
+
1
/77
上海师范大学生物医用材料方向招收2026级博士研究生
+
1
/77
锌离子混合电容器
+
1
/75
上海大学昝鹏教授、军事医学研究院伯晓晨研究员/倪铭副研究员 课题组招聘博士生
+
1
/73
浙江师范大学国家杰青杨启华教授团队招收2026年博士研究生
+
1
/68
华中科技大学郭睿副教授建筑与能源科学课题组诚招博士研究生(2026.1.19截止)
+
1
/50
江西理工大学 稀土学院(发光材料与器件研究所) 招收2026届 材料类博士研究生 2名
+
2
/44
江西理工大学稀土学院稀土发光材料研究所招收2026届材料科学与工程专业博士研究生2名
+
2
/40
上海科技大学物质科学与技术学院|王平鸾课题组联合招聘博士后
+
1
/29
厦门大学航空航天学院智能制造课题组招2026年申请审核制博士生1-2名
+
1
/27
中科院理化技术研究所张飞龙研究员/王树涛研究员团队招生(博士/硕士)
+
1
/27
博士后招聘(高薪40万+)
+
1
/20
科研党/导师看过来,强推这个自带“引文验真”的国产工具,改作业效率翻倍
+
1
/15
【博士后/科研助理招聘-北京理工大学-集成电路与电子学院-国家杰青团队】
+
1
/5
中国矿业大学(徐州)蒋江民课题组招收2026年秋季新能源储能材料与器件方向博士生
+
1
/4
中国科学院大学-杨晗课题组-诚聘-博士后、副研究员
+
1
/3
中国科学院大学-杨晗课题组-诚聘-博士后、副研究员
+
1
/1
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
金币: 7471
帖子: 1688
在线: 458.9小时
虫号: 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
金币: 7471
帖子: 1688
在线: 458.9小时
虫号: 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云盘
|
千易网盘
|
华为网盘
在新窗口页面中打开自己喜欢的网盘网站,将文件上传后,然后将下载链接复制到帖子内容中就可以了。
信息提示
关闭
请填处理意见
关闭
确定