[C++]求解三元一次方程组

/**
 * 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 ]
*/