
utmstruct = defaultm('utm');
utmstruct.zone = '42G';
wgs84 = wgs84Ellipsoid;
ellipsoid = [wgs84.SemimajorAxis wgs84.Eccentricity];
utmstruct.geoid = ellipsoid;
utmstruct  = defaultm(utmstruct);

latDecDeg = [ -47 -46.75];
lonDecDeg = [66.5 66.5122];
[east, north] = mfwdtran(utmstruct, latDecDeg, lonDecDeg);

distance1 = sqrt( (east(2)-east(1))^2 + (north(2)-north(1))^2)
