Wrap NUTS for non-linear fitting features (#635)

* Initial plan

* Add nuts feature to Cargo.toml and pin nuts-rs to 0.17.3

Co-authored-by: hombit <1784493+hombit@users.noreply.github.com>

* Add NUTS algorithm support for non-linear fitting features

Co-authored-by: hombit <1784493+hombit@users.noreply.github.com>

* Enable nuts unconditionally, rename nuts_ndraws to nuts_niter

Co-authored-by: hombit <1784493+hombit@users.noreply.github.com>

* Bump MSRV from 1.85 to 1.88 for nuts-derive compatibility

Co-authored-by: hombit <1784493+hombit@users.noreply.github.com>

* Fix clippy collapsible_if lint in features.rs

Co-authored-by: hombit <1784493+hombit@users.noreply.github.com>

* Resolve merge conflicts with main: bump pyo3 to 0.28.2, pyo3-arrow to 0.17, arrow-* to 58

Co-authored-by: hombit <1784493+hombit@users.noreply.github.com>

* Update CHANGELOG: note comfy-table constraint removal with MSRV bump

Co-authored-by: hombit <1784493+hombit@users.noreply.github.com>

---------

Co-authored-by: copilot-swe-agent[bot] <198982749+Copilot@users.noreply.github.com>
Co-authored-by: hombit <1784493+hombit@users.noreply.github.com>

Author: Copilot

CHANGELOG.md

@@ -9,11 +9,11 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 
 ### Added
 
---
+- `BazinFit`, `LinexpFit`, and `VillarFit` gained `nuts`, `nuts-ceres`, and `nuts-lmsder` algorithms using NUTS (No-U-Turn Sampler) for gradient-based Hamiltonian Monte Carlo optimization. New `nuts_ntune` and `nuts_niter` parameters control the number of tuning and drawing iterations respectively
 
 ### Changed
 
---
+- Minimum supported Rust version (MSRV) bumped from 1.85 to 1.88, required by the `nuts-rs` dependency; `comfy-table` upper-version bound removed (was `<7.2` due to old MSRV)
 
 ### Deprecated