function [easting, northing] = LL2UTM(lat, lon)
utm = struct();
        
        utmstruct = defaultm('utm');
        utmstruct.zone = '10S';
        wgs84 = wgs84Ellipsoid;
        ellipsoid = [wgs84.SemimajorAxis wgs84.Eccentricity];
        utmstruct.geoid = ellipsoid;
        utm.struct = defaultm(utmstruct);
        
        [easting, northing] = mfwdtran(utm.struct, lat, lon);

end