提交: 321b69036a418c0a078602292c90f1d0be823d17 - myprojects-hg-reborn (hg)

Commit MetaInfo

• 修訂: 321b69036a418c0a078602292c90f1d0be823d17 (tree)
• 時間: 2007-03-21 05:54:23
• 作者: iselllo
• Commiter: iselllo

Log Message

I added the code Python-codes/ode2_solver_test.py, which is taken from
the numpy tutorial. This code basically performs the integration of the
Airy equation. It is taken from the scipy tutorial, but I am using
different boundary conditions wrt the example.

Change Summary

• add: Python-codes/ode2_solver_test.py

差異

diff -r 2957f3555ae8 -r 321b69036a41 Python-codes/ode2_solver_test.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/Python-codes/ode2_solver_test.py	Tue Mar 20 20:54:23 2007 +0000
@@ -0,0 +1,38 @@
+#! /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 [t*y[1],y[0]]
+
+y0_0=1.
+y1_0=1.
+
+y0 = [y0_0, y1_0]
+
+x = arange(0,4.0, 0.01)
+
+t=x
+y = integrate.odeint(func, y0, t)
+
+
+print 'the shape of the Lorenz solution is', shape(y)
+
+
+pylab.plot(x,y[:,1])
+pylab.xlabel('Time')
+pylab.ylabel('y(t)')
+pylab.title('solution')
+pylab.grid(True)
+pylab.savefig('airy_modified')
+
+pylab.hold(False)
+
+
+
+
+
+print 'So far so good'
\ No newline at end of file