function example4_2703(m)

% -u''(x) + 0.1 u(x) = 1, x \in (0,1),
%               u(0) = 2,
%               u(1) = 3.

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

[A,b,utilde] = FEM(P,E,{'aloc_1dLaplace','aloc_1dMass'},{k,q},{'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')