提交: 647f7300fbc055886066d17c5278bf4481604e97 - myprojects-hg-reborn (hg)

Commit MetaInfo

• 修订版: 647f7300fbc055886066d17c5278bf4481604e97 (tree)
• 时间: 2007-03-30 02:45:44
• 作者: iselllo
• Commiter: iselllo

Log Message

I added the code Python-codes/ode5_solver_test.py which solves a system
of 2 coupled equations for a trivial prey-predator system. Something
along these lines could be used to study and solve numerically the eqs
for nucleation and accumulation mode presented in the paper by Vouitsis.

更改概述

• add: Python-codes/ode5_solver_test.py

差异

diff -r 977e0ce0d715 -r 647f7300fbc0 Python-codes/ode5_solver_test.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/Python-codes/ode5_solver_test.py	Thu Mar 29 17:45:44 2007 +0000
@@ -0,0 +1,50 @@
+#! /usr/bin/env python
+from scipy import *
+import pylab  # used to read the .csv file
+
+
+
+
+#the purpose of this code is now to solve a system of 1st order ode's which are necessary for a prey-predator system
+
+
+def population(y,t,a,b,c,p):
+	
+	return [a*y[0]-b*y[0]*y[1],-c*y[1]+p*y[0]*y[1]]
+	
+
+a=2.
+b=1.e-1
+c=1.e-1
+p=1.e-2
+
+y0_0=15.
+y1_0=15.
+
+y0 = [y0_0, y1_0]
+
+
+t=arange(0.,100.,0.01)
+
+print 't is', t
+
+y = integrate.odeint(population, y0, t,args=(a,b,c,p),printmessg=1)
+
+
+
+
+pylab.plot(t,y[:,0],t,y[:,1])
+pylab.xlabel('Time')
+pylab.ylabel('y(t)')
+pylab.legend(('prey population','predator population'))
+pylab.title('solution')
+pylab.grid(True)
+pylab.savefig('preyandpredator')
+
+pylab.hold(False)
+
+
+
+
+
+print 'So far so good'
\ No newline at end of file