function H = laws_kernel(k)
    
    L=[1;2;1];
    E=[1;0;-1];
    S=[1;-2;1];
    lc=1/6;
    esc=1/2;

    if k>=1 && k<=9
        if k==1
            H=lc*L*lc*L';
        elseif k==2
            H=lc*esc*L*E';
        elseif k==3
            H=lc*esc*L*S';
        elseif k==4
            H=esc*lc*E*L';
        elseif k==5
            H=esc*esc*E*E';
        elseif k==6
            H=esc*esc*E*S';
        elseif k==7
            H=lc*esc*S*L';
        elseif k==8
            H=esc*esc*S*E';
        elseif k==9
            H=esc*esc*S*S';
     
    else
        error();
        end
        
end