fact.m

求解量子力学的薛定谔方程

function f=fact(x);
%> The file <fact.m> calculates the factorials fact(x) of 
%> a matrix with integer entries, sets infinity for negative values.
%> Call: fact(X)
%> Input: X = vector or matrix, integer values
%> Output: a vector or matrix of same size.
%> Compare the standard gamma function of MATLAB. 
%> We want to avoid all NaN outputs!!!
%>
%>