Here's a version of solving inverse kinematics using CCD (Cylic Coordinate Descent). The prevoius methods were analytical methods, this one is iterative, and starts moving into solving for an N number of bones.

Here's the p5Js sketch:
https://editor.p5js.org/admiralra/sketches/chZFS9uma