README

The  DOS  DDAVE.EXE  program  will  compute any size parameter on any AT 
compatible  machine  with  computation  time  the  expense  paid for the 
benefit  of  unlimited  size parameters.  The ratios of spherical Bessel 
functions  are  recomputed as they are needed in zratio.for in blocks of 
NDIM  that  are  limited by 640K for DOS.  The An values are computed by 
backward  recursion  in  ratios  of  spherical Bessel functions from the 
single  Bessel  ratio  calculated  at  the  upper  limit  of NDIM.  When 
Dave.for  runs  out  of  An  values,  it  recomputes  the next ratio and 
computes another block by backward recursion.  

Instructions for the use of the program ddave 

The following is a sample input and output

Input file ddave < test.in > test.out

*************test.in**************************
0,10,1.5,0
5,10,1.5,0.01
1
10
0,100000,1.5,0.0001

***************test.out*********************
Ntheta=number of angles to 2001
x=size parameter
rfr=real index, rfi=(positive)imaginary index
ztheta=starting angle or its complement 180-ztheta
(dtheta + 1)*NTHETA + ztheta must be LE 90 DEGREES
INPUT NTHETA,X,RFR,RFI   
RFR=  .150000E+01  RFI=  .000000E+00  X=  .100000E+02
Working
Seconds=    0.000000E+00
QEXT=  .288200E+01  QSCAT=  .288200E+01  QABS=  .000000E+00
RFR=  .150000E+01  RFI=  .000000E+00  X=  .100000E+02
QEXT=  .288200E+01  QSCAT=  .288200E+01  QABS=  .000000E+00
Seconds=    0.000000E+00
Angles    .000000E+00    .000000E+00    .000000E+00
Complementary Angles    .180000E+03    .180000E+03    .000000E+00
Ntheta=number of angles to 2001
x=size parameter
rfr=real index, rfi=(positive)imaginary index
ztheta=starting angle or its complement 180-ztheta
(dtheta + 1)*NTHETA + ztheta must be LE 90 DEGREES
INPUT NTHETA,X,RFR,RFI   
Input angle increment in degrees OR final angle <_ 90 
Input starting angle 
Angles    .100000E+02    .140000E+02    .100000E+01
Complementary Angles    .170000E+03    .166000E+03    .100000E+01
RFR=  .150000E+01  RFI=  .100000E-01  X=  .100000E+02
DTHETA=       1.000000000000000
Working
Seconds=    5.000000E-02
QEXT=  .277070E+01  QSCAT=  .234413E+01  QABS=  .426563E+00
 10.00000    .196327E+04    .225380E+04    .210255E+04    .638859E+02
 11.00000    .161332E+04    .190899E+04    .175375E+04    .644477E+02
 12.00000    .129659E+04    .158639E+04    .143282E+04    .627032E+02
 13.00000    .101796E+04    .129175E+04    .114519E+04    .589812E+02
 14.00000    .779655E+03    .102894E+04    .894052E+03    .537549E+02
166.00000    .465448E+02    .170237E+02    .185793E+02   -.211465E+02
167.00000    .416948E+02    .132354E+02    .186219E+02   -.143203E+02
168.00000    .359800E+02    .989617E+01    .174145E+02   -.726638E+01
169.00000    .301525E+02    .741010E+01    .149349E+02   -.619230E+00
170.00000    .249411E+02    .609935E+01    .112525E+02    .505042E+01
RFR=  .150000E+01  RFI=  .100000E-01  X=  .100000E+02
QEXT=  .277070E+01  QSCAT=  .234413E+01  QABS=  .426563E+00
Seconds=    5.000000E-02
Angles    .100000E+02    .140000E+02    .100000E+01
Complementary Angles    .170000E+03    .166000E+03    .100000E+01
Ntheta=number of angles to 2001
x=size parameter
rfr=real index, rfi=(positive)imaginary index
ztheta=starting angle or its complement 180-ztheta
(dtheta + 1)*NTHETA + ztheta must be LE 90 DEGREES
INPUT NTHETA,X,RFR,RFI   
RFR=  .150000E+01  RFI=  .100000E-03  X=  .100000E+06
Working
Seconds=       22.680000
QEXT=  .200092E+01  QSCAT=  .109262E+01  QABS=  .908307E+00
RFR=  .150000E+01  RFI=  .100000E-03  X=  .100000E+06
QEXT=  .200092E+01  QSCAT=  .109262E+01  QABS=  .908307E+00
Seconds=       22.680000
Angles    .100000E+02    .140000E+02    .100000E+01
Complementary Angles    .170000E+03    .166000E+03    .100000E+01
Ntheta=number of angles to 2001
x=size parameter
rfr=real index, rfi=(positive)imaginary index
ztheta=starting angle or its complement 180-ztheta
(dtheta + 1)*NTHETA + ztheta must be LE 90 DEGREES
INPUT NTHETA,X,RFR,RFI   

**************************************************
The zratio routine has JFIX in it to produce
integer*4 output.  You may need to change this

See the file ddave.let for more detailed information
