work on linear newton solver in gsl