using GSL for root find and special functions