setup diary

_ diffinvを使って

Rで数値積分をする方法を聞かれたので、ついでにここにも書いておくことにする。数値積分をするためには、基本的にはy*dxを足していけば良い。dxを求めるには、diff(x)とすれば良いが、数が一つ減ってしまうので、そのままではyをかけられない。そこで、最初の要素を除いてy[-1]*diff(x)として、それを足し上げていくと数値積分となる。つまり、cumsum(y[-1]*diff(x))とすれば良い。

この方法では少し誤差が大きいので、せめて台形積分にしてみよう。yのところは、隣り合ったyの値の平均にするために、(head(y,-1)+tail(y,-1))/2として、cumsum((head(y,-1)+tail(y,-1))/2*diff(x))とすれば良い。この最初に0を加えるために、diffinv((head(y,-1)+tail(y,-1))/2*diff(x))とすると、xと数が同じになって良いだろう。

_ bullseye

そろそろdebianの次のversionがreleaseされるのではと調べてみたら，その翌日の8/14にリリースされることがわかった．そして，研究上で使っているMateriAppsのapt sourceがいつ公開されるかなと思って居たら，8/18に公開された．予想していたよりも早かった．新しいOSが出てから半年ぐらい待ってから更新することが多いのだが，今回は早めに更新してみた．DVDのISOをUSBに入れて，そこからbootしてインストールした．まずは英語で入れて，設定などを行った．desktopはlxdeにして，それ以外には以下のようなパッケージを入れた．

CPU: amd64-microcode
editor: emacs emacs-mozc
programming: ruby r-base r-cran-rgl r-cran-gsl r-cran-misc3d maxima
tex: texlive-lang-japanese texlive-latex-extra texlive-publishers
graphics: imagemagick gimp inkscape
mail: thunderbird-l10n-ja
tools: aptitude rsync
printer: cups
japanese: ibus-mozc

MateriAppsのsource listもhttps://sourceforge.net/projects/materiappslive/files/Debian/sources/から入れた．しかし，VESTAはまだ対応していないようなので，しばらく待たないといけない．あと，aptはcontribとnon-freeは追加しておいた． あとは，必要になったときに設定すれば良いかな．