連続関数の根

ホーム > 製品 > nAG数値計算ライブラリ > nAG Fortranライブラリ > サンプルソースコード集 > 連続関数の根

Keyword: 連続関数の根

概要

本サンプルは連続関数の根を求めるFortranによるサンプルプログラムです。本サンプルは許容誤差を1.0e-5とする、区間[1,0]の範囲内の以下の連続関数の根の近似を求めて出力します。非線形補間法、外挿法、二分法を組み合わせた手法が使用されています。

連続関数の根のデータ　

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

出力結果

（本ルーチンの詳細はc05adf のマニュアルページを参照）

この出力例をダウンロード

 C05ADF Example Program Results

 Zero at X =     0.56714

３行目に連続関数の根が出力されています。

ソースコード

（本ルーチンの詳細はc05adf のマニュアルページを参照）

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

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

!   C05ADF Example Program Text
!   Mark 23 Release. nAG Copyright 2011.

    MODULE c05adfe_mod

!      C05ADF Example Program Module:
!             Parameters and User-defined Routines

!      .. Use Statements ..
       USE nag_library, ONLY : nag_wp
!      .. Implicit None Statement ..
       IMPLICIT NONE
!      .. Parameters ..
       INTEGER, PARAMETER                  :: nout = 6
    CONTAINS
       FUNCTION f(x)

!         .. Implicit None Statement ..
          IMPLICIT NONE
!         .. Function Return Value ..
          REAL (KIND=nag_wp)                  :: f
!         .. Scalar Arguments ..
          REAL (KIND=nag_wp), INTENT (IN)     :: x
!         .. Intrinsic Functions ..
          INTRINSIC                              exp
!         .. Executable Statements ..
          f = exp(-x) - x

          RETURN

       END FUNCTION f
    END MODULE c05adfe_mod
    PROGRAM c05adfe

!      C05ADF Example Main Program

!      .. Use Statements ..
       USE nag_library, ONLY : c05adf, nag_wp
       USE c05adfe_mod, ONLY : f, nout
!      .. Implicit None Statement ..
       IMPLICIT NONE
!      .. Local Scalars ..
       REAL (KIND=nag_wp)                  :: a, b, eps, eta, x
       INTEGER                             :: ifail
!      .. Executable Statements ..
       WRITE (nout,*) 'C05ADF Example Program Results'

       a = 0.0E0_nag_wp
       b = 1.0E0_nag_wp
       eps = 1.0E-5_nag_wp
       eta = 0.0E0_nag_wp

       ifail = -1
       CALL c05adf(a,b,eps,eta,f,x,ifail)

       WRITE (nout,*)

       SELECT CASE (ifail)
       CASE (0)
          WRITE (nout,99999) 'Zero at X =', x
       CASE (2,3)
          WRITE (nout,99999) 'Final point = ', x
       END SELECT

99999  FORMAT (1X,A,F12.5)
    END PROGRAM c05adfe

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

概要

出力結果

ソースコード