File: series_sin.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 245 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
from pylab import *
from scipy import *

x = linspace(-pi,pi,40)
a = zeros(40)

plot(x,sin(x))
for n in  range(1,5):
	sign = (-1)**(n+1)
	term = x**(2*n-1) / factorial(2*n-1)
	a = a + sign * term
	print n,sign
	plot(x,term)
plot(x,a,'+')
show()