関連情報

変形ベッセル関数 K0(x)

Fortranによるサンプルソースコード
使用ルーチン名:s18acf

Keyword: 変形ベッセル関数

概要

本サンプルは変形ベッセル関数 K0(x)を求めるFortranによるサンプルプログラムです。 本サンプルは引数xを読み込み、xの各値について変形ベッセル関数 K0(x)を求めて出力します。

※本サンプルはNAG Fortranライブラリに含まれるルーチン s18acf() のExampleコードです。本サンプル及びルーチンの詳細情報は s18acf のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで

入力データ

(本ルーチンの詳細はs18acf のマニュアルページを参照)

このデータをダウンロード
S18ACF Example Program Data
             0.4
             0.6
             1.4
             1.6
             2.5
             3.5
             6.0
             8.0
             10.0
             1000.0 

  • 1行目はタイトル行で読み飛ばされます。
  • 2〜11行目に変形ベッセル関数の引数xの値を指定しています。

出力結果

(本ルーチンの詳細はs18acf のマニュアルページを参照)

この出力例をダウンロード
 S18ACF Example Program Results

      X           Y

    4.000E-01   1.115E+00
    6.000E-01   7.775E-01
    1.400E+00   2.437E-01
    1.600E+00   1.880E-01
    2.500E+00   6.235E-02
    3.500E+00   1.960E-02
    6.000E+00   1.244E-03
    8.000E+00   1.465E-04
    1.000E+01   1.778E-05
    1.000E+03   0.000E+00

  • 3〜12行目に引数xの値と変形ベッセル関数の値が出力されています。

ソースコード

(本ルーチンの詳細はs18acf のマニュアルページを参照)

※本サンプルソースコードは科学技術・統計計算ライブラリである「NAG Fortranライブラリ」のルーチンを呼び出します。
サンプルのコンパイル及び実行方法


このソースコードをダウンロード
    PROGRAM s18acfe

!      S18ACF Example Program Text

!      Mark 23 Release. NAG Copyright 2011.

!      .. Use Statements ..
       USE nag_library, ONLY : nag_wp, s18acf
!      .. Implicit None Statement ..
       IMPLICIT NONE
!      .. Parameters ..
       INTEGER, PARAMETER              :: nin = 5, nout = 6
!      .. Local Scalars ..
       REAL (KIND=nag_wp)              :: x, y
       INTEGER                         :: ifail, ioerr
!      .. Executable Statements ..
       WRITE (nout,*) 'S18ACF Example Program Results'

!      Skip heading in data file
       READ (nin,*)

       WRITE (nout,*)
       WRITE (nout,*) '     X           Y'
       WRITE (nout,*)

DATA:  DO
          READ (nin,*,IOSTAT=ioerr) x

          IF (ioerr<0) THEN
             EXIT DATA
          END IF

          ifail = -1
          y = s18acf(x,ifail)

          IF (ifail<0) THEN
             EXIT DATA
          END IF

          WRITE (nout,99999) x, y
       END DO DATA

99999  FORMAT (1X,1P,2E12.3)
    END PROGRAM s18acfe


Results matter. Trust NAG.

Privacy Policy | Trademarks