File: trapez.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 (17 lines) | stat: -rwxr-xr-x 376 bytes parent folder | download | duplicates (3) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
from math import *

def sqr(a):
        return sqrt(1.0 - a**2)

def trapez(f, a, b, n):
        h = (b-a) / n
        sum = f(a)
        for i in range (1,n):
                sum = sum + 2 * f(a + h * i)
        sum = sum + f(b)
        return 0.5 * h * sum

print( 4*trapez(sqr,0.,1.,100))
print( 4*trapez(sqr,0.,1.,1000))
print( 4*trapez(sqr,0,1,100))   # Why the error ?