function ps23(n)
% function ps23(n)
%
% random trials for problem 2.3P: 
%    n blocks, prod(h) in [0,1]

if nargin<1, n=2; end
N=10;
gg=zeros(N,1);
a=cos(pi/n)^n;
pp=linspace(a,1,N)';
for k=1:N,
    h=rand(n,1);
    h=h*(pp(k)/prod(h))^(1/n);
    gg(k)=ps23b(h);
end

ga=max(0,1-a./pp);
close(gcf);plot(pp,gg,'.',pp,ga); grid