//initialiser le vecteur b

function b = initialiseVecteurb(n)

b = zeros(n,1);
  
  b(1,1) = 1;
  b(n,1) = 1;
  
endfunction
