lSum Practice with range() -- Sage

a=0 b=5 n=10 h=(b-a)/n f(x)=x^2 sum([f(a+i*h)*h for i in range(n)]).n()

 $\newcommand{\Bold}[1]{\mathbf{#1}}35.6250000000000$

range(n)

 $\newcommand{\Bold}[1]{\mathbf{#1}}\left[0, 1, 2, 3, 4, 5, 6, 7, 8, 9\right]$

[i for i in range(n)]

 $\newcommand{\Bold}[1]{\mathbf{#1}}\left[0, 1, 2, 3, 4, 5, 6, 7, 8, 9\right]$

[i*h for i in range(n)]

 $\newcommand{\Bold}[1]{\mathbf{#1}}\left[0, \frac{1}{2}, 1, \frac{3}{2}, 2, \frac{5}{2}, 3, \frac{7}{2}, 4, \frac{9}{2}\right]$ 
__SAGE__
Traceback (click to the left of this block for traceback)
...
__SAGE__

__SAGE__
Traceback (most recent call last):
  File "c_lib.pyx", line 68, in sage.ext.c_lib.sage_python_check_interrupt (sage/ext/c_lib.c:736)
KeyboardInterrupt
__SAGE__

[a+i*h for i in range(n)]

 $\newcommand{\Bold}[1]{\mathbf{#1}}\left[0, \frac{1}{2}, 1, \frac{3}{2}, 2, \frac{5}{2}, 3, \frac{7}{2}, 4, \frac{9}{2}\right]$

[f(a+i*h) for i in range(n)]

 $\newcommand{\Bold}[1]{\mathbf{#1}}\left[0, \frac{1}{4}, 1, \frac{9}{4}, 4, \frac{25}{4}, 9, \frac{49}{4}, 16, \frac{81}{4}\right]$

[f(a+i*h)*h for i in range(n)]

 $\newcommand{\Bold}[1]{\mathbf{#1}}\left[0, \frac{1}{8}, \frac{1}{2}, \frac{9}{8}, 2, \frac{25}{8}, \frac{9}{2}, \frac{49}{8}, 8, \frac{81}{8}\right]$

sum([f(a+i*h)*h for i in range(n)])

 $\newcommand{\Bold}[1]{\mathbf{#1}}\frac{285}{8}$

sum([f(a+i*h)*h for i in range(n)]).n()

 $\newcommand{\Bold}[1]{\mathbf{#1}}35.6250000000000$

def lSum(a,b,n,f): h=(b-a)/n return sum([f(a+i*h)*h for i in range(n)]).n() f="x**2" show(lSum(0,5,10,lambda x: eval(f))) print "n\tlSum" for n in range(5): print "10^",n,"\t",lSum(0,5,10^n,lambda x: eval(f))

 $\newcommand{\Bold}[1]{\mathbf{#1}}35.6250000000000$ 
n	lSum
10^ 0 	0.000000000000000
10^ 1 	35.6250000000000
10^ 2 	41.0437500000000
10^ 3 	41.6041875000000
10^ 4 	41.6604168750000

n	lSum
10^ 0 	0.000000000000000
10^ 1 	35.6250000000000
10^ 2 	41.0437500000000
10^ 3 	41.6041875000000
10^ 4 	41.6604168750000

@interact def _(a=0, b=5, g = x, n=(0,100,1)): f(x)=g Q=plot(f,(x,a,b)) RP=plot(0,(x,a,b)) h=(b-a)/n for j in range(n): P=plot(f(a+j*h),(x,a+j*h,a+(j+1)*h),fill=True) NP=P+RP RP=NP show(RP+Q,ticks=[[a,a+h..b],None]) k=var('k') html("The Riemann sum approximating $\int_{%s}^{%s} %s \, dx$ adds up to $%s$."%(latex(a),latex(b),latex(g),sum(h*(f(a+k*h)),k,1,n)))

Click to the left again to hide and once more to show the dynamic interactive window

@interact def _(a=1, b=2, g = 4-x^2, n=(1,10,1)): f(x)=g Q=plot(f,(x,a,b)) RP=plot(0,(x,a,b)) h=(b-a)/n for j in range(n): P=plot(f(a+j*h),(x,a+j*h,a+(j+1)*h),fill=True) NP=P+RP RP=NP show(RP+Q,ticks=[[a,a+h..b],None]) k=var('k') html("The Riemann sum approximating $\int_{%s}^{%s} %s \, dx$ adds up to $%s$."%(latex(a),latex(b),latex(g),sum(h*(f(a+k*h)),k,1,n)))

lSum Practice with range()

3329 days ago by MATH4H2012

Click to the left again to hide and once more to show the dynamic interactive window

Click to the left again to hide and once more to show the dynamic interactive window