from scipy.integrate import odeint
from scipy import arange
#from array import array


tauvalues=range(3,11)



fudge=10000000
t=arange(0,10**10 + fudge,fudge)


def f(y,x):
	k1=10**-9
	k2=10**9
	return [-k1*y[0],k1*y[0] - k2*y[1],k2*y[1]]

f0=[1,0,0]

def answers(sol):
	nfe=sol[1]['nfe']
	numeval=nfe[len(nfe)-1]
	#get answer at endpoint
	b=sol[0].tolist()
	b=b[len(b)-1]
	#give it back
	return b,numeval

for i in tauvalues:
	tau=10**-i
	soln=odeint(f,f0,t,rtol=tau,full_output=1)
	print answers(soln)