c nxval=number of original data (x) values.
c ndaf=largest desired derivative of the fitted function.
c mval=order of the daf fit. This is M.
c ndelta=generate the new function at nxval*ndelta values.
c nchanl=number of channels if f(x) is a vector or matrix.
c nder=number of known derivatives in the original function.
c      Note nder serves no useful purpose at the present.
 $daf
   nxval=301,
   ndaf=1
   mval=32
   ndelta=1
   nchanl=1,
   nder=1
 $end
 $optkase
   kase=10
 $end