Added proper code for calculating theoretical map of f_red.
Also added plotting to fit to data using simple grid search.