Code covered by the BSD License  

Highlights from
slatec

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

[lun,kprint,ipass]=qc6a(lun,kprint,ipass);
function [lun,kprint,ipass]=qc6a(lun,kprint,ipass);
persistent ver versn ; 

if isempty(ver), ver=repmat(' ',1,16); end;
if isempty(versn), versn=' 4.1' ; end;
if( kprint>=3 )
writef(lun,['1', '\n ' , '\n ' ,' CODE TO TEST SLATEC ROUTINE AAAAAA' ' \n']);
end;
%format ('1'//' CODE TO TEST SLATEC ROUTINE AAAAAA');
[ver]=aaaaaa(ver);
if( strcmp(deblank(ver),deblank(versn)) )
ipass = 1;
if( kprint>=3 )
writef(lun,[ '\n ' ,' QC6A RESULTS' ' \n']);
writef(lun,[' *** Passed -- version numbermlv = ','%16s' ' \n'], ver);
%format (' *** Passed -- version numbermlv = ',a16);
end;
else;
ipass = 0;
if( kprint>=3 )
writef(lun,[ '\n ' ,' QC6A RESULTS' ' \n']);
end;
if( kprint>=2 )
writef(lun,[' *** Failed -- version numbermlv from AAAAAA = ','%16s',' but expected version numbermlv = ','%16s' ' \n'], ver , versn);
end;
%format (' *** Failed -- version numbermlv from AAAAAA = ',a16,' but expected version numbermlv = ',a16);
end;
if( kprint>=2 && ipass==1 )
writef(lun,[ '\n ' ,' **QC6A   PASSED  ALL TESTS *' ' \n']);
end;
%format [' **QC6A   PASSED  ALL TESTS *');
if( kprint>=1 && ipass==0 )
writef(lun,[ '\n ' ,' **QC6A   FAILED SOME TESTS *' ' \n']);
end;
%format [' **QC6A   FAILED SOME TESTS *');
return;
%format [' QC6A RESULTS');
end %subroutine qc6a

Contact us at files@mathworks.com