征集求解线性方程组的数学软件

wayne

哦，明白了，我最初以为是把1000个的数据存在一个文件里面，
如果已经用别的程序产生1000个m文件了，那就该按你的方式来了(1000个文件有点多。。。)

wayne

GSL的最新版 gsl-1.13 的windows版本，下载地址：
http://cid-3575dd6acf57d7ac.skyd ... t/gsl-1.13-MinGW.7z

静态库文件是 libgsl.dll.a，libgslcblas.dll.a ，
所以，要用MinGW来编译。
附带有两个例子，使用说明

wayne

gsl-1.14出来了


=================================
另外，求解线性代数方程组的最优方法至今还是一个难题
--------------------------------10000个科学难题(数学卷)   第505页有详尽的介绍

mathematica

python??????????

mathematica

fortran好像也能求解线性方程组

mathematica

parigp 应该也能求解！

haifeng

Sowya 提供了解线性方程组的求解函数 solve().  

运行 Sowya.exe,  在提示符 >> 后键入矩阵 A 和 b, 然后输入 solve(A*x==b) 回车.


>> A=[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]
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]
--------------------

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=[9;2;4;-8;15]
input> [9;2;4;-8;15]
--------------------

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

------------------------