我有以下腳本來使用odeint來計算dRho。 P_r = 10e5
rho_r = 900
L = 750
H = 10
W = 150
A = H * W
V = A * L
fi = 0.17
k = 1.2e-13
c = 12.8e-9
mu = 2e-3
N = 50
dV = V/N
dx = L/N
P_in = P_r
rho_in = rho
這是我的代碼。它是一個函數,用於評估另一個函數在某個x值下的導數。即使對於分數階導數(a),我也希望它返回有效的輸出。 from scipy.special import gamma
import scipy.integrate as integrate
import sympy as sp
import scipy as sc
import math
def f(z):
r