File: series_sin.py

package info (click to toggle)
pycode-browser 20120614%2Bgit%2Bb041dd2-7
  • links: PTS
  • area: main
  • in suites: jessie, jessie-kfreebsd
  • size: 6,236 kB
  • ctags: 1,194
  • sloc: python: 2,220; xml: 152; makefile: 39
file content (15 lines) | stat: -rwxr-xr-x 307 bytes parent folder | download | duplicates (12) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
#Example series_sin.py

from pylab import *
from scipy import factorial

x = linspace(-pi, pi, 50)
y = zeros(50)
for n in range(5):
        term = (-1)**(n) * (x**(2*n+1)) / factorial(2*n+1)
        y = y + term
        plot(x,term)

plot(x, y, '+b')
plot(x, sin(x),'xr') # compare with the real one
show()