function example1_2d()

a = 1; b = 1; h = 0.1;
[P,T,ID,IN] = discretizeRectangle(a,b,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);

uH = A\b;
u = utilde;
nonDirNodes = setdiff(1:N,ID);
u(nonDirNodes) = u(nonDirNodes) + uH;

trisurf(T',P(1,:),P(2,:),u')
