a = 6378137;
f = 1/298.257223563;
e = sqrt(2*f - f^2);
ellipsoid = [a e];
utmstruct= defaultm('utm');
utmstruct.zone = '10N';
utmstruct.geoid = ellipsoid;
utmstruct = defaultm(utm(utmstruct));


lat1 = 36 + 42.740/ 60.0;
lon1 = -(122 + 11.210/60.0);

[e, n] = mfwdtran(utmstruct, lat1, lon1)

[lat2, lon2] = minvtran(utmstruct, e, n)


lat3 = 36 + 41.439/ 60.0;
lon3 = -(122 + 09.962/60.0);

[e2, n2] = mfwdtran(utmstruct, lat3, lon3)

[lat4, lon4] = minvtran(utmstruct, e2, n2)

dist = sqrt((e-e2)^2 + (n-n2)^2)
