test_fact.m

Three dimensional hydrogen atom spherical harmonics visualization.

x = [-3 0 5]

fact(x)

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