一个关联某个未知函数本身及其各种阶数的衍生物的等式。
为什么 scipy.integrate.solve_ivp (RK45) 的准确性与我自制的 RK4 实现相比非常差？
我想解决什么问题 我使用 scipy.integrate.solve_ivp 求解了常微分方程组，并将结果与我自制的四阶龙格库塔 (RK4) 实现进行了比较。
用于求解联立 ODE 的并行化代码 (DifferentialEquations.jl) - Julia
我有以下 ODE 耦合系统（来自离散化积分微分偏微分方程）： xi 是我控制的 x 网格上的点。我可以用以下简单的 pi 来解决这个问题...
我这样做是为了使用射击方法解决BVP功能。 我不明白为什么当我的步长 h = 0.1 时，限制 x_end 超出了我在函数solve下的 while 循环中设置的值...
我尝试从头开始开发有机朗肯循环模型。热交换器是循环中唯一的动态子模型，每当我使用水作为介质时（无论是在...
我正在尝试求解这个微分方程组： 这些函数的行为应如下所示，即 r-> 无穷大 (F = g = lambda = const.)： 不久前我遇到了类似的问题，一位用户......
我正在尝试使用compute_orbit方法模拟行星的轨道，但是当我绘制结果位置时，我得到一条直线而不是预期的椭圆轨道。以下是重新...
如何在 SymPy 中按预期求解这个微分方程 (1-x^2)y'=x^2-xy-1？
我想用 Sympy 求解 (1-x^2)y'=x^2-xy-1 ，但无法按预期求解。 我期望得到以下解决方案。 -11:y(x)...
我想使用 scipy.odeint 在 Python 中求解常微分方程，同时从定义的函数返回一些值，这些值不是微分。 我尝试存储附加值...
非线性微分方程 (Korteweg–De Vries) 随着时间的推移变得不稳定
我想用初始条件 U(x,0) = n(n+1)/cosh^2(x) 求解 Korteweg-De Vries 方程。然而，当 n>1（n 是整数）时，我的系统在数值上变得不稳定，怎么办...
我有三个偏微分方程 (PDE) 和一个变量的解析解，如图所示。使用这些方程，我想求解 \phi(x,y,t)、p(x,y,t)、C_{a}(x,y,t) 和 C_{b}(x,y,t) ...
如何使用RK-4和射击法在Python中求解三个耦合微分方程？或者使用solve_bvp？
我正在尝试用Python求解三个耦合微分方程。我使用的是RK-4技术和射击方法。我正在尝试绘制 f 和 N 函数。 将 numpy 导入为 np 导入 matplotlib。
尝试根据以下模型文档模拟数学模型。 4.1.1.4 暴露于蒸气：蒸发 https://www.rivm.nl/bibliotheek/rapporten/2017-0197.pdf Python脚本： 导入
我对 Python 比较陌生，并尝试使用它来求解二阶非线性微分方程，特别是电解质中的泊松-玻尔兹曼方程。 phi''(r) + (2/r)*phi'(r) ...
我正在尝试模拟基于菲克第二定律的简单扩散。 从 pylab 导入 * 将 numpy 导入为 np 网格点 = 128 定义轮廓（x）： 范围 = 2。 散乱=.1576 剂量=1 重新...
尝试绘制图表时出现错误：TypeError：无法将表达式转换为浮点数 我还得到空图： 从 sympy 导入 * 从 IPython.display 导入显示 初始化打印() 导入matpl...
使用`solve_ivp`和`LSODA`来求解复杂的ODE
我正在尝试求解复杂的微分方程组。这些方程是刚性的，因此我需要使用一种可以处理复杂的 ODE 和刚度切换的方法。我已经降落在
DifferentialEquations.jl - 弹簧质量阻尼器系统 ODE
这是使用 DifferentialEquations.jl 在 julia 中编写弹簧质量阻尼器系统的正确方法吗？ 函数 smd(du, u, p, t) c、k、m = p du[1] = dx = u[2] du[2] = dv = -(c/m)*u[2] - (k...
我对 R 相当陌生，并且拥有这个 ODE 系统 rac{dS}{dT} = - eta(1-\mu)S(t) rac{I(t))}{1-D(t))} rac{dI}{dT} = - eta(1-\mu)S(t) rac{I(t))}{1-D(t))}-\delta I(t))-\phi I (t)) rac{dR}{...
使用 julia 的 NonlinearSolve 求解稳态时函数调用错误
我正在使用 NonlinearSolve 来解决稳态问题。当我做 prob = SteadyStateProblem(eoms!, σ₀, p) 溶胶 = 求解（ 概率， SSRrootfind(), abstol = slv.abstol, reltol = slv.relto...
TensorFlow 二阶导数计算在 Quantum PINN 求解器中始终返回零
我正在尝试实现一个简单的物理信息神经网络（PINN），用于在给定任意标量势的情况下求解与时间无关的薛定谔方程。我认为细节不是