I guess that the problem is fixed, but the matrices become so ill-conditioned that they are singular to working precision
tried rcond, doesn't seem to work. made minor changes to calculate_newton_basis, displays how many points used to get power fcn below eps*10