% EXAMPLE2_63 Routine used in Example 2.63.
%
f3='         zeros,  g=%3.4f\n';
fprintf('\n\n')
N=6; s=1; a0=1; same=0; nd=[N 2*N];
for g=[43646.1 43646.2]
  fprintf(f3,g)
  ab=r_jacobi(N); zw=gauss(N,ab);
  xw=zeros(2*N,2*(s+1));
  xw(1:N,1)=2*zw(:,1)+1; xw(1:N,2)=zw(:,1);
  xw(1:N,3)=2*zw(:,2); xw(1:N,4)=g*zw(:,2);
  xw(N+1:2*N,2)=zw(:,1)+2; xw(N+1:2*N,4)=zw(:,2);
  B=stieltjes_sob(N,s,nd,xw,a0,same);
  z=sobzeros(N,N,B)
  fprintf('\n')
end
