Comments and Ratings

Made quick test against alpha=beta=gamma=0 3jsymbols in Table 2, in D.M. Brink and R.G. Satchler "Angular Momentum" Oxford Science Publications and can not get a match

Simple. Does what it says and the answers check with Abramowitz and Stegen Table 25.9
using roots(LaguerrePoly(n)) for the
Gauss-Laguerre quadrature abcissas.

Thanks a lot.

sehr gut code!!!

recommended for anyone in early quantum mechanics classes (or teachers looking for material)

The loop in rows 59-64 can be replaced by a simple direct sum:
t = tmin : tmax;
wigner = sum( (-1).^t ./ ( factorial(t) .* factorial(t-t1) .* factorial(t-t2) ...
                .* factorial(t3-t) .* factorial(t4-t) .* factorial(t5-t) ) );

Hi,
I find your function interesting and usefull, but I'm afraid summing the vectors produced reveals an error; while the sum of the coefficients should yield 1 (as Pn(z)=1 for every n), in fact this goes wrong (at least on my computer) at n=25 because of roudof-error (I think), to deviate to something of the order of 1.0e20 at n=99. Do you know of any kind of solution?
Other than that it works very well!
Thanks.

Indeed this function does not work for symbolic computations so should be removed from this list. It works for numerical computations but for this there are better alternatives, e.g. nchoosek or other.

simple algorithm, fast performance.
Why don't I make things like it?

Please to be able to run the IS-95A CDMA Power Control Simulation file, how do get get the required products: CDMA Reference Blockset,Communications Blockset,Communications Toolbox,Signal Processing Blockset,Signal Processing Toolbox,Simulink,Stateflow. Thanks.

Very good routine

Good stuff!
If you need Hermite-Polynomials derived from the Gaussian measure
1/(sqrt(2pi))exp(-x^2/2)
then you just have delete all the '2' in the program and replace them by 1

thank you very much

very good site

there is absolutly not documentation or instruction set. Why load up a set of indeciperable files? Terrible

Nice implementation. This is a good script for series with a small number of terms, however, the transform matrix is highly ill-conditioned. For a series with k terms, cond(M(k)) = O(exp(0.85*k)). i.e. the condition number doubles for each term added to the series. Treating the odd and even terms separately will improve this.

or use Matlab:
F = factor(A) ;
X(:,1) = unique(F(:)') ;
X(:,2) = hist(F,X(:,1)) ;