Create main.m

The main program for flux estimation