PE21 solution and supporting functions in Euler library