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

#-------------------------------------------------------------------------------------
#En este modulo se interpola el valor del hamiltoniano para un tiempo dado.
#-------------------------------------------------------------------------------------

disp ('');

function interp=fi(ht,tl,t,n,nt)

for i=1: 2*n
  for j=1: 2*n
  
   H(i,j)=interp1 ( tl, ht{i,j}, t);
  endfor;
 endfor;

 interp=H;
 
endfunction;
 
disp ('');