function x=groveropt(n)
N=2^n
s=2*sqrt(N-1)/N;
c=1-2/N;
R=[c s;-s c];
d=1/sqrt(N);
u=sqrt((N-1)/N);
k=(pi/4)*sqrt(N)
for j=1:ceil(sqrt(N))
    du=[d u]';
    du=R*du;
    d=du(1);
    u=du(2);
    x(j)=d;
end
plot(x);
