算数数学教育の暗黒面

黒木玄 (Gen Kuroki)

2018-08-21~2018-08-26

このファイルは次の場所できれいに閲覧できる:

このファイルはJulia言語カーネルの Jupyter notebook である. 自分のパソコンにJulia言語をインストールしたい場合には

を参照せよ. このファイルは JuliaBox でも使用できるかもしれない. このファイル中のJulia言語のコードを理解できれば, Julia言語からSymPyを用いた数式処理や数値計算の結果のプロットの仕方を学ぶことができる.

$ \newcommand\eps{\varepsilon} \newcommand\ds{\displaystyle} \newcommand\Z{{\mathbb Z}} \newcommand\R{{\mathbb R}} \newcommand\C{{\mathbb C}} \newcommand\T{{\mathbb T}} \newcommand\QED{\text{□}} \newcommand\root{\sqrt} \newcommand\bra{\langle} \newcommand\ket{\rangle} \newcommand\d{\partial} \newcommand\sech{\operatorname{sech}} \newcommand\cosec{\operatorname{cosec}} \newcommand\sign{\operatorname{sign}} \newcommand\sinc{\operatorname{sinc}} \newcommand\arctanh{\operatorname{arctanh}} \newcommand\sn{\operatorname{sn}} \newcommand\cn{\operatorname{cn}} \newcommand\cd{\operatorname{cd}} \newcommand\dn{\operatorname{dn}} \newcommand\real{\operatorname{Re}} \newcommand\imag{\operatorname{Im}} \newcommand\Li{\operatorname{Li}} \newcommand\np[1]{:\!#1\!:} \newcommand\PROD{\mathop{\coprod\kern-1.35em\prod}} $

謝辞: このノートの作成にはインターネット上で得た知人との膨大なやりとりの結果理解できたことが含まれている. 1人ひとり名前を挙げることはしないが, 私にたくさんのことを教えてくれた人達にとても感謝している. $\QED$

In [1]:
using Printf
const e = 
linspace(start, stop, length) = range(start, stop=stop, length=length)
using Base64

using Plots
pyplot()
#gr(); ENV["PLOTS_TEST"] = "true"
#clibrary(:colorcet)
#clibrary(:misc)

function pngplot(P...; kwargs...)
    sleep(0.1)
    pngfile = tempname() * ".png"
    savefig(plot(P...; kwargs...), pngfile)
    showimg("image/png", pngfile)
end
pngplot(; kwargs...) = pngplot(plot!(; kwargs...))

showimg(mime, fn; scale="") = open(fn) do f
    base64 = base64encode(f)
    if scale == ""
        display("text/html", """<img src="data:$mime;base64,$base64">""")
    else
        display("text/html", """<img src="data:$mime;base64,$base64" width="$scale">""")
    end
end

using SymPy
#sympy[:init_printing](order="lex") # default
#sympy[:init_printing](order="rev-lex")

using SpecialFunctions
using QuadGK
using Elliptic.Jacobi: cd, sn