概要 バネマスダンパ系の運動を数値計算により解く 環境 julia v(1.4.1) "DifferentialEquations" => v"6.13.0" "Plots" => v"0.29.9" 問題設定 質量m，ばね定数k，粘性cのバネマスダンパ系を考える． 運動方程式は，となる． この式を，以下のように式変形…

概要 前回の記事（以下）で常微分方程式df/dt = cos(t)の数値積分を行った． Julia(1.4.1)~常微分方程式の計算~ - Tagebuch von Spargel 本記事ではsolverによる精度の比較を行う． 結果 解析解との誤差の平均，最大誤差，実行時間は下表のようになる． 計算…

やること Julia言語で数値積分を行う． 数値積分に使うものとしてDifferensialEquations.jlとSundials.jlがある模様だが，DifferentialEquationsを使用する． 基本的には以下の公式サイトに乗っている． Home · DifferentialEquations.jl 環境 JupyterLab "j…

やること juliaを使用して数列作成してy=sin(x)のグラフを描く． 三角関数の作成でハマったのでメモ 環境 julia v1.1 jupyter lab 引数xの作成 主に以下の三種類． linspaceは使えない模様 # matlabライクな書き方,0.1刻みで出力 x = 0:0.1:2*pi #rangeを使…