Fixed function for geodist and bug in angdist matrix

old geodesic distance function did everything except for what it was
supposed to. This should come closer now. Also fixed indice mistake in
creation of the angular distance matrix. Results seem a lot better.