使用Sowya解线性方程组
作品展示：平方根近似值计算器 - 数学软件 - 数学研发论坛 - Powered by Discuz! (emath.ac.cn)
下面使用 Sowya 进行计算.
>> A=[3200/3,-1600/3,0,0,0,0,0;
A=[3200/3,-1600/3,0,0,0,0,0;
-1600/3,4200/3,-3200/3,400/3,0,0,0;
0,-3200/3,6400/3,-3200/3,0,0,0;
0,400/3,-3200/3,5600/3,-3200/3,400/3,0;
0,0,0,-3200/3,6400/3,-3200/3,0;
0,0,0,400/3,-3200/3,5600/3,-3200/3;
0,0,0,0,0,-3200/3,6400/3]
A=[3200/3,-1600/3,0,0,0,0,0;
input> [3200/3,-1600/3,0,0,0,0,0;-1600/3,4200/3,-3200/3,400/3,0,0,0;0,-3200/3,6400/3,-3200/3,0,0,0;0,400/3,-3200/3,5600/3,-3200/3,400/3,0;0,0,0,-3200/3,6400/3,-3200/3,0;0,0,0,400/3,-3200/3,5600/3,-3200/3;0,0,0,0,0,-3200/3,6400/3]
--------------------
3200/3 -1600/3 0 0 0 0 0
-1600/3 4200/3 -3200/3 400/3 0 0 0
0 -3200/3 6400/3 -3200/3 0 0 0
0 400/3 -3200/3 5600/3 -3200/3 400/3 0
0 0 0 -3200/3 6400/3 -3200/3 0
0 0 0 400/3 -3200/3 5600/3 -3200/3
0 0 0 0 0 -3200/3 6400/3
--------------------
>> b=[-184/3;-89/3;-12/3;-6/3;-12/3;-152;-12/3]
input> [-184/3;-89/3;-12/3;-6/3;-12/3;-152;-12/3]
--------------------
-184/3
-89/3
-12/3
-6/3
-12/3
-152
-12/3
--------------------
>> solve(A*x==b)
>> 3200/3 -1600/3 0 0 0 0 0 -184/3
-1600/3 4200/3 -3200/3 400/3 0 0 0 -89/3
0 -3200/3 6400/3 -3200/3 0 0 0 -12/3
0 400/3 -3200/3 5600/3 -3200/3 400/3 0 -6/3
0 0 0 -3200/3 6400/3 -3200/3 0 -12/3
0 0 0 400/3 -3200/3 5600/3 -3200/3 -152
0 0 0 0 0 -3200/3 6400/3 -12/3
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 + C_6*E_6 + C_7*E_7
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 E_6 E_7
------------------------
-47|200 0 0 0 0 0 0 0
-71|200 0 0 0 0 0 0 0
-353|960 0 0 0 0 0 0 0
-113|300 0 0 0 0 0 0 0
-611|1600 0 0 0 0 0 0 0
-23|60 0 0 0 0 0 0 0
-575980|2976000 0 0 0 0 0 0 0
------------------------
>>