from
slatec
by Ben Barrowes The slatec library converted into matlab functions.
[sbegresult,reset]=sbeg(reset);
function [sbegresult,reset]=sbeg(reset);
sbegresult=[];
persistent i ic mi ;
;
if isempty(i), i=0; end;
if isempty(ic), ic=0; end;
if isempty(mi), mi=0; end;
% intrinsic real ::;
if( reset )
mi = 891;
i = 7;
ic = 0;
reset = false;
end;
ic = fix(ic + 1);
i = fix(i.*mi);
i = fix(i - (1000.*(fix(i./1000))));
if( ic>=5 )
ic = 0;
i = fix(i.*mi);
i = fix(i - (1000.*(fix(i./1000))));
end;
sbegresult = real(i-500)./1001.0;
csnil=dbstack(1); csnil=csnil(1).name(1)~='@';
if csnil&&~isempty(inputname(1)), assignin('caller','FUntemp',reset); evalin('caller',[inputname(1),'=FUntemp;']); end
end %function sbeg