哦，明白了，我最初以为是把1000个的数据存在一个文件里面，
如果已经用别的程序产生1000个m文件了，那就该按你的方式来了(1000个文件有点多。。。)
GSL的最新版 gsl-1.13 的windows版本，下载地址：
http://cid-3575dd6acf57d7ac.skydrive.live.com/self.aspx/.Public/mathsoft/gsl-1.13-MinGW.7z
静态库文件是 libgsl.dll.a，libgslcblas.dll.a ，
所以，要用MinGW来编译。
附带有两个例子，使用说明
gsl-1.14出来了
=================================
另外，求解线性代数方程组的最优方法至今还是一个难题
--------------------------------10000个科学难题(数学卷) 第505页有详尽的介绍
python??????????
fortran好像也能求解线性方程组
parigp 应该也能求解！
Sowya 提供了解线性方程组的求解函数 solve().
运行 Sowya.exe,在提示符 >> 后键入矩阵 A 和 b, 然后输入 solve(A*x==b) 回车.
>> A=
input>
--------------------
1 1 -5 -1 1
8 5 -1 0 4
3 -7 -5 -1 1
2 1 8 3 -1
0 6 2 7 -2
--------------------
>> b=
input>
--------------------
9
2
4
-8
15
--------------------
>> solve(A*x==b)
>> 1 1 -5 -1 1 9
8 5 -1 0 4 2
3 -7 -5 -1 1 4
2 1 8 3 -1 -8
0 6 2 7 -2 15
The solution is:
x= E_0 + C_1*E_1 + C_2*E_2 + C_3*E_3 + C_4*E_4 + C_5*E_5
where E_0 is the special solution and the others form the base of the solution:
---------------
E_0 E_1 E_2 E_3 E_4 E_5
------------------------
2509|1670 0 0 0 0 0
1671|1670 0 0 0 0 0
-3853|1670 0 0 0 0 0
118|167 0 0 0 0 0
-68009|15698 0 0 0 0 0
------------------------