n=50;
    alpha=1;
    A=gallery('poisson',n);

    K=[A -1/alpha*speye(n^2);
       speye(n^2) A];
    b=[zeros(n^2,1) -ones(n^2,1)];
    solution=Kb;
    %embedding with boundaries...