function sj = gp(j)



j = 1/norm(j)*j;

sx = [0 1; 1 0];

sy = [0 -i; i 0];

sz = [1 0; 0 -1];

sj = j(1)*sx+j(2)*sy+j(3)*sz;