24小时热门版块排行榜    

Znn3bq.jpeg
查看: 922  |  回复: 8

zouji1983

金虫 (正式写手)

[求助] Fortran error求助

Fortran error: run-time error M6201:MATH-logOMAIN error
我的Log表达式是这样的:
f(x)=1.72992961d9*log(abs(x))-1.72992961d9*log(abs(1-x))-1.182988286d10*(1-2*x)+2.75714286d9*(6*x-6*x**2-1)+3.083731429d9*(24*x**2-10*x-16*x**3+1)-1d11*(xc-phi1*x)/phi2 + 9.229804981d10
X一定不是负数或者0啊,还有什么原因啊?谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiwei43

至尊木虫 (正式写手)

一定是数据超界了,你简单手算一下,看看哪部分超界了。
2楼2012-11-27 11:10:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zouji1983

金虫 (正式写手)

引用回帖:
2楼: Originally posted by weiwei43 at 2012-11-27 11:10:43
一定是数据超界了,你简单手算一下,看看哪部分超界了。

谢谢!不好意思,想问一下,数据超界了是什么意思?
3楼2012-11-27 16:38:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zouji1983

金虫 (正式写手)

引用回帖:
2楼: Originally posted by weiwei43 at 2012-11-27 11:10:43
一定是数据超界了,你简单手算一下,看看哪部分超界了。

怎么检查呢?
4楼2012-11-27 17:01:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiwei43

至尊木虫 (正式写手)

引用回帖:
4楼: Originally posted by zouji1983 at 2012-11-27 17:01:59
怎么检查呢?...

整数,实数都是有范围的,超过这个范围就是超界。
一定是你那个计算结果太大了,或者太小了超界了。
你可以把某个值(最好是变量极值)带入表达式,手算一下是不是超界了。
5楼2012-11-28 08:39:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (著名写手)

感谢参与,应助指数 +1
本帖仅楼主可见
6楼2012-11-28 15:35:32
已阅   申请程序强帖   回复此楼   编辑   查看我的主页

zouji1983

金虫 (正式写手)

引用回帖:
6楼: Originally posted by snooper007 at 2012-11-28 15:35:32
你用iIF语句判别啊

谢谢大家的帮助!我的变量X的范围是0.2~0.999,by the way, 怎么用if来判别呢?
7楼2012-11-28 18:01:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (著名写手)

本帖仅楼主可见
8楼2012-11-28 23:33:02
已阅   申请程序强帖   回复此楼   编辑   查看我的主页

zouji1983

金虫 (正式写手)

引用回帖:
8楼: Originally posted by snooper007 at 2012-11-28 23:33:02
变量X的范围是0.2~0.999, 有可能使得log(abs(1-x)出现domain  error
你用双精度试验一下...

但实际上我定义的X就是双精度的实数啊
9楼2012-11-29 21:47:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zouji1983 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 评审有感 +7 popular289 2026-05-18 12/600 2026-05-19 05:02 by redmoonzpc
[教师之家] 上海大学实验技术岗位非升即走 +9 嘻嘻哈哈乐呵呵 2026-05-15 9/450 2026-05-19 01:47 by jurkat.1640
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[基金申请] 今年审到国自然15份,谈谈感受 +16 国自然国社科中 2026-05-17 16/800 2026-05-18 14:58 by gy116024
[硕博家园] 考博自荐 +5 科研狗111 2026-05-13 6/300 2026-05-18 11:22 by 糊糊涂涂好
[基金申请] 青C资助名额大幅增加! +12 西葫芦炒鸡蛋 2026-05-13 16/800 2026-05-18 10:02 by Equinoxhua
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 5/250 2026-05-17 18:39 by Equinoxhua
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 7hlccowb3h 2026-05-15 4/200 2026-05-17 07:46 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 8/400 2026-05-17 07:26 by 11n4dfd8yn
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:11 by 11n4dfd8yn
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 7/350 2026-05-17 06:55 by 11n4dfd8yn
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:37 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:45 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[基金申请] 请问大佬b0816评完了吗 +4 市民华南虎 2026-05-12 8/400 2026-05-16 19:54 by Equinoxhua
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
请填处理意见