I cannot solve it, even on an EC2 instance with 256 GB of RAM. If all else fails, we can use [theta]=0.I am not sure if it is acceptable, but I will acknowledge the first person to offer the correct solution in my article.