常微分方程组求数值解python+matlab

python实现

常微分函数表达式：

import scipy.integrate as spi
import numpy as np
from matplotlib import pyplot as plt
import warnings
warnings.filterwarnings('ignore')
# 参数
ATR=0.2
ATAR=0.3
SR=0.1
SAR=0.3
SRC=0.7
PRC=0.2
PNP=0.3
AN=0
VR=SR+ATR+AN
VAR=SAR+ATAR+PRC*SRC
VI=PNP
a12=(VR-VAR)/3
a21=(-VR+VAR)/3
a31=(-VR+VI)/3
a32=(-VAR+VI)/3
RS=0.25
ARS=0.1
I=0.65
y=[RS,ARS,I]
k=9
w1=(k*k+k)/((k+3)*(k-2))#2
w2=3/((k+3)*(k-2))#0.5
w3=(k*k+k-3)/((k+3)