sums

import numpy
import pylab
import time

n = pylab.arange(1,10000,1)

J_left = numpy.zeros(n.size)
J_right = numpy.zeros(n.size)
t_run= numpy.zeros(n.size)

for i in range (1,n.size):
    t0=time.time()

    w_left = numpy.arange(0,numpy.pi/2,(numpy.pi/2)/n[i])

    w_right = w_left + (numpy.pi/2)/n[i]

    area_left = numpy.sin(w_left)*(w_right-w_left)

    area_right = numpy.sin(w_right)*(w_right-w_left)

    J_left[i] = area_left
    J_right[i] = area_right

    t1 = time.time()
    t_run[i]=t1-t0

J_labs = numpy.abs(1-J_left)
J_rabs = numpy.abs(1-J_right)

pylab.figure(1)
pylab.plot(n,J_labs)
pylab.show()