fpoly.f90

来自「FORTRANvisualfortran常用数值算法集及源码」· F90 代码 · 共 10 行

SUBROUTINE fpoly(x,p,np)
INTEGER np
REAL x,p(np)
INTEGER j
p(1)=1.
do j=2,np
  p(j)=p(j-1)*x
end do
END SUBROUTINE fpoly