I've tried a few different ways to get sympy to integrate a convolution of a pulse from -0.5 to 0.5 with itself, but everything I try produces exceptions or what looks like incorrect results. I've got a gist showing a couple of my attempts at
http://nbviewer.ipython.org/gist/mwiebe/48bed4e5d52335c47e25 . I worked this out with pen+paper, but wanted to show how nicely one could use sympy to derive it, only to find that everything I try fails. Any suggestions for this, or is this something I shouldn't expect sympy to be able to do?
Thanks,
Mark