function [ ] = mef(n);
%
% Metodo dos elementos finitos para a resolucao do problema de condicoes
% de fronteira:  -u''(x) = f(x) com u(0) = u(1) = 0.
%
% Funcoes chamadas: f_P.
%
% L. N. Vicente, Novembro de 2004.
%
%
h = 1/(n+1);
%
for i=1:n,
    f_grelha(i) = f_P(h*i);
    solucao(i)  = -((h*i)^4)/12 + (h*i)/12;
end
%
e = ones(n,1);
A = (1/h) * spdiags([-e 2*e -e], -1:1, n, n);
%
b = h * f_grelha';
%
u = A \ b;
%
erro   = sqrt(h)*norm( u - solucao' );
limsup = h^2 * norm( [ 0 f_grelha 1/2 ] ) / (pi^2);
format = '\n Erro = %8.2e    Limite = %8.2e    h^2 = %8.2e \n\n';
h2     = h^2;
fprintf(format,erro,limsup,h2);
%
% Fim.