duclid.f

来自「最优化算法」· F 代码 · 共 16 行

c
c-----------------------------------------------------------------------
c
      subroutine duclid (n,x,y,ps,izs,rzs,dzs)
      integer n,izs(1)
      real rzs(1)
      double precision x(n),y(n),ps,dzs(1)
c
      integer i
c
      ps=0.d0
      do 10 i=1,n
   10 ps=ps+x(i)*y(i)
      return
      end