initial mass guess will be 0.001 (maybe add an input parameters...)
removed output