大家好,目前我使用的是XE2013 with VS2012
目前碰到了一個問題是使用imsl計算bessel function
當回傳數值大於320次方時程式直接shut down
舉例來說
MA(I,J)=DBSI0(LN1(I)*RSD)/(LN1(I)*DBSI1(LN1(I)*RSD))
其中DBSI0與DBSI1是imsl的double precision 的 bessel function指令
當LN1(I)*RSD的值大於700時,會導致DBSI0與DBSI1回傳的值分別都大於320次方
實際上在數學運算時MA(I,J)算出來差不多10的-2到-3之間,但卻會因為分母與分子回傳
的值大於320次方而導致程式終止
有沒有先進能提供方法讓數值超過320次方依然能運算?
或是有甚麼解決方法,先謝謝各位