5 DOUBLE PRECISION FUNCTION ddwnrm(NEQ,V,RWT,RPAR,IPAR)
20 IMPLICIT DOUBLE PRECISION(a-h,o-z)
26 IF(
abs(v(i)*rwt(i)) .GT. vmax) vmax =
abs(v(i)*rwt(i))
28 IF(vmax .LE. 0.0d0)
GO TO 30
31 20 sum = sum + ((v(i)*rwt(i))/vmax)**2
32 ddwnrm = vmax*sqrt(sum/neq)
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)