home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l292 / 1.ddi / INTEGRDE.FOR < prev    next >
Encoding:
Text File  |  1989-10-10  |  1.0 KB  |  40 lines
  1.  
  2.  
  3.       INCLUDE 'STDHDR.FOR'
  4.       REAL  y1( 0: maxv)     !Array of function values
  5.       REAL  xl1              !lower limit of integration
  6.       REAL  xh1              !upper limit of integration
  7.       REAL  ans1             !value of integral
  8.       INTEGER  xxl1,xxh1,i1,numseg
  9.       REAL  sp,i
  10.  
  11.  
  12.       DO i1 = 0, 10          !sets up a vector which is a piecewise
  13.                              !approximation of the curve y = e^x
  14.         i = i1/10.0          !across the interval 0 - 1}
  15.          y1(i1) = exp(i)
  16.       END DO
  17.       xl1 = 0.0
  18.       xxl1 = 0
  19.       xxh1 = 10
  20.       xh1 = 1.0
  21.       sp = 1.0/10.0
  22.       CALL IntegrateVector(y1,sp,xxl1,xxh1,ans1)
  23.       PRINT *, 'Integrate Vector'
  24.       PRINT *, 'Area under the curve = ', ans1
  25.       PRINT *
  26.       PRINT *, 'Integrate Function'
  27.       DO i1 = 1, 10
  28.         numseg = 2*i1
  29.         CALL IntegrateFunction(numseg,xl1,xh1,ans1)
  30.         PRINT *, 'Number of segments = ',numseg
  31.         PRINT *, '  Area under the curve = ',ans1
  32.       END DO
  33.       END
  34.  
  35.  
  36.  
  37.  
  38.  
  39.  
  40.