from math import sin

def f(x):
    return sin(x)/x

def integral(a,b,n):
    "Integrálközelítés téglalapszbállyal."
    dx = (b-a)/n

    x = a
    osszeg = 0.0
    for m in range(n):
        osszeg += f(x) * dx
        x += dx
    return osszeg

def test():
    n=5
    print 'A sin x/x  1 és 5 közötti integráljának közelítő értékei\ntéglalap szabály és n felosztás esetén.'
    for i in range(15):
        print 'n= %-8d: %15.13f' % (n, integral(1.0,5.0,n))
        n *= 2

test()