我试图在Sage笔记本上绘制斜坡场和几个独特的解决方案,但是我的斜坡场与解决方案不匹配。 c = var('c')
y = (ln(x)+c)/x
show(y)
p1 = plot_slope_field(y, (x,0,5), (c,-5,5))
for i in range(-5,5):
p1 = p1 + plot(y(c=i),x,0,5, ymin=-5, ymax
XMDS是一个软件包,它允许通过将它们编译为C并以本地机器体系结构允许的最大速度运行它们来解决普通,部分和随机微分方程。 通常情况下,XMDS支持带有GOMP OpenMP库的GCC编译器。 我试图安装XMDS,但每当我运行它,我得到以下错误: Checking for program mpic++ :
/Users/dg813/Downloads/XMDS.app/Conten
我第一定义函数为DY/DT = Y和dy/dt的= T: function dy=d(y):
dy=y
end
function ddy=dd(t):
ddy=t
end
然后我分别使用ODE45,: [t,y]=ode45('d',[1 10],1)
[t,y]=ode45('dd',[1 10],1)
返回以下错误:使用d 太多输入参数出错。 我的问题是:
我试图解决的4个偏微分方程耦合非线性系统,下面的代码: import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint
from odeintw import odeintw
# constants
global rH, L, q, p, Q, mu
rH = 1.0 ; L