%%%%%%%%%%
% ZrCX.m %
%%%%%%%%%%

nc = 8;
np0 = 8 * nc^3;
e0 = -7.998489;

a = [
    1.0       52888.056642  -32761.809480
    0.994629  52878.931716  -32663.510483
    0.989746  52870.606239  -32573.956087
    0.97998   52853.215893  -32394.976431
    0.969727  52834.628938  -32206.146745
    0.959961  52817.197908  -32026.683773
    0.949707  52799.367985  -31837.503996
    0.899902  52697.106422  -30914.019511
    0.799805  52461.962019  -29034.520904
    0.699707  52162.894278  -27115.761039
    0.599609  51793.716966  -25161.903944
    0.5       51331.832466  -23168.218821
    0.399902  50785.440874  -21133.436140
    0.299805  50093.529081  -19045.858821
    0.199707  49247.877743  -16893.924353
    0.0996094 48260.709394  -14700.425348
    0.0       47137.683567  -12461.815365
    ];

plot (1-a(:,1), a(:,2).^(1/3)/nc, 'o');
n = (a(:,1)+1)/2*np0;
de = (a(:,3) - n * e0) ./ (np0-n);

input ('Press a key... ');

plot (a(:,1), de, 'o');