CDF and Inverse CDF of Wrapped Cauchy Distribution

   function main
    p    = linspace(0,1,11);
    prob = linspace(0,1,100);
    x    = zeros(numel(p),numel(prob));
    for i = 1:numel(p)
      for j = 1:numel(prob)
        x(i,j) = INVCDF(p(i),prob(j));
      end
    end  
    plot(prob,x(1,:),prob,x(2,:),prob,x(3,:),prob,x(4,:),prob,x(5,:),...
         prob,x(6,:),prob,x(7,:),prob,x(8,:),prob,x(9,:),prob,x(10,:),...
         prob,x(11,:))
   end

   function x = INVCDF(p,prob)
     fun = @(x)CDF(p,x)-prob;
     x = fzero(fun,pi);  
   end

   function prob = CDF(p,x) 
     if x