add Minimizer and FitResult

Now this is possible:

```
matrix = combined.fit(print_level=-1).save().covariance.to_numpy()
```