function fd3d(L)
% Implementation d'une méthode de différences finies en trois dimensions
% pour résoudre le problème du quiz 9 
  I = speye(L,L);
  subdiag = sparse(2:L,1:L-1,1,L,L);
  D = 2*I-subdiag-subdiag';
  A=???;
  N=L*L*L;
  h = 1./(L+1);
  b = h*h*ones([N,1]);
  x=A\b;
  norm(x,inf)
end

