Skip to main content
6 events
when toggle format what by license comment
Sep 23, 2015 at 9:17 vote accept Funzies
Sep 22, 2015 at 9:22 comment added Funzies @MarcoB Thanks for your reply. Indeed, your code gives a good speed-up. However, to ensure convergence I need to take a lot more terms into account. Going to 100 terms, the results becomes 'weird', probably due to numerical errors as you mentioned. Is there any way this can be fixed?
Sep 21, 2015 at 16:27 comment added MarcoB @march I noticed that as well, while playing with the machine-precision version. I amended my answer accordingly, also including your suggestion to prevent needless recalculation of the intensity expression.
Sep 21, 2015 at 16:25 history edited MarcoB CC BY-SA 3.0
Added plotting
Sep 21, 2015 at 16:17 comment added march Note the most recent comment by the OP in the OP. There are also some numerical stability issues it seems.
Sep 21, 2015 at 16:14 history answered MarcoB CC BY-SA 3.0