私は、Schemeで代入を使うのが苦手である。 そこで、次のような練習問題を解いてみる。 問題1：10000以上で、２番目に小さい素数を求めよ。 問題2：10000以上で、ｎ番目に小さい素数を求めよ。

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

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

Schemeプログラマのレベル https://practical-scheme.net/wiliki/wiliki.cgi?Scheme%3AScheme%E3%83%97%E3%83%AD%E3%82%B0%E3%83%A9%E3%83%9E%E3%81%AE%E3%83%AC%E3%83%99%E3%83%AB10 とりあえず自己採点してみた。nibosiiwasi 自己採点結果 レベル2～3（線…

xyzzyのカスタマイズをしてみた。

私は最近、日常の事務の仕事で、 xyzzyとGauche(schemeの処理系)を使い、 最後にエクセルを少し使って書類を作っている。schemerに憧れ、仕事と趣味の両方で、無理やりxyzzyとGaucheを使うことにした。これは自分にとって良かったと思う。何よりschemeは楽し…

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

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

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

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

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

今日、仕事でSchemeを使ってたら、ときどき掛け算でおかしなことが起きました。

今日は、 近況。脱エクセルを目指して - nibosiiwasi’s blog について、詳しく書きます。

前回までのSchemeで数値計算。

なにも考えずに、Scheme で数値計算する主旨のブログを勢いで書き始めたことを、私は後悔し始めています。このままでは、明らかに非実用的なプログラムができる。 申し訳ありませんが、この計画、中断させて頂きます。もっと他のことも勉強してから、数値計…

今日は、逆行列を計算するプログラムを書いてみた。このページの後ろの方で、動作確認してみる。 まず逆行列を単純にガウスの消去法（＝掃き出し法）で計算するプログラム。ページ最後の追記をご参照下さい。

SICPのサブセクション3.3.1を自習する。

今日は、Schemeの変数周辺を試してみる。自分用の備忘録。

今日は、今後の方向性を決めるだけ。数値計算の中でも、固有値の計算は大きなテーマであり、大きな壁でもある。行列の固有値計算に関係する、次の３つの事実が広く知られている。

線形代数周辺の問題を解けるようにすることは、私の大きな目標の１つである。数値計算をする上で、線形代数、行列計算は避けて通れない。まず、行列計算の本当に基本的な部分をlist操作を使って書く。Gaucheにはarrayもあるのは知っているが、あえて使わない…

今もそうだが、学生のときもやはり、私は落ちこぼれだった。6〜8年前の学生時代、私はMATLABを使っていた。学校を卒業してMATLABを使えなくなった瞬間、私はPCで何もできなくなった。全然プログラミングできない自分が情けなかった。学生時代にMATLABで書い…