Keyword: ケルビン関数
概要
本サンプルはケルビン関数 kei xを求めるFortranによるサンプルプログラムです。 本サンプルは引数xを読み込み、xの各値についてケルビン関数 kei xを求めて出力します。
※本サンプルはnAG Fortranライブラリに含まれるルーチン s19adf() のExampleコードです。本サンプル及びルーチンの詳細情報は s19adf のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで
入力データ
(本ルーチンの詳細はs19adf のマニュアルページを参照)- 1行目はタイトル行で読み飛ばされます。
- 2〜8行目にケルビン関数の引数xの値を指定しています。
出力結果
(本ルーチンの詳細はs19adf のマニュアルページを参照)| この出力例をダウンロード |
S19ADF Example Program Results
X Y
0.000E+00 -7.854E-01
1.000E-01 -7.769E-01
1.000E+00 -4.950E-01
2.500E+00 -1.107E-01
5.000E+00 1.119E-02
1.000E+01 -3.075E-04
1.500E+01 7.963E-06
- 5〜11行目に引数xの値とケルビン関数の値が出力されています。
ソースコード
(本ルーチンの詳細はs19adf のマニュアルページを参照)
※本サンプルソースコードは科学技術・統計計算ライブラリである「nAG Fortranライブラリ」のルーチンを呼び出します。
サンプルのコンパイル及び実行方法
| このソースコードをダウンロード |
PROGRAM s19adfe
! S19ADF Example Program Text
! Mark 23 Release. nAG Copyright 2011.
! .. Use Statements ..
USE nag_library, ONLY : nag_wp, s19adf
! .. 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,*) 'S19ADF 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 = s19adf(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 s19adfe
