compare.f

一维热传导方程的求解

      PROGRAM COMPARE
      
      IMPLICIT NONE
      
      INTEGER           I, I1, I2
      DOUBLE PRECISION  EXPL,IMPL
      
      
      OPEN(90,FILE='EXPL.OUT')
      OPEN(91,FILE='IMPL.OUT')
      OPEN(50,FILE='DIFF.OUT')      
      
      DO I = 1, (100*21)
      READ(90,*) I1, EXPL
      READ(91,*) I2, IMPL
      
      WRITE(50,'(I4,1X,F14.6)') I1, (EXPL-IMPL)
      
      ENDDO
      
      END