use analytic derivative for marginal effects
add creating_settings.m and run.m