私はmatlab環境がないため、schemeで数値計算環境を作り始めた。一方、世間にはmatlab以外にもscilabやpythonのnumpy,scipy等 すぐれた数値計算ライブラリがある。それらは、すぐれている。 それらを使わずに、数値計算環境を自作するのは、とてつもなく非生…

今日は、べき乗法（power iteration）のプログラムをschemeで書く。

9.1節「対称行列の固有値」を読み、ヤコビ法を実装した。

今読んでる本、"数値計算: 洲之内著,石渡改訂,サイエンス社" 豊富な例題が並んでるのはありがたいけど、 数学の厳密性に関しては、かなり不満がある。

このエントリーは、完全に自分用の勉強ページです。

「（全てのサイズの）フィルベルト行列の逆行列は、どの成分も整数である。」 という定理がある。 Richardson "The Filbert Matrix" 1999https://arxiv.org/abs/math/9905079 これは、ヒルベルト行列 Choi "Tricks or Treats with the Hilbert Matrix" 1983…

他の方のブログを見て知った。 転置行列は、Schemeで短く書ける。

これから数値計算プログラムを書くに先立って(まだ書けてない)、 これからできるであろう自分のプログラムを「客観的」に評価する必要を感じている。

円周率の近似計算には、ラマヌジャンの公式を使うと良いそうなので、試しに自分でもコードを書いてみた。

数値計算を一から勉強をしようと思い、