How-to: sum

 1 Suppose I have defined a function f(x)=x^2  and a list x=[0,1,2,3,4]  Now I can calculate the product x[2]*f(x[2])  But: If I use the command k=var('k') sum(x[k]*f(x[k]),k,0,4)  I will get the message: TypeError: unable to convert x (=k) to an integer  How can I realize such a summation? asked Mar 26 '11 tmaxara 65 ● 2 ● 4 ● 8 kcrisman 6614 ● 13 ● 66 ● 149

 8 Here are some ways to sum the squares of the numbers from 0 to 1000, with timings: sage: var('x');f = x^2; sage: timeit('sum(f,x,0,1000)') 5 loops, best of 3: 271 ms per loop sage: timeit('sum([f(x=t) for t in [0..1000]])') 5 loops, best of 3: 39 ms per loop sage: timeit('sum([t^2 for t in [0..1000]])') 625 loops, best of 3: 1.21 ms per loop sage: timeit('reduce(lambda x,y:x+y^2,[0..1000],0)') 625 loops, best of 3: 1.13 ms per loop sage: timeit('reduce(lambda x,y:x+y^2,xrange(0,1001),0)') 625 loops, best of 3: 945 µs per loop  So if speed matters to you, you should learn a little python! For example, the syntax [EXPR for VAR in LIST] is very handy - it evaluates EXPR with VAR ranging over the list (or any iterable object) LIST. You can even use it with a filtering: sum([x^2 for x in [0..1000] if is_prime(x)])  will sum the squares of primes in this range. To understand the last two summation examples above you should read about lambda, reduce, and xrange in a python tutorial. It will pay off! posted Apr 02 '11 parzan 848 ● 3 ● 12 ● 30
 4 Given the syntax of the sum() function, which can be read by executing "sum?" (with the question mark) in the Sage notebook, you don't need to define a list to compute the sum. For example, sage: var('x') sage: f = x^2 # note that you don't write f(x) = x^2 sage: sum(f,x,0,4) 30  Also, be careful defining a variable, x, as well as a list with the same name. (You wrote f = x^2 as well as x=[0,1,2,3,4].) Hope this helps. posted Mar 26 '11 cswiercz 809 ● 5 ● 15 ● 33 http://www.cswiercz.info/
 2 The first syntax is symbolic summation: sage: var('a,b,i') sage: sum(2^i, i, a, b) -2^a + 2^(b + 1)  It will always try to evaluate the sum symbolically. If you use something that doesn't understand symbolic variables (like GF(8), which is implemented by polybori) then you get the "unable to convert to integer" error. Because polybori can only exponentiate by integers, and not by formal variables. The second syntax is Python's list creation and summation of lists: sage: [ 2^j for j in range(0,3) ] [1, 2, 4] sage: sum(_) 7  Note that the variable inside the list comprehension (which I called j here) is automatically created and will overwrite other variables. In your second example, calling the inner variable i overwrites the previously-defined symbolic variable var('i') posted May 16 '11 Volker Braun 2238 ● 5 ● 22 ● 52 This is a great thread, if one can call this a thread. Many different points of view and followups! kcrisman (May 17 '11)
 0 To add to the diversity of answers here, another valid approach is to add a single line to what you already have: f(x)=x^2 g(x)=x*f(x) k=var('k') sum(g(k),k,0,4)  This will return 100. The point is that you can use this syntax to sum the values of a proper function, but not arbitrary Python expressions. posted May 18 '11 Jason Bandlow 381 ● 1 ● 8 ● 18

[hide preview]