latlon2xy and xy2lonlat not exactly inverse of each other

[kmaterna]

One function uses haversine, one uses an assumption of the size of a degree. The swap-test returns nearly but not exactly the same number. Should fix in future version.