fixed bias to linspace(0, M, M+1). plots of different polynomials now corresponds to expected form