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

#--------------------------------------------------------------------------------------------------------
#Calculamos en esta subrutina el disipador de la evolucion de la matriz densidad del sistema.
#--------------------------------------------------------------------------------------------------------

disp('');

function dis=fD(a,d,g,n1,n2);

D=0;

for k=1:n1
 for i=1:n2
  for j=1:n2
   D=D+g(i,j,k)*(a{k,j}*d*transpose(conj(a{k,i}))-1/2*(transpose(conj(a{k,i}))*a{k,j}*d+d*transpose(conj(a{k,i}))*a{k,j}));
  endfor;
 endfor;
endfor;

dis=D;

endfunction;

disp('');