[問題] 數值太大導致run time error

[SYlin (喔)]

我有個問題

b=10E5 w=0.5

在跑到這一行的時候

s=((b*w)**2.0*exp(b*w))/(1+exp(b*w))**2.0

就會出現 run time error

請問有什麼辦法可以解決嗎?

希望各位高手幫幫忙

感恩~

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.116.91.131
推 egggg:改寫成s=((b*w)**2.0*exp(-b*w))/(1+exp(-b*w))**2.0
推 sjgau:使用 REAL*8
推 egggg:改成real*8還是有危險，不過這數根本上是個0...
→ SYlin:已經是REAL*8了
→ SYlin:感謝樓上兩位  我後來取LOG算  再用EXP  就可以算了
→ SYlin:不過不知道這樣精確度會不會比較差
推 Gwaewluin:應該是醫定變差吧，又是log又是exp的
→ Gwaewluin:      一