package matriz;

public class MatrizApp {
    public static void main(String[] args) {
//vectores y matrices
        double[][] a1={{1, 2, 3},{4,5,6},{7,8,9}};
        Matriz a=new Matriz(a1);
        double[] v1={1,2,3};
        Vector v=new Vector(v1);
        Vector r=Matriz.producto(v, a);
        System.out.println("Vector "+r);
        r=Matriz.producto(a, v);
        System.out.println("Vector "+r);
//matrices
        double[][] b1={{1, 0, -1},{2,1,3},{-1, 0, 2}};
        Matriz b=new Matriz(b1);
        Matriz re=Matriz.suma(a, b);
        System.out.println("matriz "+re);
        re=Matriz.producto(a, b);
        System.out.println("matriz "+re);
//número por una matriz
        re=Matriz.producto(a, 2.0);
        System.out.println("matriz "+re);
        re=Matriz.producto(2.0, a);
        System.out.println("matriz "+re);
//determinantes
        System.out.println("determinante "+b.determinante());
        System.out.println("determinante "+a.determinante());
        double[][] h1={{3, 1, -1, 2, 1},{-2, 3, 1, 4, 3},{1, 4, 2, 3, 1},{5, -2, -3, 5, -1},{-1, 1, 2, 3, 2}};
        Matriz h=new Matriz(h1);
        System.out.println("determinante "+h.determinante());
//matriz inversa
        double[][] c1={{1.8, -3.8, 0.7, -3.7},{0.7, 2.1, -2.6, -2.8},{7.3, 8.1, 1.7, -4.9},{1.9, -4.3, -4.9, -4.7}};
        Matriz c=new Matriz(c1);
        Matriz in=Matriz.inversa(c);
        System.out.println("matriz inversa "+in);
        System.out.println("matriz unidad "+Matriz.producto(c, in));
//sistema de ecuaciones lineales
        double[][] m1={{3, -1, 0}, {-2, 1, 1}, {2, -1, 4}};
        Matriz coef=new Matriz(m1);
        double[] n1={5, 0, 15};
        Vector ter=new Vector(n1);
        Vector solucion=Matriz.producto(Matriz.inversa(coef), ter);
        System.out.println("solución "+solucion);
//otro sistema de ecuaciones
        double[][] m2={{7.9, 5.6, 5.7, -7.2}, {8.5, -4.8, 0.8, 3.5}, {4.3, 4.2, -3.2, 9.3}, {3.2, -1.4, -8.9, 3.3}};
        coef=new Matriz(m2);
        double[] n2={6.68, 9.95, 8.6, 1};
        ter=new Vector(n2);
        solucion=Matriz.producto(Matriz.inversa(coef), ter);
        System.out.println("solución "+solucion);
//polinomio característico
        double[][] p1={{1,2,3,4},{2,1,2,3},{3,2,1,2},{4,3,2,1}};
        Matriz p=new Matriz(p1);
        double pol[]=p.polCaracteristico();
        System.out.println("Polinomio característico");
        for(int i=0; i<pol.length; i++){
            System.out.print((double)Math.round(pol[i]*1000)/1000+" , ");
        }
        System.out.println("");
//valores propios y vectores propios
        System.out.println("Valores y vectores propios ");
        double[][] val1={{7, -1, -1}, {-1, 5, 1},{-1, 1, 5}};
        Matriz matriz=new Matriz(val1);
        Matriz vectores=new Matriz(matriz.n);
        double[] valores=new double[matriz.n];
        try{
            vectores=matriz.valoresPropios(valores, 20);

        }catch(ValoresExcepcion ex){
            System.out.println("Al calcular los valores propios se ha producido una excepción\n "
                +ex.getClass()+ " con el mensaje "+ ex.getMessage());
        }
        for(int i=0; i<matriz.n; i++){            //valores propios
        	System.out.print(valores[i]+" , ");
        }
        System.out.println("");
        System.out.println(vectores);       //vectores propios

        try  {
//espera la pulsación de una tecla y luego RETORNO
            System.in.read();
        }catch (Exception e) {  }
    }
}