Fortran 90における数値積分

Question

Fortran 90では、数学関数をある制限内の1つの変数と数値的に統合したいと考えています。例えば、f（x）= x ** 2を0から10に積分する。私の持つ関数はこれよりも複雑で、積分の限界を数回変えて実行しなければならない。私はインターネット上で 'QUADPACK'ライブラリが私にこれを助けてくれるかもしれないことを知りました。しかし、どのように私はこのコードを私のコードで呼び出すことができるようにこのライブラリをインストールすることができますか？高度な指示に従うことができないため、いくつかの詳細を記入してください。

Answer 1

私は中点法を使ってx^2を統合する単純なプログラムを提供しました。

program integrate 
    implicit none 
    integer,parameter :: cp = selected_real_kind(14) 
    integer,parameter :: N = 1000 
    real(cp),dimension(N) :: f,xc 
    real(cp),dimension(N+1) :: x 
    real(cp) :: s,xmax,xmin,dx 
    integer :: i 
    xmin = 0.0_cp 
    xmax = 10.0_cp 
    dx = (xmax - xmin)/real(N,cp) 
    x = (/(xmin + dx*(i-1),i=1,N+1)/) 
    ! Define x at center 
    do i=1,N 
    xc(i) = x(i) + 0.5_cp*dx 
    enddo 
    ! Define f 
    do i=1,N 
    f(i) = xc(i)**2 
    enddo 
    ! Integrate (Midpoint method) 
    s = 0.0_cp 
    do i=1,N 
    s = s + f(i)*dx 
    enddo 
    write(*,*) 'sum = ',s 
    end program