大侠救命！g77编译后赋值语句错误！

平谷六号

请看这一段程序：
            write(*,*) 'x0 =', x0
            X1=x0-((DFCPHI(X0, THETA)-AIPVAL)/pfcphi(x0,theta))
            WRITE(*,*)'X1 =',X1
            write(*,*) 'DFCPHI(X0, THETA)-AIPVAL =',
     #        DFCPHI(X0, THETA)-AIPVAL
            write(*,*) 'DFCPHI(X0, THETA) =',DFCPHI(X0, THETA)
            write(*,*) 'PFCPHI(X0, THETA) =',PFCPHI(X0, THETA)
            write(*,*)
     #        'x0-((DFCPHI(X0, THETA)-AIPVAL)/pfcphi(x0,theta))=',
     #        x0-((DFCPHI(X0, THETA)-AIPVAL)/pfcphi(x0,theta))
显然最后三行输出的表达式就是第三行赋值语句的右侧，因此第四行的输出值应与最后三行相等，这之间也没有再对X0、X1和THETA显式或隐式地赋值，可是请看输出结果：
x0 =  0.00852136688
X1 = -1.80225983
DFCPHI(X0, THETA)-AIPVAL = -0.00117079928
DFCPHI(X0, THETA) =  0.000185418379
PFCPHI(X0, THETA) =  0.180912063
x0-((DFCPHI(X0, THETA)-AIPVAL)/pfcphi(x0,theta)) =  0.0149930154
两个输出语句得到两个毫不相干的结果！这是怎么回事？我正在算一个很急的题目，有熟悉这个的大侠请拉兄弟一把，谢谢了！
我的操作系统是redhat 9.0，平台Xeon 2.4 GHz x 2，主板E7500CW2，1G Kingston ECC Reg内存。