File: bes2.py

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

def vjn(n,xarray):
	y = []
	for  x in xarray:
		val = special.jn(n,x)        # Compute Jn(x)
		y.append(val)
	return y

a = linspace(0,10,100)
for n in range(5):
	b = vjn(n,a)
	plot(a,b)
show()