Double pendulum non-function script for Octave