f = fopen("EDC.TXT", "r");
a = fscanf(f, "%g", [5, Inf]);
fclose(f);
time = a(2,:);
temperature = a(4,:);
plot(time, temperature);
xlabel('Thousands of years ago');
ylabel('Temperature Anomaly T(t)-T(1970)');
axis tight;