我試圖在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