sympy是一个非常好用的基于python的符号计算库,科技做微分、积分、极限等一系列高等数学运算。
SymPy is a Python library for symbolic mathematics. It aims to become a full-featured computer algebra system (CAS) while keeping the code as simple as possible in order to be comprehensible and easily extensible. SymPy is written entirely in Python. https://www.sympy.org/en/index.html
例如验证等熵过程的熵变为0。本文python环境安装了anaconda全家桶,并在jupyter notebook中运行,首先导入计算库sympy:
import sympy
为了让公式更好看,使其排版如同?????LATEX,我们如下设置:
from sympy import init_printing
init_printing(use_latex=True)
好了定义一些变量,比如熵变,压强等等:
Cp,Rg,T1,T2,P1,P2,kappa = sympy.symbols('Cp Rg T1 T2 P1 P2 kappa')
dS = Cp*sympy.ln(T2/T1)-Rg*sympy.ln(P2/P1)
dS
运行程序,结果如下图,可见公式排版杠杠的:
#又根据等熵过程有P2/P1为:
ratio=sympy.Pow(T2/T1,kappa/(kappa-1))
ratio
#Cp与kappa的关系为:
kappa*Rg/(kappa-1)
稍微整理,可知熵变为0
早上喝了一碗羊汤,要求开发票!老板问:"啥单位?"我说:"中华人民共和国内蒙古自治区巴彦淖尔市乌拉特右翼后旗乌兰哈达苏木那令格嘎查那令格村那令格二社村民委员会火山历史文化遗址申报联合国教育科学文化组织世界非物质文化遗产工作领导小组办公室主任乌力图乌力吉吉尔格拉.巴彦尔太!"老板说:"羊汤钱不要了!