clear Wk*
Vk = (w' .* (V(x) + VH(x))) * ones(1,N);

Wk = zeros(N);
tic;
for k1=1:N
    Wk(k1,:) = B(k1,:) * (Vk .* B.');
end
fprintf('Referenz: %s\n', szeit(toc));

% ohne for Schleife
tic;
Wk1 = B * (Vk .* B.');
fprintf('ohne for Schleife:   %s    ', szeit(toc));
serr(Wk, Wk1);

% mit Sparse
Vks = diag(sparse(w' .* (V(x) + VH(x))));
tic;
Wk2 = B * Vks * B.';
fprintf('mit Sparse Struktur: %s    ', szeit(toc));
serr(Wk, Wk2);

whos Vk Vks;