% fator de deslocamento alfa
% pa(1)=TT1 temperatura de trabalho (K)
% pa(2)=TT0 temperatura referencia (K)
% pa(3)=teta1 coeficiente temperatura
% pa(4)=teta2 coeficiente temperatura
function alf=alfa(pa)
  % inicio
  dt=(pa(1)-pa(2));
  alf=10^(-pa(3)*dt/(pa(4)+dt)); 
  % fim
end