My only suggestion here is to check the expressions to see if they are correct, either by working them out by hand, or checking against some other CAS or WolframAlpha. That will tell you if there is a problem in your formulation or a bug in SymPy.
Aaron Meurer On Thu, Jul 16, 2015 at 10:19 AM, <[email protected]> wrote: > Hello Aaron, > > Thank you so much for the help. I only dig out SymPy about once a year so > am not as familiar with it as I would like. > > I have attached updated code with alterations that you suggested. The > formulation is slightly different (and simpler). The values I get for odd > values of N look right now - they converge in a sensible fashion and match > what I expect for this problem (You can see as N increases they tend > towards ~0.1). However the even values of N do not look right - they > oscillate about 0 as I increase the number of terms in the series. They > should all be positive and follow the even values. Adding extra precision > with evalf() did not seem to help. > > I guess at this point it's actually to do with the properties of the sinc > integral, but thought I should check that I'm not doing something stupid > (again)! > > Thanks again, > Toby > > > -- > You received this message because you are subscribed to the Google Groups > "sympy" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to [email protected]. > To post to this group, send email to [email protected]. > Visit this group at http://groups.google.com/group/sympy. > To view this discussion on the web visit > https://groups.google.com/d/msgid/sympy/2e5b13a0-bbf9-4334-80b5-6d6e4ca8b56c%40googlegroups.com > <https://groups.google.com/d/msgid/sympy/2e5b13a0-bbf9-4334-80b5-6d6e4ca8b56c%40googlegroups.com?utm_medium=email&utm_source=footer> > . > > For more options, visit https://groups.google.com/d/optout. > -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To post to this group, send email to [email protected]. Visit this group at http://groups.google.com/group/sympy. To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/CAKgW%3D6%2B1BYuRN-p20j0_970Zzantyatmm_yv%3D7scmEdU7Jzo3g%40mail.gmail.com. For more options, visit https://groups.google.com/d/optout.
