SymPy求解微分方程初值问题时平方根使用报错求助

第一个报错：TypeError: Cannot convert expression to float #

• 原因：math.sqrt是Python标准库中处理数值计算的函数，只能接受整数、浮点数这类具体数值，而你传入的y(t)是SymPy定义的符号函数（属于符号表达式，不是具体数值），math.sqrt无法将符号表达式转换成浮点数，因此抛出错误。
• 解决：改用SymPy自带的sym.sqrt()函数，它专门用于处理符号表达式的平方根运算。

正确代码示例：

import sympy as sym
from sympy import display_latex

t = sym.symbols('t')
y = sym.Function('y')
eq3c = sym.Eq(y(t).diff(t), t - 2*sym.sqrt(y(t)))
print("The equation")
display_latex(eq3c)
eq3c_sol = sym.dsolve(eq3c, y(t), ics={y(0):1})
display_latex(eq3c_sol)

第二个报错：NotImplementedError: The given ODE... cannot be solved by the factorable group method #

• 原因：你用了(y(t))**(1/2)，这里的1/2是Python的浮点数（0.5），SymPy在处理浮点数指数时，会将整个表达式转换成数值型近似形式，干扰了SymPy符号求解器对ODE类型的识别，导致求解器找不到合适的解法。
• 解决：使用SymPy的有理数类型表示1/2，比如sym.Rational(1,2)或sym.S(1)/2，或者直接用sym.sqrt(y(t))（更直观，推荐）。

正确代码示例（两种写法可选）：

import sympy as sym
from sympy import display_latex

t = sym.symbols('t')
y = sym.Function('y')
# 写法一：用sym.sqrt
eq3c = sym.Eq(y(t).diff(t), t - 2*sym.sqrt(y(t)))
# 写法二：用有理数指数
# eq3c = sym.Eq(y(t).diff(t), t - 2*(y(t))**sym.Rational(1,2))
print("The equation")
display_latex(eq3c)
eq3c_sol = sym.dsolve(eq3c, y(t), ics={y(0):1})
display_latex(eq3c_sol)

运行后可得到正确解析解：$y{\left(t \right)} = \left(\frac{t}{2} - 1\right)^{2}$

内容的提问来源于stack exchange，提问作者username_1326