はてなブックマーク

s=0d0 do ix=1,Nx cx=1d0 if(ix.eq.1.or.ix.eq.Nx)cx=0.5d0 x=(ix-1)*(xb-xa)/dble(Nx-1)+xa do iy=1,Ny cy=1d0 if(iy.eq.1.or.iy.eq.Ny)cy=0.5d0 y=(iy-1)*(yb-ya)/dble(Ny-1)+ya do iz=1,Nz cz=1d0 if(iz.eq.1.or.iz.eq.Nz)cz=0.5d0 z=(iz-1)*(zb-za)/dble(Nz-1)+za s=s+cx*cy*cz*f(x,y,z) enddo enddo enddo s = s*((xb-xa)/dble(Nx-1))*((yb-ya)/dble(Ny-1))*((zb-za)/dble(Nz-1)) 3,4桁以上の計算精度が欲しい時、台形則では間に合いません。 刻み幅を小さくすると精

正多角形とスピログラフの数式です。 正多角形(2次元) 正多角形(3次元) 星型多角形(2次元) 星型多角形(3次元) トロコイド(2次元) トロコイド(3次元) スピログラフ 2016/10/24 幾何学模様を作りたい方は幾何学模様(gnuplot)へどうぞ。 [adsense1] 正多角形の式(2次元) 正n角形の数式は[1]を参考に修正しました。 正n角形の数式は、 \( \begin{equation} \left\{ \begin{aligned} x&\displaystyle =\cos t \cdot \frac{\cos \frac{a_n}{2}}{\cos\left[a_n\left\{\frac{t}{a_n}-{\rm floor}\left(\frac{t}{a_n}\right)\right\}-\frac{a_n}{2}\right]}\\ y&\displ

ルンゲ=クッタ法(Runge-Kutta method、RK法)とは? 僕の知る限りの知識で紹介します。 特に良く使われる陽的ルンゲ=クッタ法は、 ・実装が簡単 ・良いアルゴリズムではない という手法です。 良いアルゴリズムである陰的ルンゲ=クッタ法は、 陰的ルンゲ=クッタ法 をご覧ください。 もくじ ルンゲ=クッタ法の系統 4次ルンゲ=クッタ法が使われる理由 陽的ルンゲ=クッタ法の段数と次数について 誤差について Butcher tableによるルンゲ=クッタ法の記述 高次の陽的ルンゲ=クッタ法(10,12,14次) 埋め込まれた陽的ルンゲ=クッタ法 ルンゲ=クッタ=フェールベルグ法による刻み幅の自動制御(適応刻み幅制御) ルンゲ=クッタ=フェールベルグ法による刻み幅の自動制御のプログラム(実数/複素数 で １/２ 階微分方程式を解くプログラム) 不連続な点を含む場合 刻み幅制御のベンチ