  Our aim in this paper is to construct a numerical algorithm using Taylor
expansion of exponential map to find geodesic joining two points on a
2-dimensional surface for which a Riemannian metric is defined
