function conditionNumber_2dFEM()

a = 1; hs = 0.5*2.^[0:-1:-6];

for i=1:length(hs)
    h = hs(i);
    [P,T,ID,IN] = discretizeRectangle(a,a,h);
    N = size(P,2)
    m = size(T,2);
    nD = length(ID); uhat = zeros(nD,1);
    mN = size(IN,2); tauhat = zeros(mN,1);

    k = ones(1,m);
    f = ones(1,m);
    [A,b,utilde] = FEM(P,T,{'aloc_2dLaplace'},{k},{'bloc_2dsource'},{f},ID,uhat,IN,tauhat);

    conds(i) = cond(full(A));
end
plot(log10(hs),log10(conds),'r',log10(hs),log10(hs.^(-2)),'r--',log10(hs),log10(hs.^(-1)),'b--')
