4D Poincare surface of sections

[Shahzadi]

TL;DR Summary

I'm trying to plot a 4D Poincaré surface of sections for a system with 3 degrees of freedom. I have written a code for this. The code provides results for the integration of motion, but it does not generate the data points needed to plot the Poincaré surface of section. However, the same code works very well in generating data points to plot the Poincaré surface of section when I change the dynamical system to 2 degrees of freedom.

Can anyone help me with this? I have posted my code at the below:

https://mathematica.stackexchange.com/questions/306820/4d-poincare-surface-of-sections