Julia語言 一維數值積分QuadGK

Julia語言 一維數值積分QuadGK

使用包進行數值積分接下來，讓我們使用Julia的Package

一維數值積分Package是QuadGK.jl

https://github.com/JuliaMath/QuadGK.jl

啟動Julia，按]鍵進入Pkg模式並添加QuadGK。

quadgk（f，a，b）可以積分函數f（x）和a到b的積分區間。

====程式==== 

import Pkg

Pkg.add("QuadGK")

using QuadGK

f(x) = sin(x) + x^2

fsum2 = quadgk(f,-π,π)./(2π)

exact = ((π)^3/3 -(-π)^3/3)/(2π)

println("quadgk $fsum2, exact $exact")

輸出畫面

 Resolving package versions...
  Updating `~/.julia/Project.toml`
 [no changes]
  Updating `~/.julia/Manifest.toml`
 [no changes]
quadgk (3.2898681336964524, 0.0), exact 3.2898681336964524