


// matrice tri diagonale du laplacien en dim 1, differences finies centrees
//  periodique  coeff*[-2 1; 1 -2 1;....., 1 -2] de taille M*M

function [B]=laplacien(coef,M)

 B=zeros(M,M);
 // diagonale
 for i=1:M,
  B(i,i)=-2*coef;
 end;
 //sur diagonale
 for i=1:M-1,
  B(i,i+1)= coef;
 end;
 // sous diagonale
 for i=2:M,
  B(i,i-1)=coef;
 end;

endfunction;

