利用python画微分方程组的图

13100661183

想利用odeint()函数解微分方程，做出四个图但是一直出现问题，快奔溃了，新手一枚，希望大神能指点一下，有点急。最后做出来的应该做出四个关于时间的函数图（一年内）

1. from scipy.integrate import odeint
   
2. import numpy as np
   
3. def CTANI(Y,t):
   
4.     a = 0.269
   
5.     b = 4.228
   
6.     c = 4.226
   
7.     d = 0.060
   
8.     e = 0.250
   
9.     f = 1.860
   
10.     return a*(b-c*exp(-t))+d*W-e*Y-f*Y/(Y+Q)
    
11.     def CNOXI(Q,t):
    
12.         g = 0.200
    
13.         h = 0.050
    
14.         f = 1.860
    
15.         return g*Y-h*Q-f*Y/(Y+Q)
    
16.     def CCHLI(Z,t):
    
17.         f = 1.860
    
18.         return f*Y/(Y+Q)
    
19.     def MNSEDI(W,t):
    
20.         i = 10.400
    
21.         d = 0.060
    
22.         return i*Z-d*W
    
23. 
    
24. time = linspace(0,365,1000)
    
25. Yinit = array([0.0,2.0])
    
26. Qinit = array([0.0,2.0])
    
27. Zinit = array([0.0,2.0])
    
28. Winit = array([0.0,2.0])
    
29. Y = odeint(CTANI,Yinit,time)
    
30. Q = odeint(CNOXI,Qinit,time)
    
31. Z = odeint(CCHLI,Zinit,time)
    
32. W = odeint(MNSEDI,Winit,time)
    
33. 
    
34. figure()
    
35. plot(time,Y[:,0],label='Y')
    
36. plot(time,Q[:,0],label='Q')
    
37. plot(time,Z[:,0],label='Z')
    
38. plot(time,W[:,0],label='W')
    
39. Ylabel('Y')
    
40. Qlabel('Q')
    
41. Zlabel('Z')
    
42. Wlabel('W')
    
43. legend()
    
44. show()
    

复制代码

13100661183

[img][/img]应该是这种图