忘却の微分方程式(175) Maxima、{draw}、クーロンの法則でベクトル場を描く

Joseph Halfmoon

前回、3次元のベクトル場とスカラー場をお惚け老人なりにプロット。さて何か意味のあるものを描きて~ということになりました。そこで今回は、みんな知ってるクーロンの法則で描いてみたいと思います。想定は「原点に電子1個あり、その周りに別な電子1個を置いたときに働く力」です。まあ計算するまでもなく1点から周りにドヒャーと。 “忘却の微分方程式(175) Maxima、{draw}、クーロンの法則でベクトル場を描く” の続きを読む

忘却の微分方程式(174) Maxima、{draw}、3次元ベクトル場とスカラー場表示

3DplotEC
Joseph Halfmoon

このところvectパッケージのベクトル解析演算子共を試用。しかし、お惚け老人の「プロット能力の制約」により、その適用はいつも2次元空間でした。2次元のベクトル空間への演算子の適用結果を3次元プロットするのならば簡単。でも、本来3次元空間に対して適用するのが筋ってもんじゃありませんか。今回は少しジタバタしてみましたぞ。
“忘却の微分方程式(174) Maxima、{draw}、3次元ベクトル場とスカラー場表示” の続きを読む

忘却の微分方程式(173) Maxima、{vect}、ベクトル解析、ラプラシアン

plotLaplacianEC
Joseph Halfmoon

graddivcurl(rot)と練習してきて今回は Laplacian です。過去回はスカラー場からベクトルが飛び出てきたり、あるいはその逆であったりしました。今回は地味。スカラー場をLaplacianすると出てくるのはスカラー。でも計算の奥底にはベクトルが息づいている?凸とか凹とか調べるのにお馴染みの操作です。
“忘却の微分方程式(173) Maxima、{vect}、ベクトル解析、ラプラシアン” の続きを読む

忘却の微分方程式(172) Maxima、{vect}、ベクトル解析、ローテーション

rotation2D
Joseph Halfmoon

前回はダイバージェンス(div)、今回はローテーションです。遥かな太古、お惚け老人が学校でベクトル解析を習ったとき、rot (rotation)というお名前で習いました。rotでなく、curl とお書きになる流儀もあるみたい。Maxima様はどうも curl 流であるようです。rotと書いたらエラーになりましたぜ。 “忘却の微分方程式(172) Maxima、{vect}、ベクトル解析、ローテーション” の続きを読む

忘却の微分方程式(171) Maxima、{vect}、ベクトル解析、ダイバージェンス

DivE_2D_EC
Joseph Halfmoon

前回は、スカラー場からベクトル(勾配)を求めるグラディエント(grad)を計算してみました。プロットもなんとかなる感じ?今回はベクトル場から発散(ダイバージェンス)を求めるdiv()関数を試用してみます。こんこんとベクトルが湧いてくるのよ?でも、お惚け老人には3次元空間ベクトル場の表現はムツカシイのでいつもの2次元ね。 “忘却の微分方程式(171) Maxima、{vect}、ベクトル解析、ダイバージェンス” の続きを読む

忘却の微分方程式(170) Maxima、{vect}、ベクトル解析、gradのプロット

Gradient_EC
Joseph Halfmoon

前回、vectパッケージでgradient(勾配)の計算済。しかしプロットしてなかったです。今回はまず「スカラー場」の関数を定義してプロット、そしてそのスカラー場をgrad()関数に食わせてベクトルを計算、そしてベクトル場をプロットとな。ただしお惚け老人には3次元空間ベクトル場の表現はムツカシイので2次元ね。
“忘却の微分方程式(170) Maxima、{vect}、ベクトル解析、gradのプロット” の続きを読む

忘却の微分方程式(169) Maxima、{vect}、ベクトル解析パッケージ再び

cross
Joseph Halfmoon

過去回にて「ちょっと触って」みたものの、消化不良のまま打っちゃっていたパッケージに今回戻りたいと思います。vectパッケージとな。ベクトル解析用のパッケージです。このパッケージをば極めれば、グラディエントにダイバージェンス、カールにラプラシアンと昔苦しめられた記憶の者どもを自在に使役できるようになるのだとか?ホントか?
“忘却の微分方程式(169) Maxima、{vect}、ベクトル解析パッケージ再び” の続きを読む

忘却の微分方程式(168) Maxima、for、do、next、thru他 繰り返し

for_do_002
Joseph Halfmoon

前回は Maxima様の縁の下に隠れているLisp(Common Lisp)を呼び出すためにMaxima様の「Program Flow関係」に踏み込んでしまいました。毒を食らわば皿まで、ということで(何が毒だ)、今回はMaxima中でループを構成する方法を練習してみます。forとかdoとか「ありがちな奴ら」が大挙登場。
“忘却の微分方程式(168) Maxima、for、do、next、thru他 繰り返し” の続きを読む

忘却の微分方程式(167) Maxima、lisp呼び出し、コマケー話に躓くんだ、これが。

lispMode
Joseph Halfmoon

最近、別シリーズにてマイコン用のuLispをラズパイPico2上で走らせてます。COMMON LISPのサブセット的な。一方、Maxima様もCOMMON LISP上で「走っている」ハズ。もしやMaxima様の上でちょこっとLispの練習できるんじゃあ~りませんか、と思った次第。でもね、コマケー所に躓くんだ。これが。

“忘却の微分方程式(167) Maxima、lisp呼び出し、コマケー話に躓くんだ、これが。” の続きを読む

忘却の微分方程式(166) Maxima、{ezunits}、単位と変換、次元解析、物理定数

unit004
Joseph Halfmoon

頭の抽象度が低い年寄っす。年寄になる前からそうなので、数字には何か単位がついてないと不安を感じます。単位がついていたからといって分かったというわけでもないんだけれども。でもま、そんな抽象度の低い老人にもMaxima様のご配慮は行き届いております。ezunitsパッケージとな。直接呼ぶ必要もないのだけれども。 “忘却の微分方程式(166) Maxima、{ezunits}、単位と変換、次元解析、物理定数” の続きを読む

忘却の微分方程式(165) Maxima、{fractals}、ヒルベルト曲線、コッホ曲線

snowmaps
Joseph Halfmoon

前回まで dynamics パッケージの紡ぎだすカオスの中にフラクタルを観察してまいりました。しかし、dynamicsパッケージだけではなかったです。その名もズバリの fractals パッケージというものが存在。いままで見てきたようなフラクタル図形も描画できるけれども、ちょっと立ち位置が違うみたい?なんだそれ。

“忘却の微分方程式(165) Maxima、{fractals}、ヒルベルト曲線、コッホ曲線” の続きを読む

忘却の微分方程式(164) Maxima、{dynamics}、シダの葉にみる反復関数系

Barnsley_fern
Joseph Halfmoon

前回は連立漸化式の発展を観察。今回は反復関数系です。IFS: Iterated Function System。2次元の点をこねくり回す?のだけれども、chaos game同様、そこにランダムな介入があると。そんな乱数な?といってあら不思議。繰り返していくと何やら図形が見えてきます。その代表がバーンズリー先生のシダね。
“忘却の微分方程式(164) Maxima、{dynamics}、シダの葉にみる反復関数系” の続きを読む

忘却の微分方程式(163) Maxima、{dynamics}、こんどは連立漸化式の発展だ

evo2Dexample0
Joseph Halfmoon

第161回でevolution関数を使い「漸化式の発展」の様子を観察しました。今回はその2次元バージョン? evolution2D関数です。発展の様子を観察する対象は「連立漸化式」じゃないかと。ムツカシそうだよ。でもカオスでフラクタルな処理例はちょっとさらに上を行くのさ。ホントか?よくわからんが。

“忘却の微分方程式(163) Maxima、{dynamics}、こんどは連立漸化式の発展だ” の続きを読む

忘却の微分方程式(162) Maxima、{dynamics}、階段図で探る再帰数列のゆくえ

staircase384
Joseph Halfmoon

前々回はorbits関数で「軌道図」を描き、前回は漸化式の「発展」をevolution関数で観察しました。そして今回は staircase関数で「階段図」です。前回、前々回同様、解析の対象は「再帰的に計算される漸化式」ですが表現の仕方で随分異なる様相が見えるような気がします。気がするだけですケド。。。

“忘却の微分方程式(162) Maxima、{dynamics}、階段図で探る再帰数列のゆくえ” の続きを読む