Python 当我用Sympy积分一个函数时,它返回一个零

Python 当我用Sympy积分一个函数时,它返回一个零,python,math,sympy,calculus,Python,Math,Sympy,Calculus,我写这段代码是为了找到管状化学反应器中的反应分馏过程 来自sympy.interactive导入打印 从sympy导入限制、限制、符号、S 从辛输入积分,积分 printing.init_printing(使用_latex='mathjax') f=符号('f') def压裂(fx、V、Q): gfx=[] 对于范围(0100)内的i: ff=i/100 n=[积分(fx,(f,0.0,ff)).doit(),ff] gfx.append(n) c=0 对于gfx中的i: 如果我[0] fro

我写这段代码是为了找到管状化学反应器中的反应分馏过程

来自sympy.interactive导入打印
从sympy导入限制、限制、符号、S
从辛输入积分,积分
printing.init_printing(使用_latex='mathjax')
f=符号('f')
def压裂(fx、V、Q):
gfx=[]
对于范围(0100)内的i:
ff=i/100
n=[积分(fx,(f,0.0,ff)).doit(),ff]
gfx.append(n)
c=0
对于gfx中的i:
如果我[0]
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)