function [y_predict] = k_hb_kaccum(x_train,y_train, x_predict, h)

n = size(x_predict, 1);
l=size(h,2);

y_predict = zeros(n,l);
c = 2 * h .^ 2;
 
K = 0;
for i = 1:n
    % quadratische Abweichung:
    % qd = bsxfun(@minus,x_predict(i,:),x_train)).^2;
    %
    %Summe der quadr. Abweichung, wenn x_train mehrdimensional
    % s= sum(qd,1);
    %
    %Bestimme für jedes h Wi, da quadratische Abweichung der Punkte
    %für alle h gleich bleibt, sich somit nur Wi ändert
    %Wi = bsxfun(@rdivide, -s, c);
    W = exp(-sum(bsxfun(@minus,x_predict(i,:),x_train).^2,2)/c);
    
    K = K + sum(W,1);
    
    if all(K == 0)
        y_predict(i) = 0;       
    else
        index=find(K==0);
        y_predict(i,index(:))=0;
        %wie kann ich hier sicher stellen, dass tatsächlich die operation
        %nur an den richtigen stellen durchgefürht wird??
        y_predict(i,:) = bsxfun(@rdivide,sum(bsxfun(@times, W, y_train),1) , K);
    end
end