u1 = [0;0;0];
p10 = [0;1;0];
w1 = [0;0;1];
w1daszek = [0, -1, 0; ...
            1, 0, 0; ...
            0, 0, 0];
syms t real
expm(w1daszek*t)