function simmpplot

  cd Z:/Education/cvs/uio/postdoc/luminy/matlab
  set(gca,'FontSize',12);
  
  cvalues=[0.5 0.2 0.1 0.05];
  vals = [];
  for c=cvalues
    a=(1-sqrt(c))^2;
    b=(1+sqrt(c))^2;
    X = 0.02:0.02:3;
    newvals = sqrt( max(X-a,0) .* max(b-X,0) ) ./ ((2 *pi * c) .* X);
    vals = [vals ; newvals];
  end
  
  plot(X,vals(1,:),'k-',X,vals(2,:),'k--',X,vals(3,:),'k:',X,vals(4,:),'k-.');
  xlabel('')
  ylabel('Density')
  ax = axis;
  axis([0.01 3 0 1.6]);
  legend_handle = legend('c=0.5','c=0.2','c=0.1','c=0.05');
  set(legend_handle, 'Box', 'off');
  print -deps ../images/simmpplot.eps;