Maybe trying a simple calculation (X + n*epsilon ??) in FORTRAN and FiPy would be an instructive yet empirical approach?
It is simpler than dealing with complicated equations... Dan Lewis Assistant Professor Materials Research Center, Room 110 Rensselaer Polytechnic Institute http://www.rpi.edu/~lewisd2
