輻射輸送

■輻射輸送
最初は mathematica にやらせたけど,結局 C で書くはめになった.
久しぶりに数値計算用のプログラムを書いたら結構大変だった.
特に math.h を半年ぶりに include した.
今回やったのは輻射輸送方程式を解くのだけど,自己無撞着になるように解かなきゃいけなくて,変数の管理が面倒だった.構造体使えば良かった.
1次元平行平板問題はまだいけるけど,2次元とかはしんどいな.