%
% SFIG1B
%
syms x R cr 
dig=16; digits(dig);
c=zeros(101,1); 
n=50;
for i=1:n
  x(i)=vpa(4*i/n);
end
for k=1:101
  R=vpa(1+5*(k-1)/100);
  subs(R)
  cr=sVdMsc(dig,R,n,x);
  c(k)=subs(cr);
end
c
