NAG图书馆例行文件
C05ADF(二氧化碳自由基)
1目的
C05ADF通过非线性插值、线性外推法和平分法的组合,在给定区间内定位连续函数的零点。
2规格
整数 |
国际会计师联合会 |
真实(KIND=nag_wp) |
A、 B、每股收益、预计到达时间、F、X |
外部的,外部的 |
F类 |
|
3说明
C05ADF试图获得函数简单零点的近似值给定初始间隔这样的话.使用相同的核心算法C05AZF公司应参考其规范以了解所用方法的详细信息。 近似值到零确定为至少满足以下标准之一:
(i) |
, |
(ii) |
. |
4参考文献
布伦特R P(1973)无导数最小化算法普伦蒂斯-霍尔
5参数
- 1: A–真实(KIND=nag_wp)输入
进入时:,区间的下限。
- 2: B–真实(KIND=nag_wp)输入
进入时:,区间的上界。
约束:.
- 3: EPS–真实(KIND=nag_wp)输入
约束:.
- 4: ETA–真实(种类=nag_wp)输入
进入时:一个值,如果,被接受为零。预计到达时间可以指定为(请参见第7节). - 5: F–REAL(KIND=nag_wp)FUNCTION,由用户提供。外部程序
- 1: XX–真实(KIND=nag_wp)输入
进入时:函数必须求值的点。
F类必须是调用C05ADF的(子)程序USEd的模块子程序,或在其中声明为EXTERNAL。参数表示为输入必须不按此程序进行更改。
- 6: X–真实(KIND=nag_wp)输出
退出时:零的近似值。
- 7: IFAIL–集成输入/输出
-
进入时:国际会计师联合会必须设置为,。如果您不熟悉此参数,请参阅第3.3节有关详细信息,请参见基本介绍。对于检测到错误时可能不适合停止程序执行的环境建议使用。如果错误消息的输出不受欢迎,那么建议使用。否则,如果您不熟悉此参数,建议的值为.当值用于测试国际会计师联合会退出时。
退出时:除非例程检测到错误或已标记警告(请参阅第6节).
6错误指示灯和警告
如果进入或,在当前错误消息单元上输出解释性错误消息(定义如下X04AAF型). 例程检测到的错误或警告:
进入时, | , |
或 | , |
或 | . |
计算中要求的精确度太高;也就是说,零点的定位至少达到了相对精度,其中是机器精度,但中描述的退出条件第3节不满意。C05ADF继续超过这一点是不安全的,但X(X)返回的是零的精确近似值。
符号的变化已确定发生在由最终值定义的点附近X(X)然而,有证据表明,这种符号转换对应于.
7准确度
准确度取决于每股收益和预计到达时间。如果需要全机器精度,则它们可能设置得很小,从而导致退出,尽管这可能需要更多的迭代,而不是更低的精度。建议您设置并使用每股收益控制精度,除非您对对于的值接近零点。
C05ADF所花费的时间主要取决于评估所花费的时间F类(请参见第5节). 如果确定相对长度小于的间隔很重要包含零,或者如果F类评估成本高F类将受到限制,然后使用C05AZF公司建议使用。的使用C05AZF公司当要解决的问题的结构不允许简单F类待写:反向通信设施C05AZF公司比直接沟通更灵活F类C05ADF要求。 9示例
此示例计算的近似值为间隔内使用的公差为.
9.1程序文本
程序文本(c05adfe.f90)
9.2程序数据
没有。
9.3项目结果
程序结果(c05adfe.r)