パーセクの定義に $1$ 秒角の $\tan$ （あるいは $\tan$ の逆数である $\cot$）がいらなくなって，$\displaystyle 1 \mbox{pc} = \frac{\mbox{au}}{\tan 1^{\prime\prime}} = \mbox{au}\ \cot 1^{\prime\prime}$ ではなく，新定義 (IAU 2015 Resolution B2) では
$$1 \mbox{pc} \equiv \frac{\mbox{au}}{1^{\prime\prime}} = \frac{648000}{\pi} \mbox{au}$$
となったという話。
同等の対象を
plot2d(), draw2d())でグラフにしてみるシリーズ。 Continue reading
初等関数のグラフを描く際，$y=\tan x$ や $y=1/x$ などの不連続点を無造作に線でつなぐケースがあったので，不連続点を無造作につなげないようなグラフの描き方をまとめてみた。 Continue reading
理工系の数学 B の授業で，楕円の周長や面積，回転楕円体の表面積や体積を求めているので。
Python の SymPy Plotting Backends を使って，パラメータを変えて複数の曲線を一気にグラフにする例。グラフを描く際のオプションの設定についても2通りの例をまとめた。 Continue reading
sympy.abc のコードには，
Caveats ======= 1. As of the time of writing this, the names ``O``, ``S``, ``I``, ``N``, ``E``, and ``Q`` are colliding with names defined in SymPy. If you import them from both ``sympy.abc`` and ``sympy``, the second import will "win". This is an issue only for * imports, which should only be used for short-lived code such as interactive sessions and throwaway scripts that do not survive until the next SymPy upgrade, where ``sympy`` may contain a different set of names.
とある。また，sympy.abc にはギリシャ文字の大文字が定義されていない。したがって，今後は以下のようにしたらどうか，というメモ。 Continue reading
SymPy Plotting Backends で”少し込み入った関数”を plot() するとき，思いもかけずにエラーとなった場合の対応。
以下の Issue で回答があった件：
「高さ h からの斜方投射の最大到達距離を求める準備」のページを参照。
最大水平到達距離となる打ち出し角度 $\theta$ を，陰関数定理を使い，連立方程式の形にして SymPy や Maxima を使って求めてみる。 Continue reading