function Input = Input_one(KAPPA,PHI,GAMMA,LAMBDA,ALPHAk,ALPHA,ETA,DELT,DELTT,e,SIGMAv,BETA,Dir)
   N = length(ALPHA);
   In = zeros(N,N);
   Int = Input_bis(KAPPA,PHI,GAMMA,LAMBDA,ALPHAk,ALPHA,ETA,DELT,DELTT,e,SIGMAv,BETA,Dir);
     for i=1:N
       for j=1:N
           for r=1:N
             In(j,i)=In(j,i)+GAMMA((j-1)*N+r)*(1-SIGMAv(r))*LAMBDA(r)*(-Dir(j,i)+1/(1-ALPHA(r)-BETA(r))*Int(r,i));
           end
       end
     end
  
Input = In;
 
end 