I found that it's a bug in
.../toolbox/map/map/private/geodesicinv.m
Two places of
flat = ecc2flat(ellipsoid);
should be changed to
flat = ecc2flat(ellipsoid(2));
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!