function [ sol ] = resoladj( npt,dt,TCM,TKM,COND)
%=====
sol(1,1)=COND(1);sol(2,1)=COND(2);
sol(1,2)=COND(1)+dt*COND(3);sol(2,2)=COND(2)+dt*COND(4);
%=====
for kt=3:npt
        
        for ii=1:2
            sum=0;
            for jj=1:2
            sum=sum-dt*TCM(ii,jj)*(sol(jj,kt-1)-sol(jj,kt-2))+dt^2*TKM(ii,jj)*sol(jj,kt-1);
            end
            sol(ii,kt)=2*sol(ii,kt-1)-sol(ii,kt-2)-sum;
        end  
        
    
        
end

