diff --git a/.vs/slnx.sqlite b/.vs/slnx.sqlite
new file mode 100644
index 0000000..c39f8d6
Binary files /dev/null and b/.vs/slnx.sqlite differ
diff --git a/curvspace_app.m b/curvspace_app.m
index 21544e0..c7d5d4a 100644
--- a/curvspace_app.m
+++ b/curvspace_app.m
@@ -12,7 +12,13 @@
 for k0 = 1:len-1
    dist_bet_pts(k0) = distance_app(app,p(k0,:),p(k0+1,:));
 end
-totaldist = nansum(dist_bet_pts);
+% totaldist = nansum(dist_bet_pts);
+% check version, use nansum for 2022 and earlier and sum w/omitnan option for 2023 and later
+if double(regexp(string(version('-release')), '^[0-9]+', 'match')) >= 2023
+	totaldist = sum(dist_bet_pts, "omitnan");
+else	
+	totaldist = nansum(dist_bet_pts);
+end
 
 %% interval %%
 intv = totaldist./(N-1);
diff --git a/grid_points_app.m b/grid_points_app.m
index ec93cdf..77d57c4 100644
--- a/grid_points_app.m
+++ b/grid_points_app.m
@@ -29,7 +29,13 @@
         for j=1:1:x10-1
             dist_steps(j)=deg2km(distance(border_dpa_bound(j,1),border_dpa_bound(j,2),border_dpa_bound(j+1,1),border_dpa_bound(j+1,2)));
         end
-        seg_dist=nansum(dist_steps);
+        % seg_dist=nansum(dist_steps);
+        % check version, use nansum for 2022 and earlier and sum w/omitnan option for 2023 and later
+        if double(regexp(string(version('-release')), '^[0-9]+', 'match')) >= 2023
+	        seg_dist=sum(dist_steps, "omitnan");
+        else	
+	        seg_dist=nansum(dist_steps);
+        end
         line_steps=ceil(seg_dist/(step_size))+1;
         dpa_edge_pt=curvspace_app(app,border_dpa_bound,line_steps);