plot.gravity function

Plot gravity model