%
% SVDMSC3  Example 3.3
%
syms a R x
dig=48; digits(dig);
a=vpa(1)
R=vpa(1);
c=zeros(41,1);
for n=10:50
  n
  for i=1:n
    x(i)=a*cos(vpa((2*i-1)*pi/(2*n)));
  end
  cVsc=sVdMsc(dig,R,n,x);
  c(n-9)=subs(cVsc);
end
c
