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