如何为具有奇点点的函数绘制图
我正在尝试为以下一阶微分方程绘制解决方案曲线（u〜V）：（ u/v）（dv/du）=（1-u）/（u++v-3）， 该方程具有边界条件作为： v = 0 at u = 3 和 DV/D…
SOLVE_IVP错误：＆＃x27;所需的步长尺寸小于数字之间的间距。＆＃x27;
我有15个ODESTO可以同时解决，并在deriv1函数中定义。 T，CO2和Q有5个状态。初始条件为t = 20，CO2 = 0，q = 0。当T_eval设置为[1、2、3、4、5]时，我…
scipy中的solve_ivp中有什么事件
我试图了解 Events 参数在 solve_ivp 中正在做什么，以及该功能会花费时间以及其他功能作为其参数。 另外，在功能为0时，内部步骤的数量和大小会发生…
如何强制solve_ivp仅允许集成变量的积极解决方案？
随着时间的流逝，我正在对表面上的冰冻进行建模。这需要计算“冷冻/熔化速率”，这可能大于或小于零。 我要施加的约束是，如果表面上的冷冻冰的量为零…
在采样强制微分方程时，朱莉娅/图林的typeserror
我是朱莉娅（Julia）和图灵（Turing）的新手，并且正在尝试将强制的0-D框颂歌适合数据，但是进行采样时会遇到类型错误。 遵循此页面（在外部强迫），…
如何在Sigle线中找到两个常见的切线？如何求解二进制微分方程？
我有一个名为 f 的函数，如下所示： 两孔线 我想求解方程，以获取给定范围内函数的两个常见的切线（x1，f（x1）），（x2，f（x2））。 方程是： F'(x1…
我如何简化代码以获取Jacobian Pass
我正在尝试计算一个5x12矩阵，该矩阵包含相对于12个参数的5个函数的衍生物（S，D，Beta，...）。我正在使用Jacobian来计算事物，但是我在此处遇到此错…
用Dirichlet和Neumann边界条件求解二阶微分方程
热方程式的傅立叶定律 我想解决一个孤立的电加热杆的 ：具有dirichlet边界条件的条件 和neumann边界的条件，其中 x 的长度坐标为 杆的长度是杆k的长度…
步骤加倍runge kutta实现卡住了缩小的步骤
我需要使用Audaptive RK4方法与步骤倍加倍技术集成ODES系统。 问题在于该程序将永远继续缩小到机器精确度的同时，同时又不延长时间。 这个想法是通过…
Python模块Gekko可以处理复杂数字吗？
我面对Hessenberg Index-2 DAE，我正在尝试使用Python模块Gekko解决它。 经过几天的反复试验，我认为我离有效的代码不远。但是我刚刚意识到Gekko也许…
