fact.m

Three dimensional hydrogen atom spherical harmonics visualization.

function y = fact(x);

    y = zeros(size(x));
    
    for j = 1 : length(x)
        y(j) = prod(1 : x(j));
    end
    
    y(find(x < 0)) = inf;