提交: 4f04f941fd6963a4874023ca5925fc8d5292b6b0 - myprojects-hg-reborn (hg)

Commit MetaInfo

• 修訂: 4f04f941fd6963a4874023ca5925fc8d5292b6b0 (tree)
• 時間: 2007-03-23 00:27:38
• 作者: iselllo
• Commiter: iselllo

Log Message

I added the code Python-codes/ode4_solver_test.py. It is a very easy
examples about how to use the Python integrate.odeint to solve this time
for a 1st order ODE. As shown in Python-codes/ode3_solver_test.py,
higher order ODE's are solved as a system of 1st order ODE's.

Change Summary

• add: Python-codes/ode4_solver_test.py

差異

diff -r b7bae9fe5845 -r 4f04f941fd69 Python-codes/ode4_solver_test.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/Python-codes/ode4_solver_test.py	Thu Mar 22 15:27:38 2007 +0000
@@ -0,0 +1,36 @@
+#! /usr/bin/env python
+from scipy import *
+import pylab  # used to read the .csv file
+
+
+from scipy.special import * # to be able to use the error function
+
+def func(y, t):
+	return y+3.*t
+
+
+y0 = 1.
+
+t = arange(0,4.0, 0.01)
+
+
+y = integrate.odeint(func, y0, t,printmessg=1)
+
+
+print 'the shape of the solution is', shape(y)
+
+
+pylab.plot(t,y[:,0])
+pylab.xlabel('Time')
+pylab.ylabel('y(t)')
+pylab.title('solution')
+pylab.grid(True)
+pylab.savefig('solution_1st_order_ODE')
+
+pylab.hold(False)
+
+
+
+
+
+print 'So far so good'
\ No newline at end of file