load iers.out -ascii
mjd=iers(:,3);
polx=iers(:,4);
poly=iers(:,5);
dut1=iers(:,6);
tai=iers(:,7);
plot(mjd,tai,'r')
xlabel('Modified Julian Date')
ylabel('TAI-UT1 in s')
title('TAI-UT1 from EOP97C04')
grid