grantham_equation function

Grantham distance