function example5_0304(m)

% -k*u''(x) + 2*alpha/r*u(x) = f, x \in (0,1),
%       alpha*u(0) - k*u'(0) = alpha*uhat,
%      -alpha*u(1) - k*u'(1) = -alpha*uhat,
% kde k=alpha=f=uhat=1, r=0.1.

P = linspace(0,1,m+1);
E = [1:m;2:m+1];
ID = []; uhat = [];
alpha = 1; uhat = 1;
% % alpha = 0;
IN = [1,m+1]; tauhat = alpha*uhat*[1;1];
k = ones(1,m);
r = 0.1; k0 = 2*alpha/r * ones(1,m); % k0(:) = 0;
f = ones(1,m);

[A,b,utilde] = FEM(P,E,{'aloc_1dLaplace','aloc_1dMass'},{k,k0},{'bloc_1dsource'},{f+k0*uhat},ID,uhat,IN,tauhat);
A(1,1) = A(1,1) + alpha;
A(m+1,m+1) = A(m+1,m+1) + alpha;

u = A\b; % % pinv(full(A))*b;

plot(P,u')