lat1Deg = 36.0;
lat1Min = 40.99793;
lat1DD = lat1Deg+lat1Min/60.0;

lon1Deg = 121;
lon1Min = 56.7148;
lon1DD = lon1Deg+lon1Min/60.0;
lon1DD = -1.0 * lon1DD;

[easting1, northing1] = LatLon2UTM(lat1DD, lon1DD)

lat2Deg = 36.0;
lat2Min = 40.97765;
lat2DD = lat2Deg+lat2Min/60.0;

lon2Deg = 121;
lon2Min = 56.61755;
lon2DD = lon2Deg+lon2Min/60.0;
lon2DD = -1.0 * lon2DD;

[easting2, northing2] = LatLon2UTM(lat2DD, lon2DD)

range = sqrt( (easting2 - easting1)^2 + (northing2 - northing1)^2 )
bearing  = atan2((northing2 - northing1), (easting2 - easting1))

predictedEasting = easting2 + (range * cos(bearing))
predictedNorthing = northing2 + (range * sin(bearing))

[lat, lon] = UTM2LatLon(predictedEasting, predictedNorthing)

latMin = (lat - fix(lat)) * 60.0
lonMin = (lon - fix(lon)) * 60.0

