% Testprogramm zu SubDiv

n = 5

for j=0:n-1
  a(j+1) = (j==0)/4 + (3+2*cos(2*pi*j/n))/4/n;
  %a{j+1}(1,1) = (1+cos(2*pi*j/n))/n;
  %a{j+1}(1,1) = (j==0)*.98 + (.02 - 0.02*cos(2*pi*j/n))/n;
end

A = SubMat(a);

[lambda, V] = CharMap(A);

x = V(:,1);
y = 1.5*V(:,2);
z = x.^2 - y.^2;

figure(1)
clf reset
P = [x,y,z];
RingMesh(P)

for i=1:5
  P = SubDiv(P,A);
  RingMesh(P)
end
cameratoolbar


figure(2)
clf reset
P = [x,y,z];
RingSurf(P)

for i=1:5
  P = SubDiv(P,A);
  RingSurf(P)
end
cameratoolbar