File: euler.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 (18 lines) | stat: -rwxr-xr-x 299 bytes parent folder | download | duplicates (6)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
#integrate cosine using Euler method. The result should be sine
from pylab import *

h = 0.01    # stepsize
x = 0.0     # initail values
y = 0.0

ax = []
ay = []

while x < 2*pi:
     y = y + h * math.cos(x)   # Euler method
     x = x + h
     ax.append(x)
     ay.append(y)
plot(ax,ay)
show()