icant.f

   1       integer function icant(i,j)
   2       integer i,j
   3       if (j.ge.i) then
   4         icant=i*(i-1)/2+j
   5       else
   6         icant=j*(j-1)/2+i
   7       endif
   8       return
   9       end