/* bc program secant */
/* This returns the secant number s[n]. In fact it returns s[0],...,s[n].
   See Algorithm 4.5,  https://members.loria.fr/PZimmermann/mca/pub226.html
   by Richard Brent and Paul Zimmermann.
 */
define secant(n){
auto k,j
   s[0]=1
   for(k=1;k<=n;k++){
      s[k]=k*s[k-1]
   }
   for(k=1;k<=n;k++){
      kplus1=k+1
      for(j=kplus1;j<=n;j++){
         s[j]=(j-k)*s[j-1]+(j-k+1)*s[j]
      }
   }
   for(k=0;k<n;k++){
      print "s[",k,"]=",s[k]," "
   }
   print "s[",n,"]="
   return(s[n]) 
}