Source: https://www.orientation.sutd.edu.sg/fsolve/matlab/examples/tAki13Te-qu.pdf