% SEARSKARMAN_ISTARAL  Drag coefficient of the Sears-Haack
%   and the von Karman body
%
global alpha beta abj abj1 abj2 abj3 abj4 abl1 abl2 abl3 abl4 
f0='%8.0f %17.12f %17.12f\n';
N=250; alpha=-1/2; beta=-1/2;
abj=r_jacobi01(N,alpha,beta);
abj1=r_jacobi01(N,alpha,0); abj2=r_jacobi01(N,beta,0);
abj3=r_jacobi01(N,alpha,2*beta+1); 
abj4=r_jacobi01(N,beta,2*alpha+1);
abl1=r_jacobilog(N,alpha,0); abl2=r_jacobilog(N,beta,0);
abl3=r_jacobilog(N,alpha,2*beta+1);
abl4=r_jacobilog(N,beta,2*alpha+1);
for n=50:50:250
  [Istara,Istarl]=Istaral(n);
  fprintf(f0,n,Istara,Istarl) 
end
