Keyword: ３次, スプライン
概要
本サンプルは３次スプラインの評価を行うFortranによるサンプルプログラムです。 本サンプルは以下に示されるデータについて３次スプラインの値を求めて出力します。
※本サンプルはnAG Fortranライブラリに含まれるルーチン e02bbf() のExampleコードです。本サンプル及びルーチンの詳細情報は e02bbf のマニュアルページをご参照ください。
ご相談やお問い合わせはこちらまで
入力データ
（本ルーチンの詳細はe02bbf のマニュアルページを参照）| このデータをダウンロード |
E02BBF Example Program Data
9
4
1.00
1.00
1.00
1.00
3.00
6.00
8.00
9.00
9.00
9.00
9.00
1.00
2.00
4.00
7.00
6.00
4.00
3.00
- １行目はタイトル行で読み飛ばされます。
- ２行目にスプラインを計算する点の数(m)を指定しています。
- ３行目にスプラインが定義される区間の数(ncap)を指定しています。
- ４〜１４行目に拡張されたノット(lamda)を指定しています。
- １５〜２１行目にBスプライン係数(c)を指定しています。
出力結果
（本ルーチンの詳細はe02bbf のマニュアルページを参照）| この出力例をダウンロード |
E02BBF Example Program Results J LAMDA(J) B-spline coefficient (J-2) 1 1.0000 2 1.0000 3 1.0000 1.0000 4 1.0000 2.0000 5 3.0000 4.0000 6 6.0000 7.0000 7 8.0000 6.0000 8 9.0000 4.0000 9 9.0000 3.0000 10 9.0000 11 9.0000 R Argument Value of cubic spline 1 1.0000 1.0000 2 2.0000 2.3779 3 3.0000 3.6229 4 4.0000 4.8327 5 5.0000 5.8273 6 6.0000 6.3571 7 7.0000 6.1905 8 8.0000 5.1667 9 9.0000 3.0000
- ５〜１５行目に拡張ノットとBスプライン係数が出力されています。
- １９〜２７行目に３次スプラインの値が出力されています。
ソースコード
（本ルーチンの詳細はe02bbf のマニュアルページを参照）
※本サンプルソースコードは科学技術・統計計算ライブラリである「nAG Fortranライブラリ」のルーチンを呼び出します。
サンプルのコンパイル及び実行方法
| このソースコードをダウンロード |
PROGRAM e02bbfe
! E02BBF Example Program Text
! Mark 23 Release. nAG Copyright 2011.
! .. Use Statements ..
USE nag_library, ONLY : e02bbf, nag_wp
! .. Implicit None Statement ..
IMPLICIT NONE
! .. Parameters ..
INTEGER, PARAMETER :: nin = 5, nout = 6
! .. Local Scalars ..
REAL (KIND=nag_wp) :: a, b, s, x
INTEGER :: ifail, j, m, ncap, ncap7, r
! .. Local Arrays ..
REAL (KIND=nag_wp), ALLOCATABLE :: c(:), lamda(:)
! .. Intrinsic Functions ..
INTRINSIC real
! .. Executable Statements ..
WRITE (nout,*) 'E02BBF Example Program Results'
! Skip heading in data file
READ (nin,*)
READ (nin,*) m
READ (nin,*) ncap
ncap7 = ncap + 7
ALLOCATE (lamda(ncap7),c(ncap7))
READ (nin,*) lamda(1:ncap7)
READ (nin,*) c(1:ncap+3)
a = lamda(4)
b = lamda(ncap+4)
DO r = 1, m
x = (real(m-r,kind=nag_wp)*a+real(r-1,kind=nag_wp)*b)/ &
real(m-1,kind=nag_wp)
ifail = 0
CALL e02bbf(ncap7,lamda,c,x,s,ifail)
IF (r==1) THEN
WRITE (nout,*)
WRITE (nout,*) ' J LAMDA(J) B-spline coefficient (J-2)'
WRITE (nout,*)
DO j = 1, ncap7
IF (j<3 .OR. j>ncap+5) THEN
WRITE (nout,99999) j, lamda(j)
ELSE
WRITE (nout,99999) j, lamda(j), c(j-2)
END IF
END DO
WRITE (nout,*)
WRITE (nout,*) ' R Argument Value of cubic spline'
WRITE (nout,*)
END IF
WRITE (nout,99999) r, x, s
END DO
99999 FORMAT (1X,I3,F14.4,F21.4)
END PROGRAM e02bbfe