%
%  Initial parameter estimate.
%
  param(1:2) = [ 1, 1 ];

  param_new = lsqnonlin ( @fit_woman, param );

  mu = param_new(1);
  sigma = param_new(2);

  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];

  h = 58 : 0.25 : 70;
  c = 0.5 + 0.5 * erf ( ( h - mu ) / sigma / sqrt ( 2 ) );

  plot ( h, c, 'LineWidth', 2 );
  hold on
  bar ( hdata, cdata );

  title ( 'Compare CDF data and fitted formula' )
  xlabel ( 'Height' )
  ylabel ( 'CDF' )

  hold off

  p = exp ( - ( h - mu ).^2 / 2 / sigma^2 ) / sqrt ( 2 * pi * sigma^2 );
  plot ( h, p, 'LineWidth', 2 );
  title ( 'PDF using fitted MU and SIGMA' );
  xlabel ( 'Height' )
  ylabel ( 'PDF' )
