function  qdata = sigmatr(data)
qdata = data(:,1) + data(:,2);