/**

 * author:johnny zen

 * date:2017-09-20 11:19

 * function:Calculate Ternary system of equations

 * notice:时间仓促，仅仅实现功能，方便使用，代码质量不可参考！！！

 */

#include<iostream>

using namespace std;

template<class T>

void input(T matrix[][]){

    cout<<"please input matrix element's data"<<endl;

    for(int i = ;i<;i++){

        for(int j=;j<;j++){

            cin>>matrix[i][j];

        }

    }

    cout<<"input ok";

}

template<class T>

void calc(T matrix[][]){

    T     base_D = matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][];//计算行列式

        base_D = base_D-(matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][]);

    if(base_D != ){

        T     x_D = matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][];

            x_D = x_D-(matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][]);

        T     y_D = matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][];

            y_D = y_D-(matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][]);

        T     z_D = matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][];

            z_D = z_D-(matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][] + matrix[][]*matrix[][]*matrix[][]);

        T x =  x_D/base_D;

        T y =  y_D/base_D;

        T z =  z_D/base_D;

        cout<<"[ x:"<<x<<"; y:"<<y<<"; z:"<<z<<" ]"<<endl;

    }else{

        cout<<"【无解】";

//        return DBL_MIN;

    }

}

//计算原理：行列式

int main(){

    double matrix[][] ; //三元一次方程组

    input<double>(matrix);

    calc<double>(matrix);

    system("pause");

    return ;

}

/*
demo

2x-y+z=10;

3x+2y-z=16;

x+6y-z=28;

2 -1 1 10

3 2 -1 16

1 6 -1 28

output:input ok[ x:4.18182; y:5.09091; z:6.72727 ]

*/