Code covered by the BSD License  

Highlights from
slatec

from slatec by Ben Barrowes
The slatec library converted into matlab functions.

[f1wsresult,x]=f1ws(x);
function [f1wsresult,x]=f1ws(x);
f1wsresult=[];

;
f1wsresult = abs(x-0.33e+00).^(-0.999e+00);
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 f1ws

Contact us at files@mathworks.com