奥村先生の Python のページで，出生性比のグラフの描き方が紹介されていたhttps://oku.edu.mie-u.ac.jp/~okumura/python/sexratio.html Julia でもやってみた。 CSV ファイルをダウンロードして保存した段階で，エディタで開いて，utf-8 ファイルにした（そ…

折れ線グラフと棒グラフ--どっちが見やすい? この場合なら，私は 棒グラフのほうを選ぶ # https://covid19.mhlw.go.jp/public/opendata/newly_confirmed_cases_daily.csv# 以下では，ダウンロードしたデータを参照 using DataFrames, CSV, Datesdf = CSV.rea…

あいも変わらず，「2直線の交点座標を求めるプログラム」がよく参照されているが，わけわからん。 そこで，Julia の SymPy で書いたものも掲載しておく。これが「わけわからん」という人もいるだろうが。 using SymPy @syms a b c d x function func(a, b, c…

SymPy/Julia https://www.math.csi.cuny.edu/~verzani/tmp/julia/symbolic.html 5. 数式の図示 数式は，Julia へ戻ることなく，直接プロットすることができる。 f(x) = sin(x^4) - cos(x^4) を [0,1] の範囲で描くには以下のようにする。 using Plots gr(lab…

SymPy/Julia https://www.math.csi.cuny.edu/~verzani/tmp/julia/symbolic.html 1. Julia で数式処理 コンピュータによる計算処理で，３ 通りくらいのアプローチがある。 Mathenatica，Maple, Sage に代表される数式処理 MATLAB, R による数値処理 Python, H…

「ベクトルの内積関連の問題を Julia の SymPy で解く」の問題 2 の (2) で，「a = [1, 1], b = [1-√3, 1+√3] の b の定義で直接 sqrt(3) を使うと，数値解を求めようとするので，それを避けるために変数 x を用いて sqrt(x) とする。」などと書いてしまった…

1. 高校数学の「高次方程式・恒等式」関連の問題をPythonで解く https://qiita.com/code0327/items/66efb7ac56cdea5b0b0d Julia の SymPy でやってみる。 using SymPy @syms x f = x^3 - 1; # (1) f.(x => 2) |> string "7" # (2) f.(x => 1) |> string "0" …

高校数学の「ベクトルの内積」関連の問題をPythonで解く https://qiita.com/code0327/items/b8614649757110c654b6 「 高校数学の「ベクトルの内積」関連の問題をPythonで解く 」を参考に、値を最後に代入してみた。https://qiita.com/mrrclb48z/items/2d6779…

クレジットカードのチェックディジット クレジットカード番号には規則性がある。 一見して区別はできないが，ある計算処理をした結果を見れば，それが真正な番号であるかそうでないかはすぐにわかる。 数字の各桁が特定の規則に従う場合，それぞれの桁につい…

Nov 19, 2021 Julia v1.6.4 M1 チップの場合，私がインストールしているパッケージでは PNGFiles だけがプリコンパイルに失敗しています。 Nov 15, 2021 Julia v1.7.0-rc3 (Intel or Rosetta / M-series Processor) M1 チップへの対応は，まだまだ問題が残っ…

今年 2021 年の次の 3 個の素数年は, 2027, 2029, 2039 年である。 2027 で割ると余りが 1012029 で割ると余りが 1032039 で割ると余りが 107 となるような最小の正整数を求めよ 答えは，ずっと下へスクロール function chinese_remainder_theorem(a1, r1, a…

https://www.imojp.org/archive/mo2021/jmo2021/problems/jmo31yq.html Q1. 互いに素な正の整数 m, n が m + n = 90 をみたすとき，積 mn としてありうる，最大の値を求めよ。 value = 0 mn = () for m = 1:45 n = 90 - m gcd(m, n) == 1 && (value = max(va…

2020 年日本数学オリンピック予選 問題 1，3，4，5 https://www.imojp.org/archive/mo2020/jmo2020/problems/jmo30yq.html Q1. 千の位と十の位が 2 であるような 4 桁の正の整数のうち，7 の倍数はいくつあるか。 count = 0 for i = 0:100:900, j = 0:9 (202…

2020 年日本数学オリンピック予選 問題7 https://www.imojp.org/archive/mo2020/jmo2020/problems/jmo30yq.html Q7. 2×1010 のマス目の各マスに 1 以上 5 以下の整数を 1 つずつ書き込む。辺を共有して隣りあうどの 2 マスについても書き込まれた数の差が 2 …

Julia 1.6.3R 4.1.1macOS Monterey 12.0.1M1 chip Mac mini 2020 で，RCall が使えなくなっていたが，R 4.1.2 にバージョンアップしたら使えるようになった。 julia> using RCall[ Info: Precompiling RCall [6f49c342-dc21-5d91-9882-a32aef131414] julia> …

https://www.imojp.org/archive/mo2020/jmo2020/problems/jmo30yq.html Q6. 平面上に 3 つの正方形があり，図のようにそれぞれ 4 つの頂点のうち 2 つの頂点を他の正方形と共有している。ここで，最も小さい正方形の対角線を延長した直線は最も大きい正方形…

M1 チップ搭載の Mac minimacOs Monterey version 12.0.1Julia 1.6.3で，using SymPy しようとすると， /Users/foo/.julia/packages/SymPy/jTadH/src/numbers.jlの 256 行目からの式でエラーがあるので使えないよということになった。 件の式は， sympy_core…

1. Julia の小ネタ--043 データフレームを clipboard にコピーする clipboard(df) とすれば，df が clipboard にコピーされるのではあるが，期待するようなタブ区切りのデータではない。 以下のような単純なデータフレーム df を clipboard(df) してみると，…

1. Julia の小ネタ--042 clipboard からデータフレームを読み込む 「Julia の小ネタ--038 文字列からデータフレームを読み込む」に似ているが，clipboard の内容からデータフレームを作る例を挙げておく。 IOBuffer( ) に，クリップボードから読み込んだ文字…

clipboard へ出力 clipboard(str) とすると，str の内容がクリップボードにコピーされる。そのあと，クリップボードの内容をペースト(command + v) できる。 以下を実行後， str = "0123456789\nabcdefghij\n" clipboard(str) 別のアプリの入力場面で，comma…

1. clipboard から入力 単に clipboard() とすると，その時点でのクリップボードの内容を取り出すことができる。 文字列として利用するなら，そのままでよい。 以下の 2行をコピーした後，実行すると，2 行がそっくりそのまま表示される。ただし，実行の若干…