spherical harmonics finished

1. But very slow. Don't know if it is the algorithm problem or implement
problem. 2. spharm function modified to support negative "m". 3. need to
get deep understand of spherical harmonics to understand and optimize
the algorithm