Plotting Climate Division by Latitude and Longitude on Map

Hello. I have a map of the state of Nebraska I plotted by using the following:

figure(1) 
usamap('Nebraska')
nebraskahi = shaperead('usastatehi', 'UseGeoCoords', true,...
            'Selector',{@(name) strcmpi(name,'Nebraska'), 'Name'});
geoshow(nebraskahi, 'FaceColor', [0.3 1.0, 0.675])
mlabel off; plabel off; gridm off
hold on

and on that map I want to outline one of the climate divisions in Nebraska. I have the longitude and latitude of that climate division boundary here

lat =[43           43           43      42.8359      42.8305      42.6631      42.1351      42.1123      42.0358       42.032      41.7492      41.7455      41.3734      41.3734      41.3895           41           41           41           41      42.2557      42.2557           43];
lon=[-104.03     -102.0276     -102.0276      -102.035     -101.9827     -101.9971     -102.0106     -101.9682     -101.9652     -102.0229     -102.0172     -101.9539     -101.9553     -101.9553     -102.0146     -102.0302     -102.0302       -104.02       -104.02     -104.0263     -104.0263       -104.03];

but can't figure out the right way to plot this onto a map. Anything helps, thank you in advance.

[Michaela Wood]

Using

geoplot(lat,lon)

does work to plot the these division lines however I get the error

Error using geoplot (line 98)
Adding GeographicAxes to axes is not supported. Turn
hold off.

Error in f01_divisions (line 35)
geoplot(lat1,lon1)

and Im not sure how to resolve this.

Answer 1

[Cris LaPierre]

Specifying the axes in geoshow appears to work for me.

ax = usamap('Nebraska');
nebraskahi = shaperead('usastatehi', 'UseGeoCoords', true,...
            'Selector',{@(name) strcmpi(name,'Nebraska'), 'Name'});
geoshow(nebraskahi, 'FaceColor', [0.3 1.0, 0.675])
mlabel off; plabel off; gridm off

lat =[43           43           43      42.8359      42.8305      42.6631      42.1351      42.1123      42.0358       42.032      41.7492      41.7455      41.3734      41.3734      41.3895           41           41           41           41      42.2557      42.2557           43];
lon=[-104.03     -102.0276     -102.0276      -102.035     -101.9827     -101.9971     -102.0106     -101.9682     -101.9652     -102.0229     -102.0172     -101.9539     -101.9553     -101.9553     -102.0146     -102.0302     -102.0302       -104.02       -104.02     -104.0263     -104.0263       -104.03];

geoshow(ax,lat,lon)

[Michaela Wood]

Thank you!