#	fortran equivalent of veclib smax function
real*4 function smax (n,x,incx)
integer n,incx,incxa,i,ix
real*4 x(n)

smax=-1.e+30
incxa=abs(incx)
ix=1
do i=1,n	{
	smax=max(smax,x(ix))
	ix=ix+incxa
		}
return 
end