Added options for numerical integrator

Options flag "iopt" for choosing numerical integrator