%
% FFIG2A
%
function y=fFig2a(r)
syms x R cr
dig=112; digits(dig);
R=vpa(r);
n=50;
for i=1:n
  x(i)=vpa(4*i/n);
end
cr=sVdMsc(dig,R,n,x);
c=subs(cr);
y=log10(c);
