我写这段代码是为了找到管状化学反应器中的反应分馏过程
from sympy.interactive import printing
from sympy import Limit, limit, Symbol, S
from sympy import Integral, integrate
printing.init_printing(use_latex='mathjax')
f = Symbol('f')
def frac(fx,V,Q):
gfx=[]
for i in range(0,100):
ff=i/100
n=[Integral(fx, (f,0.0,ff)).doit(),ff]
gfx.append(n)
c=0
for i in gfx:
if i[0]<=(V/Q):
c=c+0.01
return(c-.01)
V=90
Q=17
K=0.0554
cde=4.897
fx=(-cde/(-K*(cde*(1-f))**2))
fr=frac(fx,V,Q)
如果单击图像链接,您将看到所有积分值都是零
你怎么看,保存积分值的列表都是零,我试着使用SymboLab和Wolfram Alpha,它们都给其他值不同于零
为什么Symphy给了我一个错误的函数和这些值的积分
目前没有回答
相关问题 更多 >
编程相关推荐