Keyword: 変形ベッセル関数
概要
本サンプルは変形ベッセル関数 K0(x)を求めるFortranによるサンプルプログラムです。 本サンプルは引数xを読み込み、xの各値について変形ベッセル関数 K0(x)を求めて出力します。
※本サンプルはnAG Fortranライブラリに含まれるルーチン s18acf() のExampleコードです。本サンプル及びルーチンの詳細情報は s18acf のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで
入力データ
（本ルーチンの詳細はs18acf のマニュアルページを参照）1 2 3 4 5 6 7 8 9 10 11
| このデータをダウンロード |
S18ACF Example Program Data
0.4
0.6
1.4
1.6
2.5
3.5
6.0
8.0
10.0
1000.0
- １行目はタイトル行で読み飛ばされます。
- ２～１１行目に変形ベッセル関数の引数xの値を指定しています。
出力結果
（本ルーチンの詳細はs18acf のマニュアルページを参照）1 2 3 4 5 6 7 8 9 10 11 12 13 14
| この出力例をダウンロード |
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
- ３～１２行目に引数xの値と変形ベッセル関数の値が出力されています。
ソースコード
（本ルーチンの詳細はs18acf のマニュアルページを参照）
※本サンプルソースコードは科学技術・統計計算ライブラリである「nAG Fortranライブラリ」のルーチンを呼び出します。
サンプルのコンパイル及び実行方法
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44
| このソースコードをダウンロード |
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