code: function QuadSD_ak1(NN, u, v, p, q, iPar){
                     // Divides p by the quadratic 1, u, v placing the quotient in q and the remainder in a, b
                     // iPar is a dummy variable for passing in the two parameters--a and b--by reference
                     q[0] = iPar.b = p[0];
                     q[1] = iPar.a = -(u*iPar.b) + p[1];
  
                     for (var i = 2; i < NN; i++){
                         q[i] = -(u*iPar.a + v*iPar.b) + p[i];
                         iPar.b = iPar.a;
                         iPar.a = q[i];
                     } 
                     return;
                 }
docs: Find the roots --> Begin Jenkins-Traub /*
  A verbatim copy of Mr. David Binner's Jenkins-Traub port