Add pipes to all the analysis functions

[stelmo]

I want this to be a reality:

model |> with_changed_bounds(...) |> with_added_reaction(...) |> with_enzyme_constraints(...) |> flux_balance_analysis(...) |> to_dict(:enzymes)

[stelmo]

Related to #747

[codecov]

Codecov Report

Base: 88.79% // Head: 88.62% // Decreases project coverage by -0.17% ⚠️

Coverage data is based on head (773dfda) compared to base (61b2d8c).
Patch coverage: 0.00% of modified lines in pull request are covered.

Additional details and impacted files

@@            Coverage Diff             @@
##             next     #751      +/-   ##
==========================================
- Coverage   88.79%   88.62%   -0.17%     
==========================================
  Files          88       88              
  Lines        2097     2101       +4     
==========================================
  Hits         1862     1862              
- Misses        235      239       +4     

Impacted Files	Coverage Δ
src/analysis/flux_balance_analysis.jl	81.81% <0.00%> (-18.19%)	⬇️
src/analysis/parsimonious_flux_balance_analysis.jl	70.83% <0.00%> (-6.44%)	⬇️

Help us with your feedback. Take ten seconds to tell us how you rate us. Have a feature suggestion? Share it here.

☔ View full report at Codecov.
📢 Do you have feedback about the report comment? Let us know in this issue.

[exaexa]

This is good. Do we add the Result overloads here?

[stelmo]

Yes! I was thinking that each analysis function returns both the input model and the solved jump model (probably wrapped up into a type)