function FEM1dExample3(n)

% -u''(x) + u(x) = 10 v (0,1)
% u(0) = 0
% u'(1) = 2

l = 1;
tau = 2;
k = 1 * ones(1,n);
alpha = 1 * ones(1,n);
f = 10 * ones(1,n);

[A,b] = FEM1d_plusMassMatrix(l,k,alpha,f,tau);
u = A\b;

x = linspace(0,l,n+1);
plot(x,[0,u'])