lat1 = 36.0;
lon1 = -123.0;

lat2 = 36.0166667;
lon2 = -123.0;

earthR = 6378137;

[distAng heading] = distance(lat1, lon1, lat2, lon2)
range = sin(distAng*pi/180.0) * earthR