% temp = 0.95;
temp = 0.8;

% 
% Code begins:
% --

% close all;

VOL_MAX   = 5;
N_POINTS  = 1000;
DELTA_VOL = 0.2*(1/3);
vol   = linspace( (1/3)+DELTA_VOL, VOL_MAX, N_POINTS ); 
% clear VOL_MAX N_POINTS DELTA_VOL;
press = 8*temp./( 3*vol - 1 ) - 3./(vol.^2);
gibbs = -log( vol - 1/3 ) + (1/3)./( vol - 1/3);

figure; 

plot( vol, press, '-o' ); grid on;

figure; 

plot( vol, gibbs, '-o' ); grid on;