from
slatec
by Ben Barrowes The slatec library converted into matlab functions.
[dt5result,x]=dt5(x);
function [dt5result,x]=dt5(x);
dt5result=[];
persistent a b x1 y ;
;
if isempty(a), a=0; end;
if isempty(b), b=0; end;
if isempty(x1), x1=0; end;
if isempty(y), y=0; end;
a = 0.0d+00;
b = 0.1d+01;
x1 = x + 0.1d+01;
y =(b-a)./x1 + a;
dt5result =(b-a).*df5s(y)./x1./x1;
csnil=dbstack(1); csnil=csnil(1).name(1)~='@';
if csnil&&~isempty(inputname(1)), assignin('caller','FUntemp',x); evalin('caller',[inputname(1),'=FUntemp;']); end
end %function dt5