My friend and colleague Nic Fillion and I are writing another book, this one on perturbation methods using backward error analysis (and Maple). We have decided to make the supporting materials available by means of Jupyter notebooks with a Maple kernel (there are some Maple worksheets and workbooks already, but going forward we will use Jupyter).
The presentation style is meant to aid reproducibility, and to allow others to solve related problems by changing the scripts as needed.
The first one is up at
https://github.com/rcorless/Perturbation-Methods-in-Maple
Comments very welcome. This particular method is a bit advanced in theory (but it's very simple in practice, for weakly nonlinear oscillators). I haven't coded for efficiency and there may be some improvements possible ("may" he says, sheesh). Comments on that are also welcome.
-r