Documenting my first attempt to solve the program: procedure calculation
Derived formula for 126
Added problem 126