The code uses two approaches to fit the experimental data with gamma distribution .

1. find a0 + a2P2 + a4P4, then use a0 to search for mini delta, the mixing ratio
2. fit the a0 and delta directly .

#simply make

The settings and data are hard coded.