Create mej_eq.py

ejecta mass calculation. follow treatment as in jeffrey '99 and values from stritzinger 2006. Nni from command line is derived using bol_lc.py