#! /bin/octave -qf
#
#	Autor: Giana Francisco 
#	Curso: 3R4 - Universidad Tecnologica Nacional
#
#
Titulo="\n\n\n\n\n\n\n\n\n\n\n\n\n#\n#\tResolucion de circuitos por metodo de las mallas\n#\n#\tAutor: Giana Francisco\n#\tCurso: 3R4 - Universidad Tecnologica Nacional\n#\n#\n\n\n\n\n\n\n\n\n\n\n\n\n";
printf ("%s",Titulo);
N=input ("Ingrese el tamano de la matriz NxN de impedancias: N = ");
Impedancias=rand(N,N);
for k = 1:N
       for l = k:N
               printf ("Ingrese el termino Z%i%i: ",k,l);
               Impedancias(k,l)=input (" ");
               Impedancias(l,k)=Impedancias(k,l);
       endfor
endfor
printf ("%s", "La matriz de impedancias es:\n");
Impedancias
Tensiones=rand(N,1);
for k = 1:N
	printf ("Ingrese la tension V%i: ",k);
        Tensiones(k,1)=input (" ");
endfor
printf ("%s", "el vector de tensiones es:\n");
Tensiones
printf ("%s", "el vector de corrientes resultante es:\n");
#
#	Multimplico a inzquierda por la matriz inversa de las impedancias.
#
%Corriente=inverse(Impedancias)*Tensiones;
Corriente = Impedancias \ Tensiones;
for k = 1:N
	#printf ("Corriente I%i: %f + %fj\n",k,real(Corriente(k,1)),imag(Corriente(k,1)));
	printf ("Corriente I%i:\n",k);
	Corriente(k,1)
endfor

#Corriente