% FIG1A
%
n=50; c=zeros(101,1); 
x=linspace(4/n,4,n)';
for k=1:101
  R=1+5*(k-1)/100;
  c(k)=VdMsc(R,n,x);
end
c
