Cases prefixed with W are classified as warnings and do not generate an error of type NAG:error_n. See nag_issue_warnings.
Open in the MATLAB editor: s18ct_example
function s18ct_example fprintf('s18ct example results\n\n'); x = [0; 0.5; 1; 3; 6; 10; 1000; -1]; [f, ifail] = s18ct(x); fprintf(' x e^-|x| I_1(x)\n'); for i=1:numel(x) fprintf('%12.3e%12.3e\n', x(i), f(i)); end
s18ct example results x e^-|x| I_1(x) 0.000e+00 0.000e+00 5.000e-01 1.564e-01 1.000e+00 2.079e-01 3.000e+00 1.968e-01 6.000e+00 1.521e-01 1.000e+01 1.213e-01 1.000e+03 1.261e-02 -1.000e+00 -2.079e-01