数学实验与建模 第八章 微元分析法

在学习例题8.8时，有一个细节没弄懂，期间也查了资料，写在这里，与大家分享。

import sympy as sp
sp.var('h')
sp.var('t',cls=sp.Function)
g=9.8
eq=t(h).diff(h)-10000*sp.pi/0.62/sp.sqrt(2*g)*(h**(3/2)-2*h**(1/2))
t=sp.dsolve(eq,ics={t(1):0})
t=sp.simplify(t)
print(t.args[1].n(9))

关于代码中的最后一句

print(t.args[1].n(9))

有一位博主做了解释，但是我没能理解他的解释。

https://blog.csdn.net/weixin_45870904/article/details/113100973

经过分析后，我觉得最后一段代码的含义如下：

比如你运行一个程序:测试.exe test。args[0]是你程序的名称(测试)，args[1]是你运行的程序带的第一个参数test

如果我们把最后一句改成print(t)，则输出如下：

Eq(t(h), -15260.5042485347*h**1.5 + 4578.15127456042*h**2.5 + 10682.3529739743)

因此说，args[1]代表的是取出符号表达式。

n(9)的含义是系数取前9位，如果我们把最后一句改成print(t.args[1])，则输出如下：

-15260.5042485347*h**1.5 + 4578.15127456042*h**2.5 + 10682.3529739743

如果加上.n(9)，则输出如下：

-15260.5042*h**1.5 + 4578.15127*h**2.5 + 10682.353

经过对比可以看出，n(9)的意思就是保留系数的前9位。