function f = aasmaResidual(Ma,ar,k)
% aasma_residual   Residual of relationship between A/Astar and Ma
%                  used for root-finding routines

if nargin<3,  k = 1.4;  end    %  default fluid is air

f = ar - isenAAs(Ma,k);