用於數值計算的開源庫 GNU Scientific Library

一直很好奇使用C/C++語言，該如何進行復雜的數值(數學)計算：例如解方程組，C的標準數學庫就做不到了。

最近剛瞭解到, 這些問題已經有很好用的開源庫了，它就是：
GSL - GNU Scientific Libraryhttp://www.gnu.org/software/gsl/

GSL的功能有: (來自google翻譯)

複數	Roots of Polynomials 多項式的根
Special Functions 特殊功能	Vectors and Matrices 向量和矩陣
Permutations 排列	Sorting 排序
BLAS Support BLAS的支持	Linear Algebra 線性代數
Eigensystems Eigensystems	Fast Fourier Transforms 快速傅立葉變換
Quadrature 正交	Random Numbers 隨機數
Quasi-Random Sequences 擬隨機序列	Random Distributions 隨機分佈
Statistics 統計	Histograms 直方圖
N-Tuples n元組	Monte Carlo Integration 蒙特卡羅積分
Simulated Annealing 模擬退火	Differential Equations 微分方程
Interpolation 插值	Numerical Differentiation 數值微分
Chebyshev Approximation 切比雪夫逼近	Series Acceleration 系列加速度
Discrete Hankel Transforms 離散Hankel變換	Root-Finding 求根
Minimization 最小化	Least-Squares Fitting 最小二乘擬合
Physical Constants 物理常數	IEEE Floating-Point IEEE浮點
Discrete Wavelet Transforms 離散小波變換	Basis splines 基礎樣條

可見GSL的數學功能很強大，例如模擬退火算法，這是佈局佈線必需的算法之一，不知道EDA公司的開發是否也在使用這個庫？

不過，我的問題很簡單，就是解一個一元二次方程組
3x + 2y = 8;
x + 5y = 7;
解： x=2, y=1;
實例代碼請常見附件.

一元二次方程組屬於 Linear Algebra線性代數 範疇，有興趣瞭解 解方程組的過程 可以參考
http://en.wikipedia.org/wiki/System_of_linear_equations
http://en.wikipedia.org/wiki/LU_decomposition