TempC = [0,4.4,10,15.6,21.1,26.7,32.2,37.8,48.9, ...
         60,71.1,82.2,93.3];
Visc = [1.794,1.546,1.31,1.129,0.982,0.862,0.764, ...
        0.682,0.559,0.47,0.401,0.347,0.305];
plot(TempC,Visc,'sk')
grid
xlabel('Temperature - degC')
ylabel('Viscosity - cP')
title('Viscosity of Water vs. Temperature')
