area.for

著名的berkeley大学开发的结构动力非线性分析程序

C*****
      REAL*8 FUNCTION AREA (Y1,Y2,DX)
C
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
C
      IF (Y1*Y2.GE.0.0) THEN
         AREA=0.5*(Y1+Y2)*DX
      ELSE
         TEMP=DABS(Y1)+DABS(Y2)
         AREA=0.5*Y1*DX*DABS(Y1)/TEMP
         AREA=AREA+0.5*Y2*DX*DABS(Y2)/TEMP
      ENDIF
      RETURN
      END