Code covered by the BSD License  

Highlights from
slatec

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

[kprint]=test23(kprint);
function [kprint]=test23(kprint);
persistent lin lun nerr nfail ; 

if isempty(lin), lin=0; end;
if isempty(lun), lun=0; end;
if isempty(nerr), nerr=0; end;
if isempty(nfail), nfail=0; end;
[lun ]=i1mach(2);
[lin ]=i1mach(1);
nfail = 0;
xermax(1000);
[lun]=xsetun(lun);
if( kprint<=1 )
xsetf(0);
else;
xsetf(1);
end;
[lun,kprint,nerr]=cgeck(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=cgbqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=cpoqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=cppqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=cpbqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=csiqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=cspqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=chiqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=chpqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=ctrqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=cgtqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=cptqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=cchqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=cqrqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
[lun,kprint,nerr]=csvqc(lun,kprint,nerr);
nfail = fix(nfail + nerr);
if( nfail==0 )
writef(lun,['%s'], 'TEST23 passed all tests.');
else;
writef(lun,[ '\n ' ,' *** WARNING -- ','%5i',' TEST(S) FAILED IN PROGRAM TEST23 ***' ' \n'], nfail);
%format [' *** WARNING -- ',i5,' TEST(S) FAILED IN PROGRAM TEST23 ***');
end;
return;
end %subroutine test23

Contact us