Finally the code as follows was developped.
This code is based on the interpolation of basic functions with oscillations. I didn’t succeed with the ODEs.
Graphics are obtained as follows:
Not sure if the code is optimized (of course). It is my first CRNN 