Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #!/usr/bin/python
- #MrG 2016.0621 2016AB1BC1
- #1c) 49786 l
- import math
- f="2000*math.exp(-x**2/20)"
- a=float(0)
- b=float(8)
- n=10
- def lSum(f,a,b,n):
- h=(b-a)/n
- return sum([f(a+i*h)*h for i in range(n)])
- def rSum(f,a,b,n):
- h=(b-a)/n
- return sum([f(a+i*h)*h for i in range(n)])-f(a)*h+f(b)*h
- print "l\t\tr\t\tt"
- for j in range(4):
- n=10**j
- l=50000-8050+lSum(lambda x: eval(f),a,b,n)
- r=50000-8050+rSum(lambda x: eval(f),a,b,n)
- print l,"\t",r,"\t",(r+l)/2
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement