function [ zz,alphaa ] = calculetatadj(NPT,dt,PHIHUMc,Rad)
zz=linspace(1,2,NPT)*0.;
alphaa=linspace(1,2,NPT)*0.;
%==================
% donnees initiales
%==================
zz(1)=PHIHUMc(1);
zz(2)=zz(1)+dt*PHIHUMc(2);
alphaa(1)=PHIHUMc(3);
alphaa(2)=alphaa(1)+dt*PHIHUMc(4);
X1=[zz(1),alphaa(1)]';
X2=[zz(2),alphaa(2)]';
%=============
% resolution du systeme edo
%===========
for i=3:NPT
X3=2*X2-X1-dt^2*Rad*X2;
zz(i)=X3(1);
alphaa(i)=X3(2);
X1=X2;
X2=X3;
end
end

