File: diff1.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 (14 lines) | stat: -rwxr-xr-x 482 bytes parent folder | download | duplicates (8) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
from pylab import *
dx = 0.1                # value of x increment
x = arange(0,10, dx)
y = sin(x)
yprime = []             # empty list
for k in range(99):     # from 100 points we get 99 diference values
        dy = y[k+1]-y[k]
        yprime.append(dy/dx)

x1 = x[:-1]       	# A new array without the last element
x1 = x1 + dx/2          # The derivative corresponds to the middle point
plot(x1, yprime, '+')
plot(x1, cos(x1))       # Cross check with the analitic value
show()