# 1802.py
execfile('../libreria/alfa.py')

def f (t,x): return x+t

def u (t): return math.exp(t)-t-1

#######################################
# Primo esperimento.

h=0.1; x=0
for n in xrange(21):
  t=n*h; print '%.1f %.6f %.6f' %(t,u(t),x)
  x+=h*rungekuttaphi(f,t,x,h)

#######################################
# Secondo esperimento.
print

h=0.01; x=0
for n in xrange(201):
  t=n*h
  if n%10==0: print '%.1f %.6f %.6f' %(t,u(t),x)
  x+=h*rungekuttaphi(f,t,x,h)