程序运行有问题，有关子程序C06LAF

fcode

从 NAG Library 的帮助里可以看到： c06laf 是用于反拉普拉斯变换。

IFAIL = 4
The required accuracy cannot be obtained. It is possible that ALPHAB is less than 
Alternatively, the problem may be especially difficult. Try increasing TFAC, ALPHAB or both.

意思是，需要的精度无法达到。可能是 AlphAB 小于 alpha
试着增加 TFAC ALPHAB 参数，或者两个都增加。

fcode

zjl88915 发表于 2014-2-26 15:45
对啊   c6laf用于拉普拉斯逆变换，但是我按照错误提示对相关参数修改后，问题依然存在。。。。 ...

你手边有 C06laf 的帮助吗？
我建议你认真阅读一下，以便了解函数的局限。设定合适的算例来运算。

我想你的代码里是没有错误的，绝大多数都是 NAG 库里的，输入输出代码很简单，没有问题。

应该还是算例不符合函数库的需求。一般不轻易怀疑函数库的错误。

你的拉氏域函数 Ua 太复杂了，再反变换回来就更复杂了，我实在没有精力帮你分析它。请原谅这些事情我恐怕不能代劳了。我想这恐怕不是一个 Fortran 的代码问题，而是一个学术问题了。