Added routine for implicit euler
Including use of GSL libraries