Dens = [1.00747,1.00707,1.00409,0.99908; ...
        1.01509,1.01442,1.01112,1.00593; ...
        1.03038,1.02920,1.02530,1.01977; ...
        1.06121,1.05907,1.05412,1.04798; ...
        1.09244,1.08946,1.08365,1.07699; ...
        1.12419,1.12056,1.11401,1.10688; ...
        1.15663,1.15254,1.14533,1.13774; ...
        1.18999,1.18557,1.17776,1.16971; ...
        1.20709,1.20254,1.19443,1.18614];
WtPct = [1,2,4,8,12,16,20,24,26];
Temp = [0,10,25,40];
[WP,TP] = meshgrid(WtPct,Temp);
WP = WP' ; TP = TP';
contour(WP,TP,Dens,'ShowText','on')
grid
xlabel('Wt Pct NaCl')
ylabel('Temperature - degC')
title('Density of NaCl Solutions')

