cc = 4*ones(n+1);  % coefficients in Laplace
        ce = ones(n+1);
        cw = ones(n+1);
        cn = ones(n+1);
        cs = ones(n+1);
        cw(2,2) = 0.;
        cc(2,2) = 2.;
        cs(2,2) = 0.;
        cs(3:n-1,2) = 0.;
        cc(3:n-1,2) = 3.;
        cs(n,2) = 0.;
        cc(n,2) = 2.;
        ce(n,2) = 0.;
        ce(n,3:n-1) = 0.;
        cc(n,3:n-1) = 3.;
        ce(n,n) = 0.;
        cc(n,n) = 2.;
        cn(n,n) = 0.;
        cn(3:n-1,n) =0.;
        cc(3:n-1,n) = 3.;
        cw(2,n) = 0.;
        cc(2,n) = 2.;
        cn(2,n) = 0.;
        cw(2,3:n-1) = 0.;
        cc(2,3:n-1) = 3.;