原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
合成氨计算程序
#includestdio.h
#includemath.h
double react(double n43,double n31,double n32,double n33,double f3,double k)
{ double n41,n42,f4,f;
n41=(2*n31+n33-n43)/2;
n42=(2*n32+3*n33-3*n43)/2;
f4=f3+n33-n43;
f=pow(n43,2)/(n41*pow(n42,3))-k/pow(f4,2);
return f;
}
double flash(double x[5],double k[5],double a,double er)
{
int i=0;double f,df,b;
L90:f=0;df=0;
for(i=0;i5;i++)
{f+=x[i]*(1-k[i])/(a+k[i]);
df=df-x[i]*(1-k[i])/pow((a+k[i]),2);
}
b=a-f/df;
if(fabs(b-a)/aer)
{ a=b;goto L90;}
printf(b=%8.4f\t,b);
return b;
}
int main()
{
int i=0,j=0,no=0;
double n[10][5],x1[5]={0.24,0.743,0,0.006,0.011},x4[5],x5[5],x7[5]={0.25,0.75,0,0,0},x10[5]={0.25,0.75,0,0,0},f1=0,f2=0,fd,
k1[5]={105,90,0.06,100,33},k2[5]={2400,1750,0.28,1400,500},x[10][5],w[10][5],f[10],g[10],m[5]={28,2,17,40,16},xm[10],
kp=0.35,a1=0,a2=0,ap=0.02,er=1e-3,a11=0.1,a21=0.1,n43=80,z=10,
f7=100.0,f10=100.0;
f[0]=100;
printf(Read initial values:f7=%6.2f,f10=%6.2f,a21=%10.4f,a11=%10.4f,n43=%10.4f,z=%5.2f\n,f7,f10,a21,a11,n43,z);
for(i=0;i5;i++)
printf(Read initial values:x7[%d]=%7.4f,x10[%d]=%7.4f,i,i,x7[i],x10[i]);
for(i=0;i5;i++)
{x[0][i]=x1[i];
printf(%f,x[0][i]);}
L18:f[1]=f[0]+f7;
f[2]=f[1]+f10;
for(i=0;i5;i++)
{
n[0][i]=f[0]*x[0][i];
n[6][i]=f7*x7[i];
n[9][i]=f10*x10[i];
n[1][i]=n[0][i]+n[6][i];
n[2][i]=n[1][i]+n[9][i];
x[1][i]=n[1][i]/f[1];
x[2][i]=n[2][i]/f[2];
}
L6:f1=react(n43,n[2][0],n[2][1],n[2][2],f[2],kp);
n43=n43+er;
f2=react(n43,n[2][0],n[2][1],n[2][2],f[2],kp);
fd=(f2-f1)/er;
n43=n43-er;
n[3][2]=n43-f1/fd;
printf(n[3][2]=%7.4f\n,n[3][2]);
if(fabs(n[3][2]-n43)/n43er)
goto L5;
else{n43=n[3][2];goto L6;}
L5:n[3][3]=n[2][3];
n[3][4]=n[2][4];
n[3][0]=(2*n[2][0]+n[2][2]-n[3][2])/2;
n[3][1]=(2*n[2][1]+3*n[2][2]-3*n[3][2])/2;
f[3]=n[3][0]+n[3][1]+n[3][2]+n[3][3]+n[3][4];
for(i=0;i5;i++)
{x[3][i]=n[3][i]/f[3];
x4[i]=x[3][i];
printf(x4[%d]=%6.4f\t,i,x4[i]);}
printf(f[3]=%8.4f\t,f[3]);
a1=flash(x4,k1,a11,er);
printf(a1=%10
文档评论（0）