Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #!/usr/bin/python
- from matplotlib import *
- from pylab import *
- a=int(sys.argv[1])
- b=int(sys.argv[2])
- n=int(sys.argv[3])
- f=sys.argv[4]
- def findArea(o,a,b,n,f):
- h=1.0*(b-a)/n
- x=[a+(i+o)*h for i in range(n)]
- y=[f(a+(i+o)*h) for i in range(n)]
- areas=[f(a+(i+o)*h)*h for i in range(n)]
- area=sum(areas)
- return area
- print 'n%10s %10s %10s %10s %10s' %('l', 'r', 'm', 't', 'e')
- for x in range(5):
- l=findArea(0,a,b,n**x,lambda x:eval(f))
- r=findArea(1,a,b,n**x,lambda x:eval(f))
- m=findArea(0.5,a,b,n**x,lambda x:eval(f))
- t=(r+l)/2
- e=(r-l)/2
- print 'n**%2d %10.4f %10.4f %10.4f %10.4f %10.4f' % (x, l, r, m, t, e)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement