我的下面这段程序 出现了如题的错误，将Ka与Kc取绝对值后还是出现这个错误，求高手指点一下呀。
      for i=1 to N
        fmA=(1.13-0.09*af/cf+(-0.54+0.89/(0.2+af/cf))*(af/B)^2+(0.5-1/(0.65+af/cf)+14*(1-af/cf)^24)*(af/B)^4)/sqrt(1+1.464*(af/cf)^1.65)
        fbA=(1+(-1.22-0.12*af/cf)*af/B+(0.55-1.05*(af/cf)^0.75+0.47*(af/cf)^1.5)*(af/B)^2)*fmA
        fmB=((1.1+0.35*(af/B)^2)*sqrt(af/cf))*fmA
        fbB=(1-0.34*af/B-0.11*af^2/(cf*B))*fmB
        note 确定a向裂纹的应力强度因子 Ka
        Ka=abs(sqrt(pi()*af)*(σm*fmA+σB*fbA))
        note 确定c向裂纹的应力强度因子 Kc
        Kc=abs(sqrt(pi()*af)*(σm*fmB+σB*fbB))
        Kap=Ka^m                                       &&错误：“***或^定义域错误！”
        Kcp=0.9^m*Kc^m
        af=af+A*Kap
        cf=cf+A*Kcp
      endfor
       &&对负数只能取整数次幂。 使用幂操作符时，分数不能为负数，请采用整数幂计算。

当CA=0 时，M<0 时 ，将提示：“***或^定义域错误！”