Input x,y points and degree of solution Set all coefficients to 1 For each term from degree down to 0 Compute root quotients Compute high order coefficient Divide high order term from data End for Return array of coefficients