SYmpy以前曾在各个方面工作过，但我为是否有解决方案感到惊讶。我把它留了一整夜，Python还没有给我解决方案。有没有一种方法可以加快这个过程，或者，我在这里缺少什么。
from sympy import *
import numpy as np
import math
rho_l = 1352;
rho_g= 1.225;
sigma = 0.029;
nu = 0.02;
Q = rho_g/ rho_l;
u = 99.67;
h = 1.6e-5; # Half sheet thickness
# Actual Dispersion Equation for viscous sheets
k = Symbol('k', real=True)
t = tanh(k*h);
w1 = -2 * nu * k**2 * t ;
w2 = 4* nu**2 * k**4 * t**2;
w3 = - Q**2 * u**2 * k**2;
w4 = -(t + Q)
w5 = (- Q* u**2 * k**2 + (sigma*k**3/ rho_l));
w6 = -w4;
w = ( w1 + sqrt(w2 + w3 + w4*w5))/w6;
print (w)
wprime = w.diff(k)
print (wprime)
ko = solve(wprime, k) # solve fprime = 0 with respect to x
ko = ko[0]
print (ko)
答案 0 :(得分：1)
在这种情况下尝试nsolve：
>>> nsolve(wprime, 400)
416.221161431339