#!/usr/bin/env python

import sys, commands
from numpy import array

T_vasp, Ek_vasp 	= [ ], [ ]

traj 	= 	sys.argv[1]
lista	=  	commands.getstatusoutput("grep T %s/dyn/OSZICAR | awk \'{print $3, $11}\'" %traj)[1].split() 
for i in range(1, int(float(len(lista))/2-1)):
	T_vasp.append(float(lista[2*i]))
	Ek_vasp.append(float(lista[2*i+1]))
#	print T_vasp[-1], Ek_vasp[-1]

T_script = [ float(x.split()[0]) for x in open('T_Ek.dat').readlines() ]
Ek_script = [ float(x.split()[1]) for x in open('T_Ek.dat').readlines() ]

T_script = array(T_script)
T_vasp = array(T_vasp)
Ek_script = array(Ek_script)
Ek_vasp = array(Ek_vasp)

print max(abs(T_script-T_vasp))
print max(abs(Ek_script-Ek_vasp))