inclusion of the x^TDx -bx problem with a segmentation fault