/*
	Cliente
*/

import Vector , Matriz

program Cliente    |[

/*************
************** Casos de Prueba de Vector3
*************/

   var V1 := ZeroV(3) : Vector;    // Zero vector of dim = 3
   printVector("V1 = ", V1);

   V1 := MakeVector3(1.0, 2.0, 1.0);     // vector of dim = 3
   printVector("V1 = ", V1);
   
   var V2 := MakeVector3(-1.0, 2.0, 3.0) : Vector;
   printVector("V2 = ", V2);
   
   var V3 := SumaV(V1, V2): Vector;
   printVector("V3 = ", V3);
   
//   var V4 := MultEscV(1.0/Modulo(RestaV(V1, V2)), RestaV(V1, V2)): Vector;
   var V4 := NormalizeV(RestaV(V1, V2)): Vector;
   printVector("V4 = ", V4);
   
   write(Modulo(V4));
   
   write ("V4 . V5 = " + ProdEsc(V4,V3));
   
   printVector("V4 x V4 = ", ProdVect(V4,V4));
   
/*************
************** Casos de Prueba de Matriz3x3
*************/

	var M1 := MakeMatriz3x3S( 1.0, 0.0, 0.0,
							  2.0, 4.0, 0.0,
							  3.0, 2.0, 9.0 ) : Matriz;
							
	var M2 := MakeMatriz3x3V( MakeVector3(0.0, 1.0, 2.0),
							  MakeVector3(1.0, 2.0, 3.0),
							  MakeVector3(2.0, 3.0, 4.0) ) : Matriz;

	printMatriz("M1 = ", M1);
	printMatriz("M2 = ", M2);

	printMatriz("Ident = ", Ident(3) );     // Matriz Identidad
	printMatriz("ZeroM = ", ZeroM(3) );     // Matriz Ceros

	var A := MakeMatrizM( 3, %[ i,j : 3,3 : toDouble(i+j-1) ]% ) : Matriz;
	printMatriz("A = ", A);
	
	write("det(A) = "     + Determinante( A ) );              // det(A)
	
	V1 := MakeVector3(1.0, 2.0, 1.0);     // vector of dim = 3
	printVector("V1 = ", V1);
	
	var V5:= MultMV(A, V1) : Vector;    // V5 = A * V1
	printVector("V5 = ",  V5);                 
	
	var B := MakeMatrizM( 3, %[ i,j : 3,3 : (i < j) => 0.0 | 1.0 ]% ) : Matriz; 
	printMatriz("B = ", B);
	
	B := RestaM (MultEscM(0.5,A), Ident(3));
	printMatriz("B = ", B);
	write("det(B) = "     + Determinante( B ) );              // det(B)
	
	B :=  MultEscC( -1.0, 2, MultEscF( -2.0, 1, B ));
	printMatriz("B = ", B);
	write("det(B) = "     + Determinante( B ) );              // det(B)

	var C:= SumaM( MultEscM(2.0, A), B) : Matriz;
	printMatriz("C = ", C);
	write("det(C) = "     + Determinante( C ) );              // det(B)
	
	write("A = C? " + IgualM(A, C));
	
	write("es cero? " + IgualV (MultMV(A, V5), ZeroV(3)));
	
	write("es cero? " + IgualM (MultEscM(0.0, A), ZeroM(3)));
	
	printMatriz("(Cˆ2)' = ", Transpuesta(MultM(C,C)));
	
	var D:= MultEscF(1.0/Determinante(C), 2, C) : Matriz;
	printMatriz(" Norm = ", MultM(D, D));
	write("det(D) = "     + Determinante( D ) );              // det(D)
	
	
// A la matriz A multiplicada por 0.5 restarle la matriz Identidad y dar resultado en B
//  Modificar la matriz B multiplicando la 2da fila por -2 y la 3ra columna por -1
//  Multiplicar la matriz A por 2 y sumar a B, obteniendo C
//  Determinar si A y C son iguales
//  Determinar si A multiplicada a V5 es el vector cero.
//  Multiplicar A por cero y verficar si es igual a la matriz ZeroM
//  Obtener la transpuesta de la matriz C2

   
/*************
************** FIN
*************/

   skip
   
]|