Rename function 'poly' to 'polyfromroots'
Refactor code into groups of functions
Refactor function 'polyder'
Extend function 'polyval' to support bivariate and trivariate polynomials