求教个FORTRAN程序嗯,就是计算如下的式子：以下是我编的程序,在I,J均是正数的时候算出来的结果与MATLAB的计算

求教个FORTRAN程序
嗯,就是计算如下的式子：

以下是我编的程序,在I,J均是正数的时候算出来的结果与MATLAB的计算结果吻合,但是当I,J出现负数的时候结果即发散,试了很多种写法,都有类似的问题,所以来求教.
DOUBLE PRECISION FUNCTION JSC(I,J)
x05COMMON /COEF/ PI
x05DOUBLE PRECISION PI
x05DOUBLE PRECISION K
IF(I.NE.0)THEN
x05 IF(I.NE.J)THEN
x05 IF(I.NE.-J)THEN
K=(1.0D0-(-1.0D0)**(I+J))/(I+J)+(1.0D0-(-1.0D0)**(I-J))/(I-J)
x05 ELSE
x05 K=0.0D0
x05 ENDIF
x05ELSE
x05K=0.0D0
x05ENDIF
x05ELSE
x05K=0.0D0
x05ENDIF
x05JSC=K/2.0D0/PI
RETURN
END