wm = mean(w(2:10,:))
wdm = w - wm(ones(1,40),:);
plot(wdm)
plot(1:40, mean(wdm.'), 1:40, median(wdm.'), 1:40, std(wdm.')); grid on
legend('mean','median','std')
title('hale3088 W - mean of cells 2-10')
print -dwinc