function c=Galerkin_ef_1(inter,bv,n)
% Elementos finitos de ejemplo clase
%   inter, intervalo de integración; bv, valores de contorno;
%   n, número de pasos.
%
 a=inter(1); b=inter(2); ya=bv(1); yb=bv(2);
 h=(b-a)/(n+1); 
 alfa=(8/3)*h+2/h; beta=(2/3)*h-1/h;
 e=ones(n,1); 
 M=spdiags([beta*e alfa*e beta*e],-1:1,n,n);
 d=zeros(n,1); d(1)=-ya*beta; d(n)=-yb*beta;
 c=M\d;
end