大家好,目前我使用的是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次方依然能运算?
或是有什么解决方法,先谢谢各位