matthiasm@8: p1 = log(1e-5);
matthiasm@8: p2 = log(5*1e-6);
matthiasm@8: p3 = log(sum(exp([p1 p2])))
matthiasm@8: p4 = logsumexp([p1 p2],2)
matthiasm@8: p5 = logsum([p1 p2])
matthiasm@8: p6 = logsum([p1 p2])