Analyze this pseudocode

thetas=Math.acos(dotvec(deltavec, j.deltavec)/(magvec(deltavec)*magvec(j.deltavec)))

axis = scalevec(thetas, unitvec(crossvec(subvec(deltavec, j.deltavec), deltavec)))

axissum = newvec(0,0,0)
axissum = addvec(axissum, axis)
...
myaxis = scalevec(1/axis.length, axissum)

goaldeltavec = multmat(rotationmatrixbyaxis(myaxis), j.deltavec)

thetas[i]=Math.acos(dotvec(deltavec[i], j.deltavec[i])/(magvec(deltavec[i])*magvec(j.deltavec[i])))

axis[i] = scalevec(thetas[i], unitvec(crossvec(subvec(deltavec[i], j.deltavec[i]), deltavec[i])))

axissum = newvec(0,0,0)
axissum = addvec(axissum, axis[i])
...
myaxis = scalevec(1/axis.length, axissum)

goaldeltavec[i] = multmat(rotationmatrixbyaxis(myaxis), j.deltavec[i])

Maybe it's wrong, because you can't just sum rotation axes.

Analyze my anus

ChatGPT solved this shit.
https://chatgpt.com/share/68b29743-11d0-8007-8603-95ba45665906