%
%  Risolve il problema BVP
%
%  y'' + (1/epsi) y' = (1/epsi)
%  y(0) = y(1) = 0
%  
%  con diferenze in avanti per il termine di trasporto (y').
%  Il sistema risultante
% 
%  +                                            + +        +   +        +
%  |  alpha(1)  gamma(1)                        | | y(1)   |   | b(1)   |
%  |  beta(2)   alpha(2)  gamma(2)              | |        |   |        |
%  |              .......                       | |  ..    | = |        |
%  |                                            | |        |   |        |
%  |                       beta(N-1) alpha(N-1) | | y(N-1) |   | b(N-1) |
%  +                                            + +        +   +        +
%
%  E' risolto con le primitive di MATLAB
%
function [x,y]=BLsolve1(epsi,N)
  h = 1.0/N ;
  x = [h:h:1-h] ;
  % riempimento del vettore alpha
  al = -(2+h/epsi).*ones(1,N-1);
  % riempimento del vettore beta
  be = ones(N-1) ;
  % riempimento del vettore gamma
  ga = (1 + h/epsi).*ones(1,N-1) ;
  % riempimento del vettore r
  bvec = h^2.*ones(1,N-1)./epsi ;
  %
  A = sparse(diag(al(1:N-1),0)) + ... 
      sparse(diag(be(2:N-1),-1)) + ...
      sparse(diag(ga(1:N-2),+1)) ;
  y = A\bvec' ;
end
