i get zero as well. although that i dont believe that is correct. see
below.
In [884]: x,a,m = symbols('x a m')
In [885]: pprint(Integral(sin(pi*x)*sin(pi*x),(x,0,1)).doit())
1/2
In [886]: pprint(Integral(sin(pi*x)*sin(m*pi*x),(x,0,1)).doit())
2
π⋅sin(π⋅m) π⋅m ⋅sin(π⋅m)
- ────────────────────── + ──────────────────────
2 2 2 2 4 2 2 2 2 4
- π + 2⋅π ⋅m - π ⋅m - π + 2⋅π ⋅m - π ⋅m
In [887]: pprint(Integral(sin(pi*x)*sin(m*pi*x),(x,
0,1)).doit().limit(m,1))
1/2
In [888]: m = Symbol('m', integer = True)
In [889]: pprint(Integral(sin(pi*x)*sin(m*pi*x),(x,0,1)).doit())
0
i think this is inconsistent.
--
You received this message because you are subscribed to the Google Groups
"sympy" group.
To post to this group, send email to sy...@googlegroups.com.
To unsubscribe from this group, send email to
sympy+unsubscr...@googlegroups.com.
For more options, visit this group at
http://groups.google.com/group/sympy?hl=en.
