omiga欧美茄官网 2.0 破解版，有的能提供一下吗？Q,784936810,谢谢

来源：蜘蛛抓取(WebSpider) 时间：2015-03-31 08:40 标签： omiga欧米茄女款手表
&& 查看话题
【求助】有很多错误的程序，找不出来了
c& &&&此程序是运用peng-robinson方程计算气体的逸度系数faita(1),faita(2)，是利用牛顿迭代法求解
c& &&&然后通过origin拟合得到这些数据N1,N2,m1,m2,K1,K2
c& &&&求解得到f1,f2,x1的过程,首先利用牛顿迭代法求解x1, 然后利用表达式求解f1,f2
& && &PROGRAM MAIN
& & & & LOGICAL converged
& & & & REAL EXTERNAL FUNCTION MS MSD&&G dG
& & & & integer::S N SQ
& & & & parameter(EPS=1.0E-06,R=83.14,epsilon=1-2**(0.5),sigma=1+2**(0.5)
& & & & dimension:: faita(7),ns(4),ma(4),K(4),A(1000000),M(7),omiga(7),
& &&&*Tc(7),Pc(7),Z(1000000),b(7),beta(7),alpha(7),q(7),as(7),Z1(7)
& && &real:: alpha,beta,Tr,omiga,M,rho(7),q,a,b,SV,I(7)
& & & & DATA&&(ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,
& &&&*4.0/);K=(/5.0,6.0,7.0,8.0/);M=(/16.043,30.070,44.097,
& &&&*58.123,72.150, 44.01, 28.01/);omiga=(/0.012,0.100,0.152,0.2,
& &&&*0.252,0.224,0.048/);Tc=(/190.6,305.3,369.8,425.1,469.7,304.2,
& &&&*132.9/);Pc=(/45.99,48.72,42.48,37.96,33.70,73.83,34.99/)
& & & & READ(*,*)P1
& & & & P=P1/100
& && &SV=(1-(T/Tc(N))**(0.5))**2
& && &alpha(N)=(1+(0.26*omiga(N)-0.26992*omiga(N)**2)*SV)
& && &as(N)=((0.45724*R**2*Tc(N)**2)/Pc(N))*alpha(N)
& && &b(N)=0.07779*R*Tc(N)/Pc(N)
& && &beta(N)=b(N)*P/(R*T)
& && &q(N)=as(N)*alpha(N)*P/(R**2*T**2)
& && &MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
& &&&*(q*beta(N)-beta(N)**2-beta(N)**3)
& & & & MSD(Z)=3*Z**2-2*(1-beta(N))*Z+(q(N)-2*beta(N)-3*beta(N)**2)
& && &IF(N.LT.7.AND.N.GT.1)THEN
& & & & DO WHILE(.not.converged.and.ABS(MSD(Z)).GT.0)
& & & & Z=Z-MS(Z)/MSD(Z)
& & & & converged=(abs(MS(Z)/MSD(Z)).LT.EPS)
& && &END DO
& & & & ELSE
& & & & END IF
& & & & Z1(N)=Z
& && &DO 4 N=1,7
& & & & I(N)=(1/(sigma-epsilon))*log((Z1(N)+sigma*beta(N))/
& &&&*(Z1(N)+epsilon*beta(N)))
& & & & rho(N)=(P/(Z1(N)*R*T))*M(N)*EPS
& && &faita(N)=exp(Z1(N)-1-log(Z1(N)-beta)-q(N)*I(N))
4& & CONTINUE& &&&
& & & & f1=P*y1*faita(1)/X
& && &f2=P*(1-y1)*faita(2)/(1-X)
& & & & G(X)=ns(1)*LOG(1+K(1)*f1**m(1))+ns(2)*LOG(1+K(2)*f1**ma(2))-
& &&&*ns(3)*LOG(1+K(3)*f2**m(3))+ns(4)*LOG(1+K(4)*f2**ma(4))& & & && &&&
& && &dG(X)=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+K(1)*f1**ma(1))+
& &&&*(ma(2)*ns(2)*K2*f1*(ma(2)-1))/(1+K(2)*f1**m(2))-
& &&&*(ma(3)*ns(3)*K(3)*f2*(ma(3)-1))/(1+K(3)*f2**ma(3))+
& &&&*(ma(4)*ns(4)*K4*f2*(ma(4)-1))/(1+K(4)*f2**ma(4))
& && &DO WHILE(.not.converged.and.ABS(dG(X)).GT.0)
& && &X=X-G(X)/dG(X)
& && &S=S+1
& && &converged=(abs(G(X)/dG(X)).LT.EPS)
& && &END DO
& & & & x1=X
20& & format(1x,3F9.4)
& && &write(20,*)f1,f2,x1
& & & & end
& && &--------------------Configuration: bb - Win32 Debug--------------------
Compiling Fortran...
D:\vc++\MSDev98\MyProjects\22\bb.for
D:\vc++\MSDev98\MyProjects\22\bb.for(14) : Error: Syntax error, found ';' when expecting one of: , )
DATA&&(ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,
-------------------------------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(27) : Error: This name has not been declared as an array.& &
& && &MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
D:\vc++\MSDev98\MyProjects\22\bb.for(27) : Error: The shapes of the array expressions do not conform.
& && &MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
-----------------------------------------------------------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(27) : Error: The shapes of the array expressions do not conform.& &
& && &MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
D:\vc++\MSDev98\MyProjects\22\bb.for(29) : Error: This name has not been declared as an array.& &
MSD(Z)=3*Z**2-2*(1-beta(N))*Z+(q(N)-2*beta(N)-3*beta(N)**2)
D:\vc++\MSDev98\MyProjects\22\bb.for(29) : Error: The shapes of the array expressions do not conform.& &
MSD(Z)=3*Z**2-2*(1-beta(N))*Z+(q(N)-2*beta(N)-3*beta(N)**2)
D:\vc++\MSDev98\MyProjects\22\bb.for(37) : Error: The shapes of the array expressions do not conform.& &
D:\vc++\MSDev98\MyProjects\22\bb.for(43) : Error: The shapes of the array expressions do not conform.& &
& && &faita(N)=exp(Z1(N)-1-log(Z1(N)-beta)-q(N)*I(N))
D:\vc++\MSDev98\MyProjects\22\bb.for(49) : Error: This name has not been declared as an array.& &
G(X)=ns(1)*LOG(1+K(1)*f1**m(1))+ns(2)*LOG(1+K(2)*f1**ma(2))-
D:\vc++\MSDev98\MyProjects\22\bb.for(51) : Error: This name has not been declared as an array.& &
& && &dG(X)=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+K(1)*f1**ma(1))+
Error executing df.exe.
bb.obj - 10 error(s), 0 warning(s)
改成这样，还有7个错误& &
c& &&&此程序是运用peng-robinson方程计算气体的逸度系数faita(1),faita(2)，是利用牛顿迭代法求解
c& &&&然后通过origin拟合得到这些数据N1,N2,m1,m2,K1,K2
c& &&&求解得到f1,f2,x1的过程,首先利用牛顿迭代法求解x1, 然后利用表达式求解f1,f2
& && &PROGRAM MAIN
& & & & LOGICAL converged
& & & & REAL EXTERNAL FUNCTION MS MSD&&G dG
& & & & integer::S N SQ
& & & & parameter(EPS=1.0E-06,R=83.14,epsilon=1-2**(0.5),sigma=1+2**(0.5)
& & & & dimension:: faita(7),ns(4),ma(4),K(4),M(7),omiga(7),
& &&&*Tc(7),Pc(7),b(7),beta(7),alpha(7),q(7),as(7),Z1(7)
& && &real:: alpha,beta,Tr,omiga,M,rho(7),q,a,b,SV,I(7)
& & & & DATA&&(ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,4.0/);
& &&&*K=(/5.0,6.0,7.0,8.0/);M=(/16.043,30.070,44.097,
& &&&*58.123,72.150, 44.01, 28.01/);omiga=(/0.012,0.100,0.152,0.2,
& &&&*0.252,0.224,0.048/);Tc=(/190.6,305.3,369.8,425.1,469.7,304.2,
& &&&*132.9/);Pc=(/45.99,48.72,42.48,37.96,33.70,73.83,34.99/)
& & & & READ(*,*)P1
& & & & P=P1/100
& && &SV=(1-(T/Tc(N))**(0.5))**2
& && &alpha(N)=(1+(0.26*omiga(N)-0.26992*omiga(N)**2)*SV)
& && &as(N)=((0.45724*R**2*Tc(N)**2)/Pc(N))*alpha(N)
& && &b(N)=0.07779*R*Tc(N)/Pc(N)
& && &beta(N)=b(N)*P/(R*T)
& && &q(N)=as(N)*alpha(N)*P/(R**2*T**2)
& && &MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
& &&&*(q*beta(N)-beta(N)**2-beta(N)**3)
& & & & MSD(Z)=3*Z**2-2*(1-beta(N))*Z+(q(N)-2*beta(N)-3*beta(N)**2)
& && &IF(N.LT.7.AND.N.GT.1)THEN
& & & & DO WHILE(.not.converged.and.ABS(MSD(Z)).GT.0)
& & & & Z=Z-MS(Z)/MSD(Z)
& & & & converged=(abs(MS(Z)/MSD(Z)).LT.EPS)
& && &END DO
& & & & ELSE
& & & & END IF
& & & & Z1(N)=Z
& && &DO 4 N=1,7
& & & & I(N)=(1/(sigma-epsilon))*log((Z1(N)+sigma*beta(N))/
& &&&*(Z1(N)+epsilon*beta(N)))
& & & & rho(N)=(P/(Z1(N)*R*T))*M(N)*EPS
& && &faita(N)=exp(Z1(N)-1-log(Z1(N)-beta)-q(N)*I(N))
4& & CONTINUE& &&&
& & & & f1=P*y1*faita(1)/X
& && &f2=P*(1-y1)*faita(2)/(1-X)
& & & & G(X)=ns(1)*LOG(1+K(1)*f1**m(1))+ns(2)*LOG(1+K(2)*f1**ma(2))-
& &&&*ns(3)*LOG(1+K(3)*f2**m(3))+ns(4)*LOG(1+K(4)*f2**ma(4))& & & && &&&
& && &dG(X)=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+K(1)*f1**ma(1))+
& &&&*(ma(2)*ns(2)*K2*f1*(ma(2)-1))/(1+K(2)*f1**m(2))-
& &&&*(ma(3)*ns(3)*K(3)*f2*(ma(3)-1))/(1+K(3)*f2**ma(3))+
& &&&*(ma(4)*ns(4)*K4*f2*(ma(4)-1))/(1+K(4)*f2**ma(4))
& && &DO WHILE(.not.converged.and.ABS(dG(X)).GT.0)
& && &X=X-G(X)/dG(X)
& && &S=S+1
& && &converged=(abs(G(X)/dG(X)).LT.EPS)
& && &END DO
& & & & x1=X
20& & format(1x,3F9.4)
& && &write(20,*)f1,f2,x1
& & & & end
--------------------Configuration: bb - Win32 Debug--------------------
Compiling Fortran...
D:\vc++\MSDev98\MyProjects\22\bb.for
D:\vc++\MSDev98\MyProjects\22\bb.for(14) : Error: Syntax error, found ';' when expecting one of: , )
DATA&&(ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,4.0/);
-------------------------------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(27) : Error: This name has not been declared as an array.& &
& && &MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
D:\vc++\MSDev98\MyProjects\22\bb.for(27) : Error: The shapes of the array expressions do not conform.& &
& && &MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
D:\vc++\MSDev98\MyProjects\22\bb.for(29) : Error: This name has not been declared as an array.& &
MSD(Z)=3*Z**2-2*(1-beta(N))*Z+(q(N)-2*beta(N)-3*beta(N)**2)
D:\vc++\MSDev98\MyProjects\22\bb.for(43) : Error: The shapes of the array expressions do not conform.& &
& && &faita(N)=exp(Z1(N)-1-log(Z1(N)-beta)-q(N)*I(N))
D:\vc++\MSDev98\MyProjects\22\bb.for(49) : Error: This name has not been declared as an array.& &
G(X)=ns(1)*LOG(1+K(1)*f1**m(1))+ns(2)*LOG(1+K(2)*f1**ma(2))-
D:\vc++\MSDev98\MyProjects\22\bb.for(51) : Error: This name has not been declared as an array.& &
& && &dG(X)=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+K(1)*f1**ma(1))+
Error executing df.exe. 看不懂你要做什么，不过有几点建议：
你的程序似乎是 Fortran77 的固定格式，而且里面的循环什么的也是按 Fortran 77 样式写的。但里面有些写法似乎又是 Fortran 90 的样式，这样很不好，建议按一种风格来写你的程序，要么就是纯 Fortran 77 的，要么就是 Fortran90 的。
你的程序中有很多变量没有初值，比如
SV=(1-(T/Tc(N))**(0.5))**2
MS(Z)=Z**3-(1-beta(N))*Z**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-(q*beta(N)-beta(N)**2-beta(N)**3)
中的 Z。如果你这里的 MS 是一个函数（就像你前面申明的那样），那么最好用别的符号来代替 Z，否则会引起歧义。
这个程序中最大的问题也就在这些个函数的定义上面了…… Originally posted by snoopyzhao at
看不懂你要做什么，不过有几点建议：
你的程序似乎是 Fortran77 的固定格式，而且里面的循环什么的也是按 Fortran 77 样式写的。但里面有些写法似乎又是 Fortran 90 的样式，这样很不好，建议按一种风格来写你 ... 我想定义一个函数，但是错误里面好像非要我搞成数组啊，MS(Z)是一个Z的函数，MSD(Z)是MS(Z)的导数，这样可以利用牛顿迭代法求解，下面的也是，G(X)是函数本身，&&dG(X)是G(X)的导数。 $& &&&此程序是运用peng-robinson方程计算气体的逸度系数faita(1),faita(2)，是利用牛顿迭代法求解
$& &&&然后通过origin拟合得到这些数据N1,N2,m1,m2,K1,K2
$& &&&求解得到f1,f2,x1的过程,首先利用牛顿迭代法求解IAST理论得到x1, 然后利用表达式求解f1,f2
& && &PROGRAM MAIN
& & & & LOGICAL converged
& & & & REAL EXTERNAL FUNCTION MS MSD&&G dG
& & & & integer::S N SQ
& & & & parameter(EPS=1.0E-06,R=83.14,epsilon=1-2**(0.5),sigma=1+2**(0.5)
& & & & dimension:: faita(7),ns(4),ma(4),K(4),M(7),omiga(7),
& &&&*Tc(7),Pc(7),b(7),beta(7),alpha(7),q(7),as(7),Z1(7)
& && &real:: alpha,beta,Tr,omiga,M,rho(7),q,a,b,SV,I(7),ZM
& & & & DATA&&(ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,4.0/);
& &&&*K=(/5.0,6.0,7.0,8.0/);M=(/16.043,30.070,44.097,
& &&&*58.123,72.150, 44.01, 28.01/);omiga=(/0.012,0.100,0.152,0.2,
& &&&*0.252,0.224,0.048/);Tc=(/190.6,305.3,369.8,425.1,469.7,304.2,
& &&&*132.9/);Pc=(/45.99,48.72,42.48,37.96,33.70,73.83,34.99/)
& & & & READ(*,*)P1
$& &&&计算peng-robinson方程中需要的一些参量& & & &
& & & & DO 3 N=1,7
& & & & P=P1/100
& && &SV=(1-(T/Tc(N))**(0.5))**2
& && &alpha(N)=(1+(0.26*omiga(N)-0.26992*omiga(N)**2)*SV)
& && &as(N)=((0.45724*R**2*Tc(N)**2)/Pc(N))*alpha(N)
& && &b(N)=0.07779*R*Tc(N)/Pc(N)
& && &beta(N)=b(N)*P/(R*T)
& && &q(N)=as(N)*alpha(N)*P/(R**2*T**2)
3& &&&CONTINUE& && &
$& &&&运用牛顿迭代法开始求解压缩因子的方程& & & &
& & & & MS(ZM)=ZM**3-(1-beta(N))*ZM**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
& &&&*(q*beta(N)-beta(N)**2-beta(N)**3)
& & & & MSD(ZM)=3*ZM**2-2*(1-beta(N))*ZM+(q(N)-2*beta(N)-3*beta(N)**2)
& && &IF(N.LT.7.AND.N.GT.1)THEN
& & & & DO WHILE(.not.converged.and.ABS(MSD(ZM)).GT.0)
& & & & ZM=ZM-MS(ZM)/MSD(ZM)
& & & & converged=(abs(MS(ZM)/MSD(ZM)).LT.EPS)
& && &END DO
& & & & ELSE
& & & & END IF
& & & & Z1(N)=ZM
$& &&&求解密度和逸度& &
& && &DO 4 N=1,7
& & & & I(N)=(1/(sigma-epsilon))*log((Z1(N)+sigma*beta(N))/
& &&&*(Z1(N)+epsilon*beta(N)))
& & & & rho(N)=(P/(Z1(N)*R*T))*M(N)*EPS
& && &faita(N)=exp(Z1(N)-1-log(Z1(N)-beta)-q(N)*I(N))
4& & CONTINUE& &&&
$& &&&运用IAST理论以及上面得到的逸度系数，开始求解IAST的方程
& & & & f1=P*y1*faita(1)/X
& && &f2=P*(1-y1)*faita(2)/(1-X)
& & & & G(X)=ns(1)*LOG(1+K(1)*f1**m(1))+ns(2)*LOG(1+K(2)*f1**ma(2))-
& &&&*ns(3)*LOG(1+K(3)*f2**m(3))+ns(4)*LOG(1+K(4)*f2**ma(4))& & & && &&&
& && &dG(X)=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+K(1)*f1**ma(1))+
& &&&*(ma(2)*ns(2)*K2*f1*(ma(2)-1))/(1+K(2)*f1**m(2))-
& &&&*(ma(3)*ns(3)*K(3)*f2*(ma(3)-1))/(1+K(3)*f2**ma(3))+
& &&&*(ma(4)*ns(4)*K4*f2*(ma(4)-1))/(1+K(4)*f2**ma(4))
& && &DO WHILE(.not.converged.and.ABS(dG(X)).GT.0)
& && &X=X-G(X)/dG(X)
& && &S=S+1
& && &converged=(abs(G(X)/dG(X)).LT.EPS)
& && &END DO
& & & & x1=X
20& & format(1x,3F9.4)
& && &write(20,*)f1,f2,x1
& & & & end 下面这两个方程，是不能定义成这个样子的。
& && &&&MS(ZM)=ZM**3-(1-beta(N))*ZM**2+(q(N)-2*beta(N)-3*beta(N)**2)*Z-
& &&&*(q*beta(N)-beta(N)**2-beta(N)**3)
& && &&&MSD(ZM)=3*ZM**2-2*(1-beta(N))*ZM+(q(N)-2*beta(N)-3*beta(N)**2)
从你的代码上看，这两个方程中的 beta(N), q(N) 啥的都是变量，呵呵……
你这么搞搞看：
MS(ZM,betaN,qN)=ZM**3-(1-betaN)*ZM**2+(qN-2*betaN-3*betaN**2)*ZM-(qN*betaN-betaN**2-betaN**3)
然后调用的时候，记得给 ZM, betaN 和 qN 赋值。其它的函数应该类似。 在贴代码时，最好用 Discuz! 代码那块儿显示的 #。这样别人能够很方便地看或拷贝了……
另外，要自己为所有的变量赋初值，比如 MS(ZM) 中的 ZM 就没有赋初值。另外，MS(ZM) 中用到 beta(N)，这里的 N 又是哪一个值呢？ $& &&&此程序是运用peng-robinson方程计算气体的逸度系数faita(1),faita(2)，是利用牛顿迭代法求解
$& &&&然后通过origin拟合得到这些数据N1,N2,m1,m2,K1,K2
$& &&&求解得到f1,f2,x1的过程,首先利用牛顿迭代法求解IAST理论得到x1, 然后利用表达式求解f1,f2
& && &PROGRAM MAIN
& & & & LOGICAL converged
& & & & integer S N SQ
& & & & parameter(EPS=1.0E-06,R=83.14,epsilon=1-2**(0.5),sigma=1+2**(0.5)
& & & & dimension faita(7),ns(4),ma(4),K(4),M(7),omiga(7),
& &&&*Tc(7),Pc(7),b(7),beta(7),alpha(7),q(7),as(7),Z1(7)
& && &real alpha,beta,Tr,omiga,M,rho(7),q,a,b,SV,I(7),ZM,SSA,SSB,SSC,
& & & & DATA&&M=(/16.043,30.070,44.097, 58.123,72.150, 44.01, 28.01/);
& &&&*& & & & omiga=(/0.012,0.100,0.152,0.2,
& &&&*0.252,0.224,0.048/);Tc=(/190.6,305.3,369.8,425.1,469.7,304.2,
& &&&*132.9/);Pc=(/45.99,48.72,42.48,37.96,33.70,73.83,34.99/)
& & & & READ(*,*)P1
$& &&&计算peng-robinson方程中需要的一些参量& & & &
& & & & DO 3 N=1,7
& & & & P=P1/100
& && &SV=(1-(T/Tc(N))**(0.5))**2
& && &alpha(N)=(1+(0.26*omiga(N)-0.26992*omiga(N)**2)*SV)
& && &as(N)=((0.45724*R**2*Tc(N)**2)/Pc(N))*alpha(N)
& && &b(N)=0.07779*R*Tc(N)/Pc(N)
& && &beta(N)=b(N)*P/(R*T)
& && &q(N)=as(N)*alpha(N)*P/(R**2*T**2)
3& &&&CONTINUE& && &
$& &&&运用牛顿迭代法开始求解压缩因子的方程& & & &
& && &IF(N.LT.7.AND.N.GT.1)THEN
& & & & ZM=1
& & & & DO WHILE(.not.converged.and.ABS(MSD(ZM)).GT.0)
& & & & SSA=MS(ZM,beta(N),q(N))
& & & & SSB=MSD(ZM,beta(N),q(N))
& & & & ZM=ZM-SSA/SSB
& & & & converged=(abs(SSA)/(SSB))).LT.EPS)
& && &END DO
& & & & ELSE
& & & & END IF
& & & & do 4 N=1,7
& & & & Z1(N)=ZM
4& &&&CONTINUE
$& &&&求解密度和逸度& &
& && &DO 5 N=1,7
& & & & I(N)=(1/(sigma-epsilon))*log((Z1(N)+sigma*beta(N))/
& &&&*(Z1(N)+epsilon*beta(N)))
& & & & rho(N)=(P/(Z1(N)*R*T))*M(N)*EPS
& && &faita(N)=SVF(Z1(N),beta,q(N),I(N))
5& & CONTINUE& &&&
$& &&&运用IAST理论以及上面得到的逸度系数，开始求解IAST的方程
& & & & S=0
& & & & X=1
& & & & SSC=iast(X,P,y1,faita(1),faita(2))
& & & & SSD=diast(X,P,y1,faita(1),faita(2))
& && &DO WHILE(.not.converged.and.ABS(SSD).GT.0)
& && &X=X-SSC/SSD
& && &S=S+1
& && &converged=(abs(SSC/SSD).LT.EPS)
& && &END DO
& & & & x1=X
20& & format(1x,3F9.4)
& && &write(20,*)f1,f2,x1
& & & & end
& & & & function MS(SA,SB,SC)
& && &MS(SA,SB,SC)=SA**3-(1-SB)*SA**2+(SC-2*SB-3*SB
& &&&***2)*SA-(SC*SB-SB**2-SB**3)
& & & & end
& & & & function MSD(SD,SE,SF)
& && &MSD(SD,SE,SF)=3*SD**2-2*(1-SE)*SD+(SF-2*SE-
& &&&*3*SE**2)
& & & & end
& & & & function & & & & iast(MADD,P,YY,SG1,SG2)
& && &dimension&&ns(4),ma(4),K(4)
& & & & DATA&&(ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,4.0/);
& &&&*K=(/5.0,6.0,7.0,8.0/)
& & & & f1=P*YY*SG1/MADD
& && &f2=P*(1-YY)*SG2/(1-MADD)
& & & & iast(MADD,P,YY,SG1,SG2)=ns(1)*LOG(1+K(1)*f1**ma(1))+ns(2)*LOG
& &&&*(1+K(2)*f1**ma(2))-ns(3)*LOG(1+K(3)*f2**m(3))+ns(4)*LOG(1+
& &&&*K(4)*f2**ma(4))
& & & & function&&diast(MAD,S,Y,SG3,SG4)
& && &dimension&&ns(4),ma(4),K(4)
& & & & DATA&&(ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,4.0/);
& &&&*K=(/5.0,6.0,7.0,8.0/)
& & & & f1=S*Y*SG3/MAD
& && &f2=S*(1-Y)*SG4/(1-MAD)
& & & & diast(MAD,S,Y,SG3,SG4)=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+
& &&&*K(1)*f1**ma(1))+(ma(2)*ns(2)*K2*f1*(ma(2)-1))/(1+K(2)*f1**m(2))-
& &&&*(ma(3)*ns(3)*K(3)*f2*(ma(3)-1))/(1+K(3)*f2**ma(3))+
& &&&*(ma(4)*ns(4)*K4*f2*(ma(4)-1))/(1+K(4)*f2**ma(4))
& & & & END
& & & & FUNCTION SVF(SG5,SG6,SG7,SG8)
& && &SVF(SG5,SG6,SG7,SG8)=exp(SG5-1-log(SG5-SG6)-SG7*SG8)
& & & & END
改成这样，还是有9个错误啊。。。 -----------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(76) : Error: This global name is invalid in this context.& &
& && &MSD(SD,SE,SF)=3*SD**2-2*(1-SE)*SD+(SF-2*SE-
D:\vc++\MSDev98\MyProjects\22\bb.for(75) : Warning: The return value of this FUNCTION has not been defined.& &
function MSD(SD,SE,SF)
-----------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(83) : Error: Syntax error, found ';' when expecting one of: , )
DATA&&(ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,4.0/);
-------------------------------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(87) : Error: This global name is invalid in this context.& &
iast(MADD,P,YY,SG1,SG2)=ns(1)*LOG(1+K(1)*f1**ma(1))+ns(2)*LOG
D:\vc++\MSDev98\MyProjects\22\bb.for(81) : Warning: The return value of this FUNCTION has not been defined.& &
function&&iast(MADD,P,YY,SG1,SG2)
------------------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(95) : Error: Syntax error, found ';' when expecting one of: , )
DATA&&(ns=(/2.0,3.0,4.0,5.0/);ma=(/1.0,2.0,3.0,4.0/);
-------------------------------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(99) : Error: This global name is invalid in this context.& &
diast(MAD,S,Y,SG3,SG4)=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+
D:\vc++\MSDev98\MyProjects\22\bb.for(93) : Warning: The return value of this FUNCTION has not been defined.& &
function&&diast(MAD,S,Y,SG3,SG4)
------------------^
D:\vc++\MSDev98\MyProjects\22\bb.for(106) : Error: This global name is invalid in this context.& &
& && &SVF(SG5,SG6,SG7,SG8)=exp(SG5-1-log(SG5-SG6)-SG7*SG8)
D:\vc++\MSDev98\MyProjects\22\bb.for(105) : Warning: The return value of this FUNCTION has not been defined.& &
FUNCTION SVF(SG5,SG6,SG7,SG8)
-----------------^
Error executing df.exe.
bb.obj - 9 error(s), 5 warning(s) beta(N)和q(N)本来就是计算得到的数组数值啊，怎么给初值呢，前面已经算了啊。 下面这个是可以编译通过的，但不保证准确性。建议你以后写程序时更认真一些，不要把 K(4) 写成 K4，把 ma(2) 写成 m(2) 等……
C& &&&此程序是运用peng-robinson方程计算气体的逸度系数faita(1),faita(2)，是利用牛顿迭代法求解
C& &&&然后通过origin拟合得到这些数据N1,N2,m1,m2,K1,K2
C& &&&求解得到f1,f2,x1的过程,首先利用牛顿迭代法求解IAST理论得到x1, 然后利用表达式求解f1,f2
& && &PROGRAM MAIN
& && &LOGICAL converged
& && &integer S N SQ
& && &parameter(EPS=1.0E-06,R=83.14,epsilon=1-2**(0.5),sigma=1+2**(0.5)
& && &dimension faita(7),ns(4),ma(4),K(4),M(7),omiga(7),
& &&&*Tc(7),Pc(7),b(7),beta(7),alpha(7),q(7),as(7),Z1(7)
& && &real alpha,beta,Tr,omiga,M,rho(7),q,a,b,SV,I(7),ZM,SSA,SSB,SSC,
& && &DATA&&M/16.043,30.070,44.097, 58.123,72.150, 44.01, 28.01/
& && &DATA omiga/0.012,0.100,0.152,0.2,0.252,0.224,0.048/
& && &DATA Tc/190.6,305.3,369.8,425.1,469.7,304.2,132.9/
& && &DATA Pc/45.99,48.72,42.48,37.96,33.70,73.83,34.99/
& && &READ(*,*)P1
C& &&&计算peng-robinson方程中需要的一些参量
& && &DO 3 N=1,7
& && &P=P1/100
& && &SV=(1-(T/Tc(N))**(0.5))**2
& && &alpha(N)=(1+(0.26*omiga(N)-0.26992*omiga(N)**2)*SV)
& && &as(N)=((0.45724*R**2*Tc(N)**2)/Pc(N))*alpha(N)
& && &b(N)=0.07779*R*Tc(N)/Pc(N)
& && &beta(N)=b(N)*P/(R*T)
& && &q(N)=as(N)*alpha(N)*P/(R**2*T**2)
3& &&&CONTINUE
C& &&&运用牛顿迭代法开始求解压缩因子的方程
& && &IF(N.LT.7.AND.N.GT.1)THEN
& && &&&ZM=1
& && &&&DO WHILE(.not.converged.and.ABS(MSD(ZM)).GT.0)
& && &&&SSA=MS(ZM,beta(N),q(N))
& && &&&SSB=MSD(ZM,beta(N),q(N))
& && &&&ZM=ZM-SSA/SSB
& && &&&converged=(abs(SSA)/(SSB)).LT.EPS
& && &&&END DO
& && &ELSE
& && &END IF
& && &do 4 N=1,7
& && &Z1(N)=ZM
4& &&&CONTINUE
C& &&&求解密度和逸度
& && &DO 5 N=1,7
& && &I(N)=(1/(sigma-epsilon))*log((Z1(N)+sigma*beta(N))/
& &&&*(Z1(N)+epsilon*beta(N)))
& && &rho(N)=(P/(Z1(N)*R*T))*M(N)*EPS
& && &faita(N)=SVF(Z1(N),beta,q(N),I(N))
5& &&&CONTINUE
C& &&&运用IAST理论以及上面得到的逸度系数，开始求解IAST的方程
& && &SSC=iast(X,P,y1,faita(1),faita(2))
& && &SSD=diast(X,P,y1,faita(1),faita(2))
& && &DO WHILE(.not.converged.and.ABS(SSD).GT.0)
& && &X=X-SSC/SSD
& && &S=S+1
& && &converged=(abs(SSC/SSD).LT.EPS)
& && &END DO
& && &x1=X
20& & format(1x,3F9.4)
& && &write(*,20)f1,f2,x1
& && &function MS(SA,SB,SC)
& && &MS=SA**3-(1-SB)*SA**2+(SC-2*SB-3*SB
& &&&***2)*SA-(SC*SB-SB**2-SB**3)
& && &function MSD(SD,SE,SF)
& && &MSD=3*SD**2-2*(1-SE)*SD+(SF-2*SE-
& &&&*3*SE**2)
& && &function& && && &iast(MADD,P,YY,SG1,SG2)
& && &dimension&&ns(4),ma(4),K(4)
& && &DATA&&ns/2.0,3.0,4.0,5.0/
& && &DATA ma/1.0,2.0,3.0,4.0/
& && &DATA K/5.0,6.0,7.0,8.0/
& && &f1=P*YY*SG1/MADD
& && &f2=P*(1-YY)*SG2/(1-MADD)
& && &iast=ns(1)*LOG(1+K(1)*f1**ma(1))+ns(2)*LOG
& &&&*(1+K(2)*f1**ma(2))-ns(3)*LOG(1+K(3)*f2**ma(3))+ns(4)*LOG(1+
& &&&*K(4)*f2**ma(4))
& && &function&&diast(MAD,S,Y,SG3,SG4)
& && &dimension&&ns(4),ma(4),K(4)
& && &DATA&&ns/2.0,3.0,4.0,5.0/
& && &data ma/1.0,2.0,3.0,4.0/
& && &data K/5.0,6.0,7.0,8.0/
& && &f1=S*Y*SG3/MAD
& && &f2=S*(1-Y)*SG4/(1-MAD)
& && &diast=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+
& &&&*K(1)*f1**ma(1))+(ma(2)*ns(2)*K2*f1*(ma(2)-1))/(1+K(2)*f1**ma(2))-
& &&&*(ma(3)*ns(3)*K(3)*f2*(ma(3)-1))/(1+K(3)*f2**ma(3))+
& &&&*(ma(4)*ns(4)*K(4)*f2*(ma(4)-1))/(1+K(4)*f2**ma(4))
& && &FUNCTION SVF(SG5,SG6,SG7,SG8)
& && &SVF=exp(SG5-1-log(SG5-SG6)-SG7*SG8)
& && &END Originally posted by zyj8119 at
beta(N)和q(N)本来就是计算得到的数组数值啊，怎么给初值呢，前面已经算了啊。 但你那个 MS(ZM) 函数中的 beta(N) 已不在循环中了啊，这个时候 N 就不知道该取什么值了…… Originally posted by snoopyzhao at
但你那个 MS(ZM) 函数中的 beta(N) 已不在循环中了啊，这个时候 N 就不知道该取什么值了…… 谢谢！！！！！！跪谢！！！！ C& &&&此程序是运用peng-robinson方程计算气体的逸度系数faita(1),faita(2)，是利用牛顿迭代法求解
C& &&&然后通过origin拟合得到这些数据N1,N2,m1,m2,K1,K2
C& &&&求解得到f1,f2,x1的过程,首先利用牛顿迭代法求解IAST理论得到x1, 然后利用表达式求解f1,f2
& && &PROGRAM MAIN
& && &LOGICAL converged
& && &integer S N SQ
& && &parameter(EPS=1.0E-06,R=83.14,epsilon=1-2**(0.5),sigma=1+2**(0.5)
& && &dimension faita(7),ns(4),ma(4),K(4),M(7),omiga(7),
& &&&*Tc(7),Pc(7),b(7),beta(7),alpha(7),q(7),as(7),Z1(7)
& && &real alpha,beta,Tr,omiga,M,rho(7),q,a,b,SV,I(7),ZM,SSA,SSB,SSC,
& &&&*SSD,T,f1,f2
& && &DATA M/16.043,30.070,44.097, 58.123,72.150, 44.01, 28.01/
& && &DATA omiga/0.012,0.100,0.152,0.2,0.252,0.224,0.048/
& && &DATA Tc/190.6,305.3,369.8,425.1,469.7,304.2,132.9/
& && &DATA Pc/45.99,48.72,42.48,37.96,33.70,73.83,34.99/
& && &READ(*,*)P1
C& &&&计算peng-robinson方程中需要的一些参量
& && &DO 3 N=1,7
& && &P=P1/100
& && &SV=(1-(T/Tc(N))**(0.5))**2
& && &alpha(N)=(1+(0.26*omiga(N)-0.26992*omiga(N)**2)*SV)
& && &as(N)=((0.45724*R**2*Tc(N)**2)/Pc(N))*alpha(N)
& && &b(N)=0.07779*R*Tc(N)/Pc(N)
& && &beta(N)=b(N)*P/(R*T)
& && &q(N)=as(N)*alpha(N)*P/(R**2*T**2)
3& &&&CONTINUE
C& &&&运用牛顿迭代法开始求解压缩因子的方程
& && &T=298.0
& && &IF(N.LT.7.AND.N.GT.1)THEN
& && &&&ZM=1
& && &&&DO WHILE(.not.converged.and.ABS(MSD(ZM,beta(N),q(N))).GT.0)
& && &&&SSA=MS(ZM,beta(N),q(N))
& && &&&SSB=MSD(ZM,beta(N),q(N))
& && &&&ZM=ZM-SSA/SSB
& && &&&converged=(abs(SSA)/(SSB)).LT.EPS
& && &&&END DO
& && &ELSE
& && &END IF
& && &do 4 N=1,7
& && &Z1(N)=ZM
4& &&&CONTINUE
C& &&&求解密度和逸度
& && &DO 5 N=1,7
& && &I(N)=(1/(sigma-epsilon))*log((Z1(N)+sigma*beta(N))/
& &&&*(Z1(N)+epsilon*beta(N)))
& && &rho(N)=(P/(Z1(N)*R*T))*M(N)*EPS
& && &faita(N)=SVF(Z1(N),beta,q(N),I(N))
5& &&&CONTINUE
C& &&&运用IAST理论以及上面得到的逸度系数，开始求解IAST的方程
& && &SSC=iast(X,P,y1,faita(1),faita(2))
& && &SSD=diast(X,P,y1,faita(1),faita(2))
& && &DO WHILE(.not.converged.and.ABS(SSD).GT.0)
& && &X=X-SSC/SSD
& && &S=S+1
& && &converged=(abs(SSC/SSD).LT.EPS)
& && &END DO
& && &x1=X
& && &f1=P*y1*faita(1)/x1
& && &f2=P*(1-y1)*faita(2)/(1-x1)
20& & format(1x,3F9.4)
& && &write(20,*)f1,f2,x1
& && &function MS(SA,SB,SC)
& && &MS=SA**3-(1-SB)*SA**2+(SC-2*SB-3*SB
& &&&***2)*SA-(SC*SB-SB**2-SB**3)
& && &function MSD(SD,SE,SF)
& && &MSD=3*SD**2-2*(1-SE)*SD+(SF-2*SE-
& &&&*3*SE**2)
& && &function& &iast(MADD,P,YY,SG1,SG2)
& && &dimension&&ns(4),ma(4),K(4)
& & & & REAL MADD
& && &DATA&&ns/2.0,3.0,4.0,5.0/
& && &DATA ma/1.0,2.0,3.0,4.0/
& && &DATA K/5.0,6.0,7.0,8.0/
& && &f1=P*YY*SG1/MADD
& && &f2=P*(1-YY)*SG2/(1-MADD)
& && &iast=ns(1)*LOG(1+K(1)*f1**ma(1))+ns(2)*LOG
& &&&*(1+K(2)*f1**ma(2))-ns(3)*LOG(1+K(3)*f2**ma(3))+ns(4)*LOG(1+
& &&&*K(4)*f2**ma(4))
& && &function&&diast(MAD,S,Y,SG3,SG4)
& && &dimension&&ns(4),ma(4),K(4)
& & & & REAL MAD
& && &DATA&&ns/2.0,3.0,4.0,5.0/
& && &data ma/1.0,2.0,3.0,4.0/
& && &data K/5.0,6.0,7.0,8.0/
& && &f1=S*Y*SG3/MAD
& && &f2=S*(1-Y)*SG4/(1-MAD)
& && &diast=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+
& &&&*K(1)*f1**ma(1))+(ma(2)*ns(2)*K(2)*f1*(ma(2)-1))/(1+K(2)*f1**ma(2))
& &&&*-(ma(3)*ns(3)*K(3)*f2*(ma(3)-1))/(1+K(3)*f2**ma(3))+
& &&&*(ma(4)*ns(4)*K(4)*f2*(ma(4)-1))/(1+K(4)*f2**ma(4))
& && &FUNCTION SVF(SG5,SG6,SG7,SG8)
& && &SVF=exp(SG5-1-log(SG5-SG6)-SG7*SG8)
此程序没有错误和warning了，但是就是运行不出结果，我是COMPAQ VISUAL FORTRAN. 运行不出结果是什么意思？没有输出，还是别的？
你的这一部分有问题，我在前面的程序中其实改过了……
20& & format(1x,3F9.4)
& && &write(20,*)f1,f2,x1
& && &IF(N.LT.7.AND.N.GT.1)THEN
& && &&&ZM=1
& && &&&DO WHILE(.not.converged.and.ABS(MSD(ZM,beta(N),q(N))).GT.0)
& && &&&SSA=MS(ZM,beta(N),q(N))
& && &&&SSB=MSD(ZM,beta(N),q(N))
& && &&&ZM=ZM-SSA/SSB
& && &&&converged=(abs(SSA)/(SSB)).LT.EPS
& && &&&END DO
& && &ELSE
& && &END IF
也有问题，一个是 N 在这里未定义，因为不在任何循环中，所以这部分代码根本就没有运行；第二，ELSE 之后没有语句了，如果的确没有，那么可以去掉。
其它的东西我就不懂了，那你的算法的问题了…… 另外，
& && &SV=(1-(T/Tc(N))**(0.5))**2
这里的 T 未定义。程序中 converged 也未赋初始值。
Warning 还是有一些的，比如一些未使用的变量…… 大哥，最后再说一次，以后再贴程序，麻烦用 code 模式好不？
code 模式对于程序的阅读和拷贝有很大的帮助…… Originally posted by snoopyzhao at
大哥，最后再说一次，以后再贴程序，麻烦用 code 模式好不？
code 模式对于程序的阅读和拷贝有很大的帮助…… C& &&&此程序是运用peng-robinson方程计算气体的逸度系数faita(1),faita(2)，是利用牛顿迭代法求解
C& &&&然后通过origin拟合得到这些数据N1,N2,m1,m2,K1,K2
C& &&&求解得到f1,f2,x1的过程,首先利用牛顿迭代法求解IAST理论得到x1, 然后利用表达式求解f1,f2
& && &PROGRAM MAIN
& && &LOGICAL converged
& && &integer S N SQ
& && &parameter(EPS=1.0E-06,R=83.14,epsilon=1-2**(0.5),sigma=1+2**(0.5)，N=7
& && &dimension faita(7),ns(4),ma(4),K(4),M(7),omiga(7),
& &&&*Tc(7),Pc(7),b(7),beta(7),alpha(7),q(7),as(7),Z1(7)
& && &real alpha,beta,Tr,omiga,M,rho(7),q,a,b,SV,I(7),ZM,SSA,SSB,SSC,
& &&&*SSD,T,f1,f2
& && &DATA M/16.043,30.070,44.097, 58.123,72.150, 44.01, 28.01/
& && &DATA omiga/0.012,0.100,0.152,0.2,0.252,0.224,0.048/
& && &DATA Tc/190.6,305.3,369.8,425.1,469.7,304.2,132.9/
& && &DATA Pc/45.99,48.72,42.48,37.96,33.70,73.83,34.99/
& && &READ(*,*)P1
C& &&&计算peng-robinson方程中需要的一些参量
& && &DO 3 N=1,7
& && &P=P1/100
& && &SV=(1-(T/Tc(N))**(0.5))**2
& && &alpha(N)=(1+(0.26*omiga(N)-0.26992*omiga(N)**2)*SV)
& && &as(N)=((0.45724*R**2*Tc(N)**2)/Pc(N))*alpha(N)
& && &b(N)=0.07779*R*Tc(N)/Pc(N)
& && &beta(N)=b(N)*P/(R*T)
& && &q(N)=as(N)*alpha(N)*P/(R**2*T**2)
3& &&&CONTINUE
C& &&&运用牛顿迭代法开始求解压缩因子的方程
& && &T=298.0
& && &IF(N.LT.7.AND.N.GT.1)THEN
& && &&&ZM=1
& && &&&DO WHILE(.not.converged.and.ABS(MSD(ZM,beta(N),q(N))).GT.0)
& && &&&SSA=MS(ZM,beta(N),q(N))
& && &&&SSB=MSD(ZM,beta(N),q(N))
& && &&&ZM=ZM-SSA/SSB
& && &&&converged=(abs(SSA)/(SSB)).LT.EPS
& && &&&END DO
& && &ELSE
& && &END IF
& && &do 4 N=1,7
& && &Z1(N)=ZM
4& &&&CONTINUE
C& &&&求解密度和逸度
& && &DO 5 N=1,7
& && &I(N)=(1/(sigma-epsilon))*log((Z1(N)+sigma*beta(N))/
& &&&*(Z1(N)+epsilon*beta(N)))
& && &rho(N)=(P/(Z1(N)*R*T))*M(N)*EPS
& && &faita(N)=SVF(Z1(N),beta,q(N),I(N))
5& &&&CONTINUE
C& &&&运用IAST理论以及上面得到的逸度系数，开始求解IAST的方程
& && &SSC=iast(X,P,y1,faita(1),faita(2))
& && &SSD=diast(X,P,y1,faita(1),faita(2))
& && &DO WHILE(.not.converged.and.ABS(SSD).GT.0)
& && &X=X-SSC/SSD
& && &S=S+1
& && &converged=(abs(SSC/SSD).LT.EPS)
& && &END DO
& && &x1=X
& && &f1=P*y1*faita(1)/x1
& && &f2=P*(1-y1)*faita(2)/(1-x1)
20& & format(1x,3F9.4)
& && &write(20,*)f1,f2,x1
& && &function MS(SA,SB,SC)
& && &MS=SA**3-(1-SB)*SA**2+(SC-2*SB-3*SB
& &&&***2)*SA-(SC*SB-SB**2-SB**3)
& && &function MSD(SD,SE,SF)
& && &MSD=3*SD**2-2*(1-SE)*SD+(SF-2*SE-
& &&&*3*SE**2)
& && &function& &iast(MADD,P,YY,SG1,SG2)
& && &dimension&&ns(4),ma(4),K(4)
& && &&&REAL MADD
& && &DATA&&ns/2.0,3.0,4.0,5.0/
& && &DATA ma/1.0,2.0,3.0,4.0/
& && &DATA K/5.0,6.0,7.0,8.0/
& && &f1=P*YY*SG1/MADD
& && &f2=P*(1-YY)*SG2/(1-MADD)
& && &iast=ns(1)*LOG(1+K(1)*f1**ma(1))+ns(2)*LOG
& &&&*(1+K(2)*f1**ma(2))-ns(3)*LOG(1+K(3)*f2**ma(3))+ns(4)*LOG(1+
& &&&*K(4)*f2**ma(4))
& && &function&&diast(MAD,S,Y,SG3,SG4)
& && &dimension&&ns(4),ma(4),K(4)
& && &&&REAL MAD
& && &DATA&&ns/2.0,3.0,4.0,5.0/
& && &data ma/1.0,2.0,3.0,4.0/
& && &data K/5.0,6.0,7.0,8.0/
& && &f1=S*Y*SG3/MAD
& && &f2=S*(1-Y)*SG4/(1-MAD)
& && &diast=(ma(1)*ns(1)*K(1)*f1*(ma(1)-1))/(1+
& &&&*K(1)*f1**ma(1))+(ma(2)*ns(2)*K(2)*f1*(ma(2)-1))/(1+K(2)*f1**ma(2))
& &&&*-(ma(3)*ns(3)*K(3)*f2*(ma(3)-1))/(1+K(3)*f2**ma(3))+
& &&&*(ma(4)*ns(4)*K(4)*f2*(ma(4)-1))/(1+K(4)*f2**ma(4))
& && &FUNCTION SVF(SG5,SG6,SG7,SG8)
& && &SVF=exp(SG5-1-log(SG5-SG6)-SG7*SG8)
& && &END Originally posted by zyj8119 at
C& &&&此程序是运用peng-robinson方程计算气体的逸度系数faita(1),faita(2)，是利用牛顿迭代法求解
C& &&&然后通过origin拟合得到这些数据N1,N2,m1,m2,K1,K2
C& &&&求解得到f1,f2,x1的过程,首先利用牛 ... 会了，谢谢提醒。。
var cpro_id = 'u1216994';
欢迎监督和反馈：本帖内容由
提供，小木虫仅提供交流平台，不对该内容负责。欢迎协助我们监督管理，共同维护互联网健康，如果您对该内容有异议，请立即发邮件到
联系通知管理员，也可以通过QQ周知，我们的QQ号为：8835100
我们保证在1个工作日内给予处理和答复，谢谢您的监督。
小木虫，学术科研第一站，为中国学术科研研究提供免费动力
广告投放请联系QQ： &
违规贴举报删除请联系邮箱： 或者 QQ:8835100
Copyright &
eMuch.net, All Rights Reserved. 小木虫 版权所有

我要回帖

更多关于 omiga欧米茄女款手表 的文章

 

随机推荐