java解四元一次方程

问题：求解四元一次方程，例如

       x+y-z=9;
       x-y+2z=8;
       2x+y-z=7;
       x+y+z=m;
这个解是 x=-2.0 y=32.0 z=21.0 m=51.0
 
如果你想利用迭代暴力计算的话，不用试了，我替你试过了，太慢，慢的一天你也解不出来，比如x，y，z，m的范围是从-2.0到51.0
 
那么解决不了了吗？这时脑海中依稀想起来了大学学过线性代数，于是乎找美丽的度娘亲热了大半个下午，终于把线性代数这部分拾起来一部分o(╯□╰)o

顺便也帮大家温习下，先来个入门简单的

PS：或许好多客官认为这个温习没必要，傻逼才看呢，可是我想说，等你想拓展这个程序，想要解5元一次方程，n元一次方程的时候，你绝对用得着。但是如果你精通线代的话，算我没说。。。

然后再看下克拉默法则：

 
然后再看看代数余子式：

 
汗，最后步入正题，代码就简单了，其实就是套公式，看代码吧。。。
（注意要解的方程组把变量都放到左边，常亮放到右边，然后再写strss这个二维数组）
public
class ArithmeticOfCalcFour {
   
    public
static void main(String[] args) {
//     例子
//     x+y-z=9;
//     x-y+2z=8;
//     2x+y-z=7;
//     x+y+z=m;
//     解x=-2,y=32,z=21,m=51
       double[][] strss={
              {1,1,-1,0,9},
            {1,-1,2,0,8},
            {2,1,-1,0,7},
            {1,1,1,-1,0},
                      };
       calculate(strss);
       //test();
    }
   
    public
static void calculate(double[][] strss) {
      
       double[][] dss={
                         {strss[0][0],strss[0][1],strss[0][2],strss[0][3]},
                         {strss[1][0],strss[1][1],strss[1][2],strss[1][3]},
                         {strss[2][0],strss[2][1],strss[2][2],strss[2][3]},
                         {strss[3][0],strss[3][1],strss[3][2],strss[3][3]}
                     };
       double[][] d1ss={
                         {strss[0][4],strss[0][1],strss[0][2],strss[0][3]},
                         {strss[1][4],strss[1][1],strss[1][2],strss[1][3]},
                         {strss[2][4],strss[2][1],strss[2][2],strss[2][3]},
                         {strss[3][4],strss[3][1],strss[3][2],strss[3][3]}
                     };
       double[][] d2ss={
                         {strss[0][0],strss[0][4],strss[0][2],strss[0][3]},
                         {strss[1][0],strss[1][4],strss[1][2],strss[1][3]},
                         {strss[2][0],strss[2][4],strss[2][2],strss[2][3]},
                         {strss[3][0],strss[3][4],strss[3][2],strss[3][3]}
                     };
       double[][] d3ss={
                         {strss[0][0],strss[0][1],strss[0][4],strss[0][3]},
                         {strss[1][0],strss[1][1],strss[1][4],strss[1][3]},
                         {strss[2][0],strss[2][1],strss[2][4],strss[2][3]},
                         {strss[3][0],strss[3][1],strss[3][4],strss[3][3]}
                     };
       double[][] d4ss={
                         {strss[0][0],strss[0][1],strss[0][2],strss[0][4]},
                         {strss[1][0],strss[1][1],strss[1][2],strss[1][4]},
                         {strss[2][0],strss[2][1],strss[2][2],strss[2][4]},
                         {strss[3][0],strss[3][1],strss[3][2],strss[3][4]}
                     };
      
       double d=calculateFour(dss);
       if(d==0){
           System.out.println("无解或多个解");
           return ;
       }
       double d1=calculateFour(d1ss);
       double d2=calculateFour(d2ss);
       double d3=calculateFour(d3ss);
       double d4=calculateFour(d4ss);
      
       double x=d1/d;
       double y=d2/d;
       double z=d3/d;
       double m=d4/d;
       System.out.println("结果：x="+x+"y="+y+"z="+z+"m="+m);
       System.out.println("计算完毕!");
    }
   
    /**
     * 计算4阶线性方程组
     * @param strss
     * @return
     */
    private
static double calculateFour(double[][] strss) {
      
//     {strss[0][0],strss[0][1],strss[0][2],strss[0][3]},
//     {strss[1][0],strss[1][1],strss[1][2],strss[1][3]},
//     {strss[2][0],strss[2][1],strss[2][2],strss[2][3]},
//     {strss[3][0],strss[3][1],strss[3][2],strss[3][3]}
      
       double[][] t00={
                         {strss[1][1],strss[1][2],strss[1][3]},
                         {strss[2][1],strss[2][2],strss[2][3]},
                         {strss[3][1],strss[3][2],strss[3][3]}
                     };
       double[][] t01={
                         {strss[1][0],strss[1][2],strss[1][3]},
                         {strss[2][0],strss[2][2],strss[2][3]},
                         {strss[3][0],strss[3][2],strss[3][3]}
                     };
       double[][] t02={
                         {strss[1][0],strss[1][1],strss[1][3]},
                         {strss[2][0],strss[2][1],strss[2][3]},
                         {strss[3][0],strss[3][1],strss[3][3]}
                     };
       double[][] t03={
                         {strss[1][0],strss[1][1],strss[1][2]},
                         {strss[2][0],strss[2][1],strss[2][2]},
                         {strss[3][0],strss[3][1],strss[3][2]}
                     };
       double result=Math.pow(-1, 0+0)*strss[0][0]*calculateThree(t00)+ Math.pow(-1, 0+1)*strss[0][1]*calculateThree(t01) +
              Math.pow(-1, 0+2)*strss[0][2]*calculateThree(t02)+ Math.pow(-1, 0+3)*strss[0][3]*calculateThree(t03) ;
      
       return result;
    }
 
    /**
     * 计算3阶行列式
     * @param strss
     * @return
     */
    private
static double calculateThree(double[][] strss) {
//     {strss[0][0],strss[0][1],strss[0][2]
//     {strss[1][0],strss[1][1],strss[1][2]
//     {strss[2][0],strss[2][1],strss[2][2]
      
       double result=strss[0][0]*strss[1][1]*strss[2][2] +  strss[0][1]*strss[1][2]*strss[2][0] +  strss[0][2]*strss[1][0]*strss[2][1] + 
              -strss[0][2]*strss[1][1]*strss[2][0] -strss[0][1]*strss[1][0]*strss[2][2] - strss[0][0]*strss[2][1]*strss[1][2] ;
      
       return result;
    }