1""" 2 Copyright 2011 Mario Mulansky 3 Copyright 2012 Karsten Ahnert 4 5 Stochastic euler stepper example and Ornstein-Uhlenbeck process 6 7 Distributed under the Boost Software License, Version 1.0. 8(See accompanying file LICENSE_1_0.txt or 9 copy at http://www.boost.org/LICENSE_1_0.txt) 10""" 11 12 13from pylab import * 14from scipy import special 15 16data1 = loadtxt("elliptic1.dat") 17data2 = loadtxt("elliptic2.dat") 18data3 = loadtxt("elliptic3.dat") 19 20sn1,cn1,dn1,phi1 = special.ellipj( data1[:,0] , 0.51 ) 21sn2,cn2,dn2,phi2 = special.ellipj( data2[:,0] , 0.51 ) 22sn3,cn3,dn3,phi3 = special.ellipj( data3[:,0] , 0.51 ) 23 24semilogy( data1[:,0] , abs(data1[:,1]-sn1) ) 25semilogy( data2[:,0] , abs(data2[:,1]-sn2) , 'ro' ) 26semilogy( data3[:,0] , abs(data3[:,1]-sn3) , '--' ) 27 28show() 29 30 31 32