my.write(str(t) + '\t' +str(x) + '\t' + str(y) + '\t' + str(z) + '\t' + str(vx) + '\t' + str(vy) + '\t' + str(vz) +'\n')

renge_kutte()
my=open('motion3','r')
def motion():
	xs=[]
	ys=[]
	zs=[]
	for line in my:
		stolbci=line.split()
		x=float(stolbci[1])
		y=float(stolbci[2])
		z=float(stolbci[3])
		xs.append(x)
		ys.append(y)
		zs.append(z)
	return xs, zs, ys
import pylab
from mpl_toolkits.mplot3d import Axes3D
xs, ys, zs = motion()
fig = pylab.figure()
axes = Axes3D(fig)
axes.plot(xs, ys, zs)
pylab.show()

def write(*args):
    print "\t".join(args) + "\n"

print >>fp, s