function r = cross(L,r0)
r = zeros(3,1);
r(1) = L(2)*r0(3)-L(3)*r0(2);
r(2) = L(3)*r0(1)-L(1)*r0(3);
r(3) = L(1)*r0(2)-L(2)*r0(1);