%% Exercise 1: Advanced Plotting of ODP 659 Benthic Oxygen Isotope Record
clear
data = load('odp659data.txt');
plot(data(:,1),data(:,2))

%%
% Values/Properties of current figure window and axis system
get(gcf)
get(gca)

%%
% BlaetterwaldDesign Settings: Figure - White background, 550 px wide,
% 260 pixels high; Axis - x-axis 16 cm wide, y-axis 6 cm high,
% font Times 12 pt, color as you like but not default, line width 0.6 pt;
% Line - Line thickness of curve 0.5 pt, axis labels Helvetica 12 pt
% title Helvetica 12 pt, axis labels Helvetica 12 pt.
figure1 = figure('Color',[1 1 1],...
    'Position',[150 400 550 260]);
axes1 = axes('XLim',[0 6000],...
    'XDir','reverse',...
    'YLim',[2 5.5],...
    'YDir','reverse',...
    'Box','on',...
    'Units','Centimeters',...
    'Position',[2 2 16 6],...
    'LineWidth',0.6,...
    'FontName','Helvetica',...
    'FontSize',12);
line1 = line(data(:,1),data(:,2),...
    'Color',[0.8 0.3 0.1],...
    'LineWidth',0.5);
set(get(axes1,'Title'),...
    'String','ODP 659 Benthic Oxygen Isotope Record',...
    'FontName','Helvetica',...
    'FontSize',12)
set(get(axes1,'XLabel'),...
    'String','Age- Thousands of Years Before the Present',...
    'FontName','Helvetica',...
    'FontSize',12)
set(get(axes1,'YLabel'),...
    'String','d18O- Permille v. PDB',...
    'FontName','Helvetica',...
    'FontSize',12)

%%
print -depsc2 -cmyk odp659figure.eps