function da = dpl(l1, l2)
ca = (l1'*l2)/(sqrt(sum(l1.^2))*sqrt(sum(l2.^2)));
da = abs(acos(ca));