function diff = women_fit ( param )
%
%  The parameter values are input in PARAM.
%
  mu = param(1);
  sigma = param(2);
%
%  Here is the data we are trying to match.
%
  cdata = [ 0.05, 0.10, 0.15, 0.25, 0.50, 0.75, 0.85, 0.90, 0.95];
  hdata = [ 59.9, 60.5, 61.4, 62.3, 64.1, 65.8, 67.0, 67.5, 68.4];
%
%  We evaluate our formula at the HDATA values.
%
  cformula = 0.5 + 0.5 * erf ( ( hdata - mu ) / sigma / sqrt ( 2 ) );
%
%  We compare our CDATA values to the output of the formula.
%
  diff = cdata - cformula; 

  return
end

  
