%EXAMPLE4
%
f0='%10.0f %6.0f %18.12f %14.4e\n';
f1='%16.0f %18.12f %14.4e\n';
dig=52; digits(dig);
%dig=48; digits(dig);
N=30;
ab0=sr_jacobi(dig,N,-1/2);
for k=1:2:47
  z(k)=vpa(1+2/(k+1)); z(k+1)=-z(k);
%  z(k)=vpa(i/k); z(k+1)=-z(k);
end
for m=2:2:48
  ab=mod_md(N,m,z,ab0,@rho0);
  ab=subs(real(ab));
  fprintf('m=%2.0f',m)
  [ab(1:3+m/2,2);ab(N,2)]
end
