clear all 
close all
clc

load data.mat


        figure(1001)
        h = pcolor(nn,RR,-dc');
        set(h, 'EdgeColor', 'none');
        set(h,'alphadata',double(~isnan(dc)))
        colorbar
        naslov = (['\Delta X = X_1 - X_2, alpha = ' num2str(alpha)]);
        title(naslov)
        xlabel 'n'
        ylabel 'R = Q/M'
        axis equal
        axis([0 1 0 1])
        
print -deps fig1.eps
print -depsc fig1.epsc
print -depsc2 fig1.epsc2
