function M=plotNS(W)
%% W=3D-Matrix von Vortizitaet Schnappschuss
%% M=Matlab movie

function M=plotNS(W)
vecMax=16; % # auf Punkten um Pfeile zu zeichnen
N=size(W,3);
[n,m]=size(W(:,:,1));
if n<=vecMax,
  % Dx und Dy in Fourier
  % Laplace-Operator in Fourier                             
  % Laplace-Loeser
end
x=0:2*pi/n:2*pi*(n-1)/n;

for j=1:N,
  image(x,x,W(:,:,j),'CDataMapping','scaled'); 
  if n<=vecMax, % um Pfeile zu zeichnen      
    % Berechnung von w-hut:
    what=fft2(W(:,:,j));
    % Berechnung von u und v
    hold on;
    % Plot von Pfeilen:
    quiver(x,x,v,u);
    hold off;
  end
  xlabel('y');ylabel('x');
  M(j)=getframe; % um einen Film zu machen
  % Benutzen Sie movie(M), um den Film zu spielen
end
