User interface of the `polynomial_model()` function

[choisy]

Let's make the user interface and drop the type and the link arguments. And let's allow k the possibility to take a range of values. In the help of the function you can still cite the names of the models.