Implement 1D `Profile` class with orthogonal polynomials like Jacobi or Chebyshev

[dpanici]

For fitting data to a profile, to give an option other than spline (which can be prone to oscillations) or PowerSeriesProfile (which cannot be used past a certain degree for fitting due to numerics), an option to use an orthogonal polynomial would be useful. Jacobi is probably the easiest given we already have derivatives implemented for them (unlike Chebyshev)

[unalmis]

chebyshev would be better done with the DCTs on #1119