対象 数値計算を使って勉強している物理,化学系の学生向け はじめに 物理現象を記述する微分方程式を数値的に解く,みたいなのは理系専攻のカリキュラムなら,どこもやっていると思います. とりあえず差分化してコードに直せば,解を出してくれて,グラフに直せばそのイメージが掴める,とても便利な方法です. 一方で面倒なこともあります. 一つがコードを書く手間です. 数値計算の講義ではCが主流かと思いますが,C言語は手軽に扱うには向きません.デバッグに時間もかかります. また,刻み幅の調整も問題です.この値が大きすぎると正しい解が得られません,かといって,小さすぎると計算時間がかかりすぎます. この二つの面倒さに対し,自分なりに解決を行いましたので,紹介することにしました pythonとscipy.integrate.solve_ivp まず,言語としてPythonを採用します.pythonといえば遅い
![scipy.integrate.solve_ivpで微分方程式をたくさん解いてみた - Qiita](https://cdn-ak-scissors.b.st-hatena.com/image/square/d2a30452f29c48d4a86cb999bb9811369f25b278/height=288;version=1;width=512/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn.qiita.com%252Fassets%252Fpublic%252Fadvent-calendar-ogp-background-7940cd1c8db80a7ec40711d90f43539e.jpg%3Fixlib%3Drb-4.0.0%26w%3D1200%26mark64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTk3MiZoPTM3OCZ0eHQ9c2NpcHkuaW50ZWdyYXRlLnNvbHZlX2l2cCVFMyU4MSVBNyVFNSVCRSVBRSVFNSU4OCU4NiVFNiU5NiVCOSVFNyVBOCU4QiVFNSVCQyU4RiVFMyU4MiU5MiVFMyU4MSU5RiVFMyU4MSU4RiVFMyU4MSU5NSVFMyU4MiU5MyVFOCVBNyVBMyVFMyU4MSU4NCVFMyU4MSVBNiVFMyU4MSVCRiVFMyU4MSU5RiZ0eHQtY29sb3I9JTIzM0EzQzNDJnR4dC1mb250PUhpcmFnaW5vJTIwU2FucyUyMFc2JnR4dC1zaXplPTU2JnR4dC1hbGlnbj1sZWZ0JTJDdG9wJnM9YjA3NzFmMjVlMTYwMjI2YWNhN2FlODc5NjQxOWJhYTA%26mark-x%3D142%26mark-y%3D100%26blend64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZoPTc2Jnc9OTcyJnR4dD0lNDBNaXlhNDImdHh0LWNvbG9yPSUyMzNBM0MzQyZ0eHQtZm9udD1IaXJhZ2lubyUyMFNhbnMlMjBXNiZ0eHQtc2l6ZT0zNiZ0eHQtYWxpZ249bGVmdCUyQ3RvcCZzPWE0OTgxZmUwYzUwN2I2MjcxZTRkZjdiOTVmNzE1Mzk2%26blend-x%3D142%26blend-y%3D504%26blend-mode%3Dnormal%26s%3D8e116d03b543e4c9c09932c886da7e98)