function functn(x,i,a) implicit double precision (a-h,o-z) dimension x(100),a(10) xi=x(i) functn=((((a(1)*xi+a(2))*xi+a(3))*xi+a(4))*xi+a(5))*xi+a(6) return end