Quellcode-Bibliothek

^© Kompilation durch diese Firma

[Weder Korrektheit noch Funktionsfähigkeit der Software werden zugesichert.]

Datei: int_sin.f90 Sprache: Fortran

!
! Copyright (C) 2005 Intel Corporation.  All rights reserved.
!
! The information and source code contained herein is the exclusive property
! of Intel Corporation and may not be disclosed, examined, or reproduced in
! whole or in part without explicit written authorization from the Company.
!
! This program computes the integral (area under the curve) of a user-supplied
! function over an interval in a stepwise fashion. The interval is split into
! segments, and at each segment position the area of a rectangle is computed
! whose height is the value of sine at that point and the width is the segment
! width.  The areas of the rectangles are then summed.
!
! The process is repeated with smaller and smaller width rectangles, more
! closely approximating the true value.
!
! The source for this program also demonstrates recommended Fortran
! coding practices.
!
program int_sin
implicit none

! Create a value DP that is the "kind" number of a double precision value
! We will use this value in our declarations and constants.
integer, parameter :: DP = kind(0.0D0)

! Declare a named constant for pi, specifying the kind type
real(DP), parameter :: pi = 3.141592653589793238_DP

! Declare interval begin and end
real(DP), parameter :: interval_begin = 0.0_DP
real(DP), parameter :: interval_end   = 2.0_DP * pi

real(DP) :: step, sum, x_i
integer :: N, i, j
real clock_start, clock_finish

write (*,'(A)') " "
write (*,'(A)') " Number of | Computed Integral |"
write (*,'(A)') " Interior Points | |"
call cpu_time (clock_start)

do j=2,26
  write (*,'(A)') "--------------------------------------"
  N = 2**j
  ! Compute stepsize for N-1 internal rectangles
  step = (interval_end - interval_begin) / real(N,DP);

  ! Approximate 1/2 area in first rectangle: f(x0) * (step/2)
  sum = INTEG_FUNC(interval_begin) * (step / 2.0_DP)

  do i=1,N-1
    x_i = real(i,DP) * step
    ! Apply midpoint rule:
    ! Given length = f(x), compute the area of the
    ! rectangle of width step
    sum = sum + (INTEG_FUNC(x_i) * step)
    end do

  ! Add approximate area in last rectangle for f(xN) * (step/2)
  sum = sum + (INTEG_FUNC(interval_end) * (step / 2.0_DP))

  write (*,'(T5,I10,T18,"|",2X,1P,E14.7,T38,"|")') N, sum
  end do

call cpu_time(clock_finish)
write (*,'(A)') "--------------------------------------"
write (*,'(A)') " "
write (*,*) "CPU Time = ",(clock_finish - clock_start), " seconds"

contains

! Function to integrate
real(DP) function INTEG_FUNC (x)
real(DP), intent(IN) :: x

INTEG_FUNC = abs(sin(x))
return
end function INTEG_FUNC

end program int_sin

¤ Dauer der Verarbeitung: 0.1 Sekunden (vorverarbeitet) ¤

Download des Quellennavigators
Download des sprechenden Kalenders
in der Quellcodebibliothek suchen

Haftungshinweis

Die Informationen auf dieser Webseite wurden nach bestem Wissen sorgfältig zusammengestellt. Es wird jedoch weder Vollständigkeit, noch Richtigkeit, noch Qualität der bereit gestellten Informationen zugesichert.

Bemerkung:

Die farbliche Syntaxdarstellung ist noch experimentell.