File: series_sc.py

package info (click to toggle)
pycode-browser 1%3A1.03-1
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 9,860 kB
  • sloc: python: 2,919; xml: 152; makefile: 76
file content (20 lines) | stat: -rwxr-xr-x 394 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 
from pylab import *

def factorial(n): # a recursive function
    if n == 0: 
         return 1 
    else: 
         return n * factorial(n-1) 

NP = 100
x = linspace(-pi, pi, NP)
sinx = zeros(NP)
cosx = zeros(NP)

for n in range(10):
        sinx += (-1)**(n) * (x**(2*n+1)) / factorial(2*n+1)
        cosx += (-1)**(n) * (x**(2*n))   / factorial(2*n)

plot(x, sinx)
plot(x, cosx, 'r')
show()