Cannot produce simple Talbot carpet

[matt8s]

Cannot produce simple Talbot carpet

[phoebe-p]

This is very strange. I downloaded your code and tried to run it and had the same issue, but if I change the number of Fourier orders from 100 to 9 I no longer get NaN (although I don't know if the results are reasonable). I realise this does not solve the issue but perhaps it gives a starting point for further investigation...

P.S. Thank you very much for making me aware of Inkstone, I did not know about it before

[matt8s]

Thank you for your prompt response. Yes, the code does run for me with up to 12 Fourier orders, but not for 13 or higher. The code seems to converge with this number of orders, so while it would be nice to enable arbitrarily large numbers of Fourier orders, I think this is good enough for my purposes!