2023年12月25日にJulia v1.10.0がリリースされました。おめでとうございます。
ちょうどJulia Advent Calendar 2023 https://qiita.com/advent-calendar/2023/julia に 積分のパッケージの紹介 https://zenn.dev/ohno/articles/440234fbb2adec がありました。
そのうちの一つ、QuadGKパッケージ https://juliamath.github.io/QuadGK.jl/stable/ を使って 『岩波数学公式I』のp.240の例の数値積分を精度保証しようとおもいます。 永らくゼロだとされ、1988年に訂正され、2018年にTwitterで少し話題になっていた積分です。
| ### A Pluto.jl notebook ### | |
| # v0.14.2 | |
| using Markdown | |
| using InteractiveUtils | |
| # ╔═╡ 9dcadaad-5e5e-494f-931d-1ab9269ed8c1 | |
| using Plots, LinearAlgebra | |
| # ╔═╡ 077d65e8-a167-4527-9bb8-8fa9eb614c46 |
| using UnicodePlots | |
| using OrdinaryDiffEq | |
| using Parameters | |
| struct Numerov <: OrdinaryDiffEqAlgorithm end | |
| mutable struct NumerovCache <: OrdinaryDiffEq.OrdinaryDiffEqMutableCache | |
| uprev # last u value | |
| g | |
| gprev |
| #= | |
| This julia script converts fortran 90 code into julia. | |
| It uses naive regex replacements to do as much as possible, | |
| but the output WILL need further cleanup. | |
| Known conversion problems such as GOTO are commented and marked with FIXME | |
| Most variable declaration lines are entirely deleted, which may or | |
| may not be useful. |
| ! MT19937, Mersenne Twister Random Number Generator ([0,1) Real Number) | |
| ! in Fortran90 | |
| ! | |
| ! Usage: | |
| ! 1) When you use mt19937, add the sentence "use mt19937" | |
| ! above the implicit sentence. | |
| ! 2) To set an initial seed, call sgrnd(seed). (The "seed" is an integer.) | |
| ! If you do not call this, the seed is 4357. | |
| ! 3) Use the function grnd(). | |
| ! (Do not declare "real(8) :: grnd".) |