I can summarize the different possibilities/proposals: (A) Create new properties: add a `P.coef_natural` property, with a suitable documentation ; maybe also add a `P.coef_internal` property. There would be no change to the existing code (only addition of properties). (B) Change `P.coef` attribute into a property, with a suitable documentation. Hide `P.coef` attribute into `P._coef` (change existing code). Do not create more properties (unlike A). - About (A), I don't think that adding `P.coef_natural` would add a risk. - About (B), it may be appreciated that the API does not change (does not occupy more namespace) - Both (A) and (B) would help basic users to get out of the `P.coef` attribute confusion. Side remark (not important):
"natural" coefficients make very little if any sense for some of the other polynomial subclasses, such as Chebyshev -- for those, there's nothing natural about them! Are you sure? Can they not be the weights at different order of approximation of a solution?