Scilab Website | Contribute with GitLab | Mailing list archives | ATOMS toolboxes
Scilab Online Help
2023.1.0 - 日本語


dawson

Dawson (スケーリングされた虚数の誤差) 関数を計算.

呼び出し手順

y = dawson(x)

引数

x
ベクトルまたは行列
y
(xと同じ大きさ)ベクトルまたは行列

説明

dawson はスケーリングされた虚数の誤差関数を計算します:

Dawson(z)=exp(-z^2).erfi(z).√π /2

i = %i;
x = [0, 0.01, 1, 100, i, 2*i, 4*i, 8*i,  1-i, 1+i, 2+i, 1e-4+2e-4*i].';
[x, dawson(x), dawson(-x)]
--> [x, dawson(x), dawson(-x)]
 ans  =
   0.     + 0.i       0.        + 0.i          0.        + 0.i
   0.01   + 0.i       0.0099993 + 0.i         -0.0099993 + 0.i
   1.     + 0.i       0.5380795 + 0.i         -0.5380795 + 0.i
   100.   + 0.i       0.0050003 + 0.i         -0.0050003 + 0.i
   0.     + i         0.        + 2.0300785i   0.        - 2.0300785i
   0.     + 2.i       0.        + 48.160012i   0.        - 48.160012i
   0.     + 4.i       0.        + 7875110.3i   0.        - 7875110.3i
   0.     + 8.i       0.        + 5.526D+27i   0.        - 5.526D+27i
   1.     - i         0.9903731 + 0.6388731i  -0.9903731 - 0.6388731i
   1.     + i         0.9903731 - 0.6388731i  -0.9903731 + 0.6388731i
   2.     + i         0.1635394 - 0.1531246i  -0.1635394 + 0.1531246i
   0.0001 + 0.0002i   0.0001    + 0.0002i     -0.0001    - 0.0002i

clf
subplot(1,2,1)
x = logspace(-20,20,100);
y = dawson(x);         // the imaginary part is zero
plot("ll", x, y)
xlabel x fontsize 3
title dawson(x) fontsize 3.5
xgrid(color("grey60"))

subplot(1,2,2)
x = logspace(-7,1,100);
y = imag(dawson(x*%i)); // the real part is zero
plot("ll", x, y)
xlabel x fontsize 3
title imag(dawson(x.i)) fontsize 3.5
xgrid(color("grey60"))

アルゴリズム

この関数は Faddeevaパッケージ ライブラリにもとづいています.

参照

  • erf — 誤差関数.
  • erfc — 相補誤差関数.
  • calerf — 誤差関数を計算する.

履歴

バージョン記述
5.5.0 関数dawsonが追加されました
Report an issue
<< calerf Special Functions delip >>

Copyright (c) 2022-2024 (Dassault Systèmes)
Copyright (c) 2017-2022 (ESI Group)
Copyright (c) 2011-2017 (Scilab Enterprises)
Copyright (c) 1989-2012 (INRIA)
Copyright (c) 1989-2007 (ENPC)
with contributors
Last updated:
Mon May 22 12:43:10 CEST 2023