1: $result$ – double scalar: The result of the function.
2: $ifail$ – int64int32nag_int scalar: $ifail = 0$ unless the function detects an error (see Error Indicators and Warnings).

Error Indicators and Warnings

There are no failure exits from this function. The argument ifail is included for consistency with other functions in this chapter.

Accuracy

Because of its close relationship with

erfc

the accuracy of this function is very similar to that in nag_specfun_erfc_real (s15ad). If

ε

and

δ

are the relative errors in result and argument, respectively, then in principle they are related by

|ε| ≃ |\frac{x e^{- x^{2} / 2}}{\sqrt{2 π} Q (x)} δ| .

For

x

negative or small positive this factor is always less than one and accuracy is mainly limited by machine precision. For large positive

x

we find

ε \sim x^{2} δ

and hence to a certain extent relative accuracy is unavoidably lost. However the absolute error in the result,

E

, is given by

|E| ≃ |\frac{x e^{- x^{2} / 2}}{\sqrt{2 π}} δ|

and since this factor is always less than one absolute accuracy can be guaranteed for all

x

Further Comments

None.

Example

This example reads values of the argument

x

from a file, evaluates the function at each value of

x

and prints the results.

Open in the MATLAB editor: s15ac_example

function s15ac_example


fprintf('s15ac example results\n\n');

x = [-20   -1     0    1    2    20];
n = size(x,2);
result = x;

for j=1:n
  [result(j), ifail] = s15ac(x(j));
end

disp('      x          Q(x)');
fprintf('%12.3e%12.3e\n',[x; result]);

s15ac example results

      x          Q(x)
  -2.000e+01   1.000e+00
  -1.000e+00   8.413e-01
   0.000e+00   5.000e-01
   1.000e+00   1.587e-01
   2.000e+00   2.275e-02
   2.000e+01   2.754e-89

PDF version (NAG web site, 64-bit version, 64-bit version)

Chapter Contents

Chapter Introduction

NAG Toolbox