[dhvnrmresult,v,ncomp]=dhvnrm(v,ncomp);

function [dhvnrmresult,v,ncomp]=dhvnrm(v,ncomp);
dhvnrmresult=[];
persistent k ; 

;
%***BEGIN PROLOGUE  DHVNRM
%***SUBSIDIARY
%***PURPOSE  Subsidiary to DDEABM, DDEBDF and DDERKF
%***LIBRARY   SLATEC
%***TYPE      doubleprecision (HVNRM-S, DHVNRM-D)
%***AUTHOR  Watts, H. A., (SNLA)
%***DESCRIPTION
%
%     Compute the maximum norm of the vector V(*) of length NCOMP and
%     return the result as DHVNRM
%
%***SEE ALSO  DDEABM, DDEBDF, DDERKF
%***ROUTINES CALLED  (NONE)
%***REVISION HISTORY  (YYMMDD)
%   820301  DATE WRITTEN
%   890531  Changed all specific intrinsics to generic.  (WRB)
%   890831  Modified array declarations.  (WRB)
%   891024  Changed references from DVNORM to DHVNRM.  (WRB)
%   891024  Changed routine name from DVNORM to DHVNRM.  (WRB)
%   891214  Prologue converted to Version 4.0 format.  (BAB)
%   900328  Added TYPE section.  (WRB)
%   910722  Updated AUTHOR section.  (ALS)
%***end PROLOGUE  DHVNRM
%
if isempty(k), k=0; end;
v_shape=size(v);v=reshape(v,1,[]);
%***FIRST EXECUTABLE STATEMENT  DHVNRM
dhvnrmresult = 0.0d0;
for k = 1 : ncomp;
dhvnrmresult = max(dhvnrmresult,abs(v(k)));
end; k = fix(ncomp+1);
v_shape=zeros(v_shape);v_shape(:)=v(1:numel(v_shape));v=v_shape;
csnil=dbstack(1); csnil=csnil(1).name(1)~='@';
if csnil&&~isempty(inputname(1)), assignin('caller','FUntemp',v); evalin('caller',[inputname(1),'=FUntemp;']); end
if csnil&&~isempty(inputname(2)), assignin('caller','FUntemp',ncomp); evalin('caller',[inputname(2),'=FUntemp;']); end
end
%DECK DINTP