File: sph1.py

package info (click to toggle)
pycode-browser 1%3A1.02%2Bgit20181006-3
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 6,088 kB
  • sloc: python: 2,779; xml: 152; makefile: 71
file content (15 lines) | stat: -rwxr-xr-x 439 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
from pylab import *
from scipy import special

a_th = []	# list to store polar angle theta from -90 to + 90 deg
a_sph = []	# list to store absolute values if sperical harminics
phi = 0.0	# Fix azimuth, phi at zero

theta = -pi/2    	# start theta from -90 deg
while theta < pi/2:
	h = special.sph_harm(0,10, phi, theta)	# (m, l , phi, theta)
	a_sph.append(abs(h))
	a_th.append(theta * 180/pi)
	theta = theta + 0.02
plot(a_th,a_sph)
show()