load ych.dat;
data = log(ych);
oy = data(:,1);
oc = data(:,2);
oh = data(:,3);