function example3_2703(m)

% -(k(x) u'(x))' = 1, x \in (0,1)
%     k(0) u'(0) = 1,
%           u(1) = 2,
% kde k(x) = 1 pro x \in (0,1/2),
%     k(x) = 5 pro x \in (1/2,1).

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); k(round(m/2+1):end) = 5;
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')