#Arkaitz Rodas Bilbao
#© Derechos reservados excepto normas de uso y distribución de la UPV-EHU.

#---------------------------------------------------------------------------------------------------------------
#En este programa se realiza la traza parcial del caso 2n a dimension 2.
#Se emplea en el metodo runge-kutta para calcular la evolucion del qubit.
#---------------------------------------------------------------------------------------------------------------

disp ('');

function par=fu(B,n)

A=[0,0;0,0];

  for z=1: n
  
   A(1,1)=A(1,1)+B(z,z);
   A(1,2)=A(1,2)+B(z,n+z);
   A(2,1)=A(2,1)+B(n+z,z);
   A(2,2)=A(2,2)+B(n+z,n+z);
	
  endfor;

  par=A;
  
endfunction;

disp('');