Function generation in Octave from scratch part 4