利用R和Octave求解线性方程组

方法一、采用R

（1）利用行列式来计算

（2）利用矩阵来计算

方法二、采用Octave

（1）利用行列式来计算

（2）利用矩阵来计算

相比之下，利用octave求解显得更简洁一点。

>> A = [3 -2;2 1];

>> B = [12; 1];

>> X = A\B;

>> printf("x1=%f\nx2=%f\n",X(1),X(2))

x1=2.000000

x2=-3.000000

其实可以编写脚本来计算。

利用R计算：

利用Octave计算：

求解方程组

利用Octave计算：

你看，结果一下就出来了：x = 1, y = 2, z = 3。

利用R计算：

> A <- matrix(c(1,2,3,2,-4,5,3,1,-2),3,3)

> B <- matrix(c(14,-3,7),3,1)

> X <- solve(A,B)

> X

[,1]

[1,] 1

[2,] 2

[3,] 3