function Y = harmonics(L, M, x) 

    MM = [0 : L]';
    %MM

    wsp = (1/2) * sqrt((2 * L + 1) .* fact(L - MM) ./ (fact(L + MM)* pi));
    %wsp
    wsp = wsp * ones(size(x));
    %wsp
    
    Y = legendre(L, x) .* wsp;
    %Y
    % size(Y)
    
    Y2 = Y(2 : L + 1, :);
    %Y2
    
    % co drugi * (-1)
    Y2 = (cos([1 : L]' * pi) * ones(size(x))) .* Y2; 
    
    Y = flipud(Y);
    Y = [Y; Y2];
    %Y
    Y = Y(L + M + 1, :);
    %Y