icnd(ndim,par,icp,nint,u,uold,ufort,upold,fi,ijac,dint)

      function []= icnd(ndim,par,icp,nint,u,uold,ufort,upold,fi,ijac,dint)
%                    
%
% local
      parameter (ndm=2);
%
      call ffff(ndm,u   ,icp,par,1,f ,dfdu);
      call ffff(ndm,uold,icp,par,0,f0,dummy);
%
       fi(1)= ( f(1) - f0(1) ) * ( dfdu(1,1)*f(1) + dfdu(1,2)*f(2) );
     *      + ( f(2) - f0(2) ) * ( dfdu(2,1)*f(1) + dfdu(2,2)*f(2) );
%
      return