1function S = skew(vec)
2
3S = [0, -vec(3), vec(2); vec(3), 0, -vec(1); -vec(2), vec(1), 0];
4
5end
6