function K=genmatrix(M,N)
% Funktion generiert die Matrix aus der Vorlesung mit M x N Gitterpunkten
% Nur Nichtnulleintraege werden gespeichert
size=(M-1)*(N-1);
K=sparse(size,size); 
v=sparse(M-1,1); % Vektor aus lauter Nullen
v(1)=2;
v(2)=-1; % v=[2,-1,0,...,0]
T1=M^2*toeplitz(v); % Matrix K_h, Skript Seite 30
w=sparse(N-1,1); %
w(1)=2;
w(2)=-1;
T2=N^2*toeplitz(w);
I1=speye(N-1); % Einheitsmatrix
I2=speye(M-1); % Einheitsmatrix
K=kron(T1,I1)+kron(I2,T2);


