function example2_2703(m)

P = linspace(0,1,m+1);
E = [1:m;2:m+1];
ID = [m+1]; uhat = [2];
IN = [1]; tauhat = [-1];
k = ones(1,m);
f = ones(1,m);

[A,b,utilde] = FEM(P,E,{'aloc_1dLaplace'},{k},{'bloc_1dsource'},{f},ID,uhat,IN,tauhat);

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

plot(P,u')