function plotSol(N,Z,EQ)
  [X,Y] = meshgrid(linspace(0,1,N+1),linspace(0,1,N+1)) ;
  surf(X,Y,Z) ;
