/*  xevtch.c    CCMATH mathematics library source code.
 *
 *  Copyright (C)  2000   Daniel A. Atkinson    All rights reserved.
 *  This code may be redistributed under the terms of the GNU library
 *  public license (LGPL). ( See the lgpl.license file for details.)
 * ------------------------------------------------------------------------
 */
#include "xpre.h"
struct xpr xevtch(struct xpr z,struct xpr *a,int m)
{ struct xpr *p,f,t,tp,w;
  w=xpr2(z,1); t=zero; tp=zero;
  for(p=a+m; p>a ;){
    f=xadd(*p--,xadd(xmul(w,t),tp,1),0);
    tp=t; t=f;
   }
  return xadd(*p,xadd(xmul(z,t),tp,1),0);
}


syntax highlighted by Code2HTML, v. 0.9.1