ソフトウェアトーク理工サイド交流祭に向けて作った動画です。こだわってるといつまでも投稿できなさそうなので、思い切って投稿します。祭りを知るのが遅すぎたのじゃ…。

5分くらいの動画を想定して、内容を削りに削りました。削った内容については、注釈や補足資料などに詰まっています。

間違ってるところがあれば(優しく)指摘いただけると助かります。Maximaについては素人なので、もっと良い使い方があるかもしれません。

今回使うMaximaのダウンロード
https://maxima.sourceforge.io/download.html

ファインマン物理学(英語版)
https://www.feynmanlectures.caltech.edu/ 

ロジスティック写像
https://ja.wolframalpha.com/input?i=ロジスティック写像

Maximaコピペ用

[m,k,r0,v0,t0,te,dt]:[0.04,1.4*10^-4,0,30,0,50,0.1];
F(x,v,t):=(-k*v^2)/m;
rk_sol:rk([v,F(r,v,t)],[r,v],[r0,v0],[t,t0,te,dt])$
plot2d([discrete, makelist([first(i),second(i)],i,rk_sol)]);