[pimachresult,dum]=pimach(dum);

function [pimachresult,dum]=pimach(dum);
pimachresult=[];
persistent pimach ; 

if isempty(pimachresult), pimachresult=0; end;
%***BEGIN PROLOGUE  PIMACH
%***SUBSIDIARY
%***PURPOSE  Subsidiary to HSTCSP, HSTSSP and HWSCSP
%***LIBRARY   SLATEC
%***TYPE      SINGLE PRECISION (PIMACH-S)
%***AUTHOR  (UNKNOWN)
%***DESCRIPTION
%
%     This subprogram supplies the value of the constant PI correct to
%     machine precision where
%
%     PI=3.1415926535897932384626433832795028841971693993751058209749446
%
%***SEE ALSO  HSTCSP, HSTSSP, HWSCSP
%***ROUTINES CALLED  (NONE)
%***REVISION HISTORY  (YYMMDD)
%   801001  DATE WRITTEN
%   891214  Prologue converted to Version 4.0 format.  (BAB)
%   900402  Added TYPE section.  (WRB)
%***end PROLOGUE  PIMACH
%
%***FIRST EXECUTABLE STATEMENT  PIMACH
pimachresult = 3.14159265358979;
csnil=dbstack(1); csnil=csnil(1).name(1)~='@';
if csnil&&~isempty(inputname(1)), assignin('caller','FUntemp',dum); evalin('caller',[inputname(1),'=FUntemp;']); end
end
%DECK PINITM