REAL*4functionsamaxf(n,dx,incx)
c
cfindstheindexofelementhavingmax.absolutevalue.
cjackdongarra,linpack,3/11/78.
cmodified3/93toreturnifincx.le.0.
cmodified12/3/93,array(1)declarationschangedtoarray(*)
c
REAL*4dx(*)
integeri,incx,ix,n
c
samaxf=0.
if(n.lt.1.or.incx.le.0)return

if(incx.eq.1)go to20
c
ccodeforincrementnotequalto1
c
ix=1
samaxf=abs(dx(1))
ix=ix+incx
do10i=2,n
if(abs(dx(ix)).le.samaxf)go to5
samaxf=abs(dx(ix))
5ix=ix+incx
10continue
return
c
ccodeforincrementequalto1
c
20samaxf=abs(dx(1))
do30i=2,n
if(abs(dx(i)).le.samaxf)go to30
samaxf=abs(dx(i))
30continue
return
end